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13,  ABSTRACT 

This  report  suomar Izes  vork  in  reliability  analysis  In  VLSI  CMOS  circuits,  particularly 
electromigration  in  metal  lines,  power  busses  and  signal  lines,  hot  carrier  effects  in 
devices,  and  voltage  drop  in  power  busses.  In  addition  to  process  parameters,  all  of 
these  effects  can  be  shown  to  depend  on  current  flow  in  the  circuit.  More  specifically, 
e 1 ectromlgrat ion  and  hot  carrier  Induced  degradation  are  both  long-term  effects  and  art 
related  to  the  average  current  flow  over  time  under  all  possible  input  signals  that  the 
design  experiences.  Maximum  voltage  drop,  on  the  other  hand,  depends  on  finding  the 
maximum  current  flow  in  the  power  busses,  caused  by  a  specific  few  Inputs.  Thus,  our 
emphasis  in  this  work  la  on  developing  fast  and  reliable  methods  for  estimating  average 
and  maximum  currents.  In  addition,  we  have  also  improved  on  our  bus  extractor  to 
extract  accurate  RC  models  of  power  busses  from  layout  Information  for  electromigration 
and  voltage  drop  estimation.  This  work  is  part  of  our  goal  of  developing  CAD  tools  for 
estimating  physical  reliability  effects  in  VLSI  deslgn.s,  and  to  automatically  modify  the 
design,  or  at  least  suggest  design  changes,  in  order  to  Improve  the  short  and  long  term 
reliability  of  th(»  design  under  realistic  operating  conditions. 
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FINAL  REPORT  EVALUATION 


Current  Estimation  For  Electromigratiori,  Hot  Carriers,  and  Voltage  Drop. 
Contract  F30602-92-C-0059 


Reliability  Analysis  of  VLSI  circuit  designs  is  necessary  to  reduce  the  time  from  concept  to 
delivery  without  compromising  reliability.  The  inclusion  of  reliability  analysis  into  the 
design  process  for  designs  with  thousands  of  gates  and  hundreds  of  inputs  poses  significant 
difficulty  because  of  the  computation  times  involved,  and  the  /ery  large  input  space. 
Deterministic  simulations  are  inadequate  to  simulate  the  effects  of  months  or  years  of 
circuit  use  when  the  number  of  inputs  is  large. 

Probabilistic  and  Statistical  techniques  offer  a  way  to  obtain  the  averaged  long  term  effects 
on  circuits  with  many  inputs  and  thousands  of  gates.  This  effort  investigated  statistical 
techniques  applicable  to  digital  CMOS  designs.  It  applied  probabilistic  methods  to  hot 
electron  degradation  and  incorporated  input  siew  rate  and  signal  correlations.  It 
incorporated  signal  correlation  effects  into  the  statistical  simulation  for  maximum  current 
estimarion. 

A  companion  effort.  "Prototype  Rule  Based  Reliability  Analysis  for  VLSI  Circuit  Designs" 
developed  layout  based  reliability  rules  for  hot  electron  degradation,  and  implemented  them 
into  a  prototype  reliability  analysis  tool.  This  approach  has  been  applied  to  circuit  designs 
with  more  than  1  million  transistors  in  about  1  CPU  hour.  The  rule  based  approach  is  not 
probabilistic,  but  is  a  worst  case  analysis  assuming  that  every  gate  switches  in  every  cycle. 

Areas  where  future  research  should  be  directed  include  applying  probabilistic  and 
statistical  methods  to  sequential  circuits,  reliability  analysis  at  the  gate  level  or  at  the 
hardware  description  language  so  that  reliability  becomes  and  integral  part  of  the  design 
process,  and  to  extend  the  area  of  reliability  issues  being  addressed  from  hot  electrons  and 
electromigration  to  other  failure  mechanisms;  long  and  short  term  reliability  problems. 

'MARTIN  J  WALTER 
Project  Engineer 
RI./ERDD 


1.  INTRODUCTION 


This  research  is  concerned  with  the  development  of  computer-aided  methodologies  and 
tools  for  assessing  and  improving  the  reliability  of  chip-level  VLSI  circuit  designs.  Many 
e.xisting  computer-aided  design  systems  tend  to  ignore  physical  reliability  constraints  during 
the  design  phase,  and  .rely  on  assessing  circuit  reliability  by  accelerated  bum-in  tests  after 
manufacture.  However,  with  the  increase  in  chip  complexity  and  the  decrease  in  feature  size, 
reliability  issues  can  no  longer  be  ignored  during  the  Jesign  process  and  need  to  be  addressed 
rigorously  throughout  the  design  cycle. 

Reliability  considerations  include  both  long-tenn  and  short-term  effects.  Long-term 
effects  are  related  to  the  physical  aging  of  integrated  circuit  elements  and  lines.  They  include 
electromigratiou  effects  in  the  lines,  and  hot-carrier  degradation  and  oxide  breakdown  in  the 
devices.  Short-term  :Tect.s  include  voltage-drop  and  noi.se  in  the  lines,  electrostatic  discharge 
through  input  protection  circuits,  system-noise  induced  latchup,  and  charge-particle  induced 
soft  errors. 

In  (his  task,  we  have  concentrated  our  efforts  on  three  reliability  issues;  namely,  elec- 
trumigration  detection  and  prevention,  voltage-drop  reduction  in  the  power  and  ground  busses, 
and  hot -carrier  induced  degradation  in  the  devices.  These  reliability  issues  can  be  shown  to 
btt  related  to  the  current  How  in  the  circuit  during  switching  activities.  Electromigration  is 
related  to  the  average  (as  weh  as  the  variance)  of  the  current  density  waveform  in  a  metal  line 
(averaged  over  all  possible  mputs).  Hot-electron  effects  arc  related  to  the  average  current 
flowing  through  the  devices,  while  maximum  tliagc  drop  is  related  to  the  maximum  current 
waveform  flowing  in  the  bus 
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Before  presenting  our  accomplishment  during  the  past  year,  we  give  a  brief  overview  of 
the  CAD  system  that  we  are  developing  to  implement  and  test  our  methodologies  and  algo¬ 
rithms.  Figure  1  shows  a  block  diagram  of  the  CAD  system.  Given  a  design  layout  specified 
in  CIF  (Caltech  Intermediate  Format),  which  could  be  specified  hierarchically,  program 
iCHARM  extracts  SPICE-file  from  it,  including  interconnect  parasitics,  and  sorts  out  the  CIF 
subfiles  of  the  power  and  ground  busses,  with  bus  contacts  labeled  to  match  the  corresponding 
node  connections  in  the  extracted  SPICE  file.  The  extracted  SPICE-file  has  the  same 
hierarchical  structure  as  the  original  CIF  file. 

The  SPICE  file  forms  the  input  to  the  simulators,  iPROBE-c,  iPROBE-d,  and  iMAX 
(The  SPICE  file  can  also  be  used  as  the  input  file  to  our  other  simulators  to  perform  timing 
and  fault  simulation  using  input  test  vector  sets;  these  other  'imulators  are  not  indicated  in  the 
Figure).  iPRQBE-c  is  a  probabilistic  simulator  which  computes  the  average  and  variance 
current  waveforms  at  declared  contact  points.  iPROBE-d  estimates  the  average  relative  dam¬ 
age  due  to  hot-carrier  effects  within  the  circuit  devices.  iMAX  computes  an  upp>er  bound 
current  waveform  envelope  at  contact  points  JET  is  a  bus  extractor  dedicated  to  extracting 
the  RC  mixlels  of  the  bus.  The  output.s  of  JET.  iMAX  and  iPROBF,  c  form  the  inputs  to  a 
bus  analyzer  and  optimizer  that  recommends  changes  to  the  bus  line  widths,  if  necessary  ,  in 
order  to  meet  both  electromigration  and  voltage-drop  constraints.  The  outputs  from 
iPROBE  c  and  iPROBE-d  can  be  used  as  pan  of  a  multiobjectivc  function  in  a  design  optimi¬ 
zation  program,  which  we  are  developing,  that  aims  at  optimizing  the  design  with  respect  to 
area,  timing,  power,  and  reliability  measures.  We  are  also  developing  a  program,  ACCORD, 
for  automatically  thangmg  the  design  at  the  logic  and  functional  levels  when  such  a  change  is 
recommended  in  the  o'^.omization  process. 
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In  Section  2,  we  present  the  accomplishments  achieved  in  iPROBE~c,  and  in  Section  3, 
the  accomplishments  in  iPROBE-d.  In  Section  4,  we  report  on  the  results  obtained  in  iMAX. 
In  Section  5,  we  describe  the  improvements  done  in  the  bus  extraction,  JET2.  Finally  in  Sec¬ 
tion  6,  we  give  a  summary  and  conclusions. 

Part  of  this  work  was  also  supported  by  the  Semiconductor  Research  Corporation  and  by 
Texas  Instruments,  Inc.  The  authors  would  like  to  thank  Professor  Farid  Najm  for  his  contri¬ 
butions  to  this  work;  in  particular  his  contribution  to  the  maximum  current  estimation 
approach  implemented  in  iMAX.  They  also  would  like  to  thank  Carolyn  Genzel  for  her  help 
in  preparing  this  report. 
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2.  ESTIMATION  OF  AVERAGE  CURRENT  WAVEFORMS 


2.1  Introduction 

Average  current  estimation  has  many  applications  in  VLSI  circuit  reliability  analysis  and 
design.  One  application  is  average  power  estimation.  Another  important  application  is  elec¬ 
tromigration  estimation.  Electromigration  in  metal  lines  has  been  shown  to  be  experimentally 
related  to  the  current  density  waveform  flowing  in  the  lines,  as  well  as  parameters  related  to 
the  metal  material  and  dimensions.  In  order  to  estimate  electromigration  effects  in  the  power 
busses  of  a  chip  design  under  expected  operating  conditions,  it  is  necessary  to  compute  the 
average,  as  well  as  the  variance,  of  the  current  density  vavefomi  in  each  section  of  the  bus 
under  all  possible  inputs  that  the  design  might  experience.  To  do  this,  one  has  to  first  com¬ 
pute  the  average  variance  current  waveforms  drawn  by  the  circuit  at  contact  points  to  the  bus. 
A  brute  force  approach  is  to  perform  exhaustive  simulation  on  the  design,  collect  all  the 
necessary  data,  and  then  compute  average  waveforms.  Such  an  approach  is  prohibitively 
expensive,  if  not  impossible,  since  the  input  space  is  usually  very  large.  Another  approach, 
which  is  sometimes  used  in  industry,  is  to  apply  Monte  Carlo  methods  to  high-level  functional 
description  of  the  design  to  obtain  toggle  rates  at  all  the  circuit  nodes,  which  in  turn,  can  be 
u.sed  to  obtain  average  waveform  of  each  subcircuit  or  gate.  However,  such  an  approach  usu¬ 
ally  docs  not  include  accurate  delay  and  timing  information  since  a  large  number  of  simula¬ 
tions  are  carried  out.  Our  approach  is  to  use  high-level  simulation  methods  to  obtain  signal 
statistics  at  the  output  of  flip-flops  in  sequential  designs.  Using  these  signal  statistics  at  inputs 
to  the  combinational  blocks  of  the  design,  we  apply  probabilistic  simulation  techniques  to 
obtain  average  and  variance  current  waveform  drawn  by  each  gate  or  subcircuil  including 
accurate  timing  information. 


During  the  past  year,  we  have  carried  out  a  number  of  improvements  on  our  probabilistic 
simulation  approach,  both  in  accuracy  and  in  speed.  The  improvements  are  carried  out  at  the 
subcircuit  or  gate  level,  where  the  statistics  of  the  current  and  voltage  waveforms  and  the 
delays  are  computed  more  accurately,  and  at  the  global  level,  where  signal  correlations  are 
computed  and  their  effects  on  the  subcircuit  analysis  are  taken  into  account.  In  the  following 
for  completeness,  we  review  our  approach.  We  then  explain  the  improvements  we  have  done 
and  include  simulation  results. 

2.2  Single  Gate  Probabilistic  Simulation 

We  consider  MOS  digital  circuits  that  have  been  partitioned  into  channel-connected  sub¬ 
circuits  or  gates.  We  assume  that  the  statistics  of  the  signals  at  inputs  to  each  gate  are 
specified  as  a  primary  input  or  computed  from  the  outputs  of  the  driving  gates  in  the  form  of 
probabilistic  waveforms.  Figure  2.1  shows  an  example  of  a  typical  probabilistic  waveform. 
The  statistics  of  the  cunent  waveform  drawn  by  a  CMOS  gate  after  switching  consist  of  the 
expected  value  and  the  variance  waveforms  where  the  expected  value  waveform  £[i(f)] 
represents  the  average  waveform  over  all  possible  power  supply  waveforms  drawn  by  the  sub- 
circuit,  and  the  variance  waveform  represents  the  variance  of  the  distribution  of  the  currents. 
The  variance  wavefomi  is  calculated  as  V(i(f))  =  £[(i(t)  -  £'[t  (?  )])’)  at  every  time  point  t 

The  expected  value  and  variance  waveforms  are  assumed  to  be  triangular  in  shape  luid 
each  is  specified  by  a  peak  value  and  a  time  span.  To  compute  the  peak  value  and  the  time 
span  of  the  expected  value  and  variance  current  pulses  each  gate  is  reduced  to  an  equivalent 
inverter  shown  in  Fig.  2.2.  The  p-pan  and  the  n-pan  of  the  gate  arc  rei.luced  to  one 
equivalent  edge  each  by  a  series-parallel  reduction  as  follows:  Fach  transistor  m  a  gate  has  a 
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conductance  gg„  when  it  is  on  and  zero  when  it  is  ojf.  T  -  probability  of  an  nMOS  transistor 
being  on  is  the  probability  of  its  gate  node  being  "high,"  while  the  probability  of  a  pMOS 
transistor  being  on  is  the  probability  of  its  gate  being  "low."  Thus  the  conductance  of  the 
transistor  c;in  be  viewed  as  a  random  variable.  The  expected  value  £[g]  and  the  variance 
V(^)  of  th<;  conductance  are  defined  respectively  at  any  time  t  ether  than  the  time  points 
where  switching  occurs  as: 

E\g]=^gon  (2.1) 

Vig)  =  gl  X  P^it)  -  (g,,  X  F  ))2  (2.2) 

where  P/,(r)  the  probability  that  the  edge  will  conduct. 

The  aim  of  the  reduction  process  is  to  find  an  equivalent  edge  between  the  output  node 
and  the  ground  node  (V^  node  for  the  p-part)  so  that  the  current  drawn  at  the  ground  contact 
{Vdd  contact  for  the  p-part)  can  be  calculated.  It  is  further  assumed  that  the  circuit  can  be 
reduced  in  a  series-parallel  fashion,  which  is  the  most  prevalent  case.  Thus,  at  every  step  of 


Figure  2. 1  Example  of  a  probabilistic  waveform. 
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Figure  2.2  Model  of  a  CMOS  gate. 

the  reduction  procedure  two  edges  are  combined  into  one.  The  analysis  presented  in  the  fol¬ 
lowing  sections  will  be  for  the  n-part  only,  since  the  p-part  is  its  complement.  The  following 
notation  will  be  used  in  the  sequel: 

X,  y:  the  conductances  of  the  two  edges  (transistors)  to  be  reduced  to  one. 

g:  the  conductance  of  the  resulting  equivalent  edge. 

:  probability  that  x  is  conducting, 

Pif^^ :  probability  that  x  is  switching  from  off  to  on. 

•  Parallel  reduction:  Assuming  independent  inputs,  the  parallel  case  is  straightforward. 
The  equivalent  conductance  is  given  by: 

g  =  x  +  y  (2.3) 

where  x,  y  and  g  are  random  variables. 

The  stochastic  quantities  describing  the  equivalent  edge,  g,  are  given  by  the  following 


9 


expressions: 


^ hp  ~  ^hx  +  ^ hy  ■  ^ hy 

(2.4) 

^Ihp  -  P  Ihx  ^ly  thy  ^  lx  ~  ^Itxx^  Ihy 

(2.5) 

E[g]^E[x]^E\y] 

(2.0) 

V(g)-  V(.T)-H/(y) 

(2.7) 

where  the  subscript  p  (for  parallel)  denotes  those  quantities  associated  with  the 
equivalent  edge.  These  lesults  follow  from  simple  probability  theory  and  are  not  approxi¬ 
mations  but  exact  representations  of  the  stochastic  parameters  of  the  equivalent  edge. 


Series  reduction:  For  the  series  case,  the  equivalent  edge,  g ,  is  high  or  is  switching  from 
low  to  high  are  calculated  by  the  following  expressions. 


Phs  =  f^hx^hy 

(2.8) 

Ihx^hy  ^Ihy^kx  ~  ^ Ihx  ^ Ihy 

(2.9) 

the  quantities  associated  with  the  equivalent  edge  are  here  denoted  by  the  subscript  5 . 


In  computing  the  expected  value  and  the  variance  of  the  series  combination,  we  follow  a 
different  approach  from  the  one  proposed  in  [I]  and  [2],  w'hich  leads  to  more  accurate  results. 
We  consider  the  equivalent  conductance  to  be  expressed  as: 


8  ix,y)= 

x+y 


(2.10) 


Following  [3],  a  Taylor  series  expansion  of  (2.10)  gives: 


^ L?  1  =  g (-«  J )  +  ^( cfi'  +  2  ~-|-r  o,  a  +  aj 

Bx ''  dv 


Bx 


a}  =  +  2(^)(|«-)r  ,,  o,  +  (-?->  a’ 


By^ 
2 

”“) 

h  ^ 


(2.11) 

(2.12) 
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where  o*  i»  »hc  vanance  V'd  i  and  r  is  Jhe  coirclation  ccwjfttieni  between  i  and  i  i/en>  d 
wc  assume  independence >  The  nghi  hand  side  is  cvaliiuied  ai  if;  j;  l.f  fs  5)  Applying  *2  1  i  > 
and  (2,12)  lo  (2.!()),  yieids 


£1,1 .  £UI£b-l  _  gbl’vu)  . 


E[x]  +  E[y]  (/-(ji)  +  (E{x]  +  Ely])^ 

=  1^4rTrv(^) 


(2  13) 


(2.14) 


E[x]-¥E[y]‘  '  ‘E[j:]  +  £Lyl 

We  have  found,  however,  that  (2.13)  and  (2,14)  do  not  give  accurate  results  espi?;cially  when 

P;,  is  rather  small.  By  using  conditional  probabilities  we  obtain  very  accurate  results.  In  the 
conditional  probability  approach,  the  expected  value  and  the  variance  of  the  equivalent  edge 
are  calculated  given  that  the  equivalent  edge  g  conducts.  With  this  condition  (2.13)  and  (2.14) 
become: 


Ff  I  ?i0  V  ?t01  =  gU  l-c  yiQ^E[y  ly  ?t0]  _  E[y  \y  ^Q]‘^V{x  Ijc  ?i0) 

+  £[3;ly^i01  (E[jc  ijc  ^tOl  +  £[y  ly  ^0])^ 


£[.rU?^0]^V(yly  ?t0) 
(£:[xLx:  ^0]  +  E[y  ly  ?^0])-'’ 


V(g  Ia  9i0,y  9t0)  =  1 


_ £  [yly  ^0] _ 

£  [x  !x  ^  0]  +  £  [y  ly  0] 


J'*V(x  lx  ^0) 


+  [ 


_ £[x  lx  ^0] _ 

£[x  lx  ;e0]  +  £[y  iy  7*0] 


]V(yly  ;t0) 


(2.15) 


(2.16) 


where 


£[.dx7i0j  = 

£[x) 

Pfa 

(2.17) 

£[x^lx  55:())  = 

E{x^] 

Ph. 

(2.18) 

o 

!l 

t>3 

to 

0]  -  £'[x  l.v  -/  Ol 

(2.19) 

E[gl  and  V(g)  are  calculated  from  (2,15)  and  (2.16),  using  equations  similar  to  (2. 17)-(2. 19). 
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(2  and  <2  16)  ha%v  been  letted  on  a  large  number  of  CMOS  subi:'rcuits 
vsiih  the  ntimher  <>1  franst*ior'>  ranging  in  (he  Nubcsfcud  l>eivveen  4  and  ?3,  Jirrangcd  in  various 
■veixs-parallci  comhinasions  The  resiiUs  are  compared  with  those  of  (Ij,  as  well  as  with 
exfiaustivc  SPICE-like  simulations.  For  comparison,  the  equations  derived  in  [1]  for  series 
combination  are: 

E[x]P^y  ^  Ey]P^ 

li&L  ^  IM.  ^  ,221) 

Klg]‘'  fUl**  E[yf 

It  is  observed  that  in  calculating  the  expected  vaiue,  the  proposed  method  has  3%  max¬ 
imum  error  compared  to  exhaustive  simulation,  while  the  previous  method  has  50%  maximum 
error;  while  in  calculating  the  variance,  the  new  method  has  a  12%  maximum  error  while  the 
old  one  92%. 

Comparing  our  new  results  with  results  obtained  using  the  original  approach,  we  con¬ 
cluded  that: 


•  the  new  series  reduction  method  produces  consistently  accurate  results, 

•  the  error  in  the  reduction  process  seems  to  be  only  slightly  dependent  on  the  input  proba¬ 
bilities,  with  higher  error  occurring  when  the  input  probabilities  are  close  to  zero  or  one, 
a  result  corroborated  by  all  the  results  to  date,  and 

•  the  nev/  method  con.sistently  outperforms  the  previous  method. 
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23  Elimination  During  Switching 


T*(e  basic  circuit  model  of  a  CMOS  gate  is  shown  in  F>g.  2.2.  CMOS  circuits  draw 
current  from  the  bus  only  when  their  outputs  change  states.  By  the  nature  of  probabilistic 
simulation,  a  probabilistic  event  at  the  input  always  produces  a  probabilistic  event  at  the  out¬ 
put.  Iherefore,  the  calculation  of  a  current  pulse  must  be  done  whenever  an  input  is  switch¬ 
ing.  We  first  assume  that  the  inputs  to  a  subcircuit  are  independent  and  switch  instantane¬ 
ously.  The  effects  of  signal  dependence  on  the  elimination  process  is  considered  in  Section 
2.5  below.  Since  the  signals  are  assumed  to  switch  instantaneously,  we  make  the  following 
observation: 

•  Observation:  Since  the  probability  that  two  input  signals  to  a  gate  switch  at  the  same 
time  is  zero,  if  the  edge  that  switches  is  in  parallel  with  an  edge  that  is  conducting,  then 
switching  will  have  no  effect  on  the  output  since  the  edge  that  is  ON  effectively  shorts 
the  edge  that  switches;  as  a  result,  the  switching  transistor  will  not  cause  the  gate  output 
to  switch.  It  follows  then  that  for  the  output  to  switch,  one  of  the  edges  must  switch  and 
all  the  conducting  paths  between  the  output  node  and  the  power  node  mu\t  go  through 
the  edge  that  has  switched, 

Ba.scd  on  the  above  observation,  the  calculation  of  the  equivalent  conductance  at  switch¬ 
ing  time  is  simplified.  Suppose  that  edge  k  is  switching;  then  for  all  edges  i  that  fonn 
conducting  paths  with  edge  k ,  it  follows  from  the  independence  assumption  that: 

Pii  is  ON  \k  is  switching  =  P{i  is  ON)  ~  Pf^j  (2.22) 

Since  at  switching  time  all  the  current  rnu.st  pass  through  the  transistor  that  is  switching,  then 
all  paths  that  bypass  that  transistor  should  not  conduct,  as  shown  in  Fig.  2.3,  which  leads  to 
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the  following  graph  reduction  algorithm: 


1.  Assign  probability  1  to  the  edge  k  that  actually  switches,  before  reduction  begins. 

2.  Separate  the  edges  into  switching  and  non-switching  edges.  A  switching  edge  is  defined 
as  the  one  that  actually  switches  and  every  equivalent  edge  that  contains  it  during  the 
process  of  graph  reduction.  Otherwise  an  edge  is  non-switching. 

3.  Proceed  to  graph  elimination  as  described  in  Section  2.2. 

(a)  Series  reduction  is  performed  according  to  (2.15)  and  (2.16). 

(b)  In  parallel  reduction  we  distinguish  between  two  cases: 

•  If  two  non-switching  edges  are  in  parallel  then  elimination 
proceeds  as  usual. 

•  If  one  of  the  edges,  which  could  be  an  equivalent  edge,  is 
switching  then  the  edge  in  parallel  with  it  is  assigned 
P),  =  0  and  the  probability  that  this  edge  is  not  conducting 
is  stored.  Let  a  be  the  set  of  these  parallel  edges. 

4.  Calculate  ]  and  V'(4)  for  the  peak  current  4  flowing  through  the  switching  edge  k, 
which  is  also  the  total  peak  current  drawn  by  the  gate,  the  following  expressions  have 
been  derived: 


t  €  a 

=  Kii  X  EIO,L,I  X  Em  X  n  E„ 

I  6  a 


(2.23) 

(2.24) 


where  is  the  conductance  of  the  equivalent  edge  between  the  output  node  and  the 

power  node. 
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It  is  evident  that,  because  of  independence,  at  switching  time  there  is  only  one  switching 
edge  in  the  graph  that  is  being  reduced  and  that  the  only  edge  remaining  after  elimination  is  a 
switching  edge. 

The  time  spans  of  the  current  waveforms  are  computed  from  the  basic  gate  model  shown 
in  Fig.  2.2  following  the  method  proposed  in  [1],  [2].  Let  E[q]  be  the  expected  value  of  the 
charge  stored  in  the  output  capacitances  C„  and  : 

Cp  (2.25) 

where  P(^  and  Pf^i^  are  the  probabilities  that  the  output  will  go  from  low  to  high  and  vice 
versa.  Also 


...  . 


EUn] 


(2.26) 


where  E[iq]  is  the  expected  value  of  the  charge-current  product  and  £[/),],  E[ip]  are  the 
expected  vuues  of  the  current  drawn  by  the  n-part  and  the  p-part  respectively.  The  time 
spans  (Tf;)  for  the  expected  value  and  (Xv)  for  the  variance  of  the  current  waveforms  are 
given  by  the  following  relations. 


Edges  that  are 
assign^ 
prob^d)ili^  zero 
during  elimination 


Figure  2.3  Elimination  example  when  one  edge  switches. 


‘  £(<! 

T  1  X  (  ^(*^1  ~  gt<  [^1  ■ 
^3  V(i) 

where  /  =  i„  +  i-,  £[j]  its  expected  value  and  V(i)  its  variance. 


(2.27) 

(2.2S) 


Example  1:  This  example  illustrates  the  accuracy  of  calculating  the  expected  value  and 
variance  current  waveforms  using  the  new  method.  Figure  2.4(a)  is  a  two-input  NAiS’D  gate, 
with  two  inputs  A  and  B.  Figures  2.4(b-d)  show  the  expected  value  of  the  cunent  waveforms 
using  the  new  approach  for  (Fig.  2.4(b)),  where  and  tf,  are  the  transition  times  for 

inputs  A  and  B  respectively,  -Q4ns  (Fig.  2.4(c))  and  tf,  -  =  4ns  (Fig.  2.4(d)). 

compared  each  time  with  exhaustive  SPICE  runs.  The  proposed  method  successfully  tracks 
SPICE  results  in  all  cases. 


2.4  Delay  Calculation 

After  having  calculated  the  peak  and  the  time  span  of  the  current  drawn  by  the  gate,  we 
need  to  calculate  the  output  voltage  waveform  that  drives  the  fanout  gates.  This  includes  the 
determination  of  the  probability  that  the  output  node  is  high  after  the  transition,  the  probabil¬ 
ity  that  it  will  transition  from  low  to  high,  and  the  time  when  the  output  transition  will  take 
place  (i.e.,  the  delay  of  the  output  event). 

The  stochastic  quantities  (steady-state  and  transition  probabilities)  can  be  derived  from 
the  following  expressions: 


^  1  h  .out  ~  .p  ^  ^  hi  ,k 

(2,29) 

1  .out  ~  ji  ^  ^  th  .k 

(2.30) 

^h.oui^^  )  ~  ^Ih.out  ~  ^hl.out 

(2.31) 
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Figure  2.4  Simulation  of  transitions  within  a  single  interval, 
where  Pih,oui  the  probability  that  the  output  will  go  from  low  to  high,  P^i  is  the  proba¬ 
bility  that  the  output  will  go  from  high  to  low.  is  the  probability  that  the  output  is 

high  before  the  transition,  is  the  probability  that  the  output  is  high  after  the  transi¬ 

tion,  Pf^  p  is  the  probability  that  the  p-part  of  the  gate  conducts  (as  calculated  by  the  graph 
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is  the 


elimination  process).  „  is  the  probability  that  the  n-part  of  the  gale  conducts,  ^ 
probability  that  the  input  node  that  is  causing  this  transition  is  going  from  low  to  high,  and 
f.  the  probability  that  this  input  node  is  going  from  high  to  low. 

For  the  calculation  of  the  delay  of  the  output  event  a  new  method  is  presented  based  on 
the  deterministic  gate  delay  model,  which  applies  to  both  switching  from  low  to  high  and  high 
to  low,  namely: 

rrf=/n2x--^~^  (2.32) 

"  G 

where  is  the  gate  delay.  Expression  (2.32)  is  derived  from  a  simple  RC  circuit  driven  by  a 
step  voltage  input.  Using  simple  probability  theory,  we  arrive  at  the  following  expression  for 
the  delay  of  each  of  the  transitions. 

For  the  transition  from  low  to  high 


\  =  ln2x{C  +  Q)  X  ( 


E\Gr 


and  for  the  transition  from  high  to  low 


^  [tjn  1  =  /n  2  X  (C  +  C„  )  X  ( 


1 


EiG, 


V(G^) 

iFAG,])^ 


(2.33) 


(2.34) 


where  G„  and  are,  respectively,  the  equivalent  conductances  for  the  n-part  and  the  p-part 
when  the  output  is  switching.  For  the  total  average  delay,  the  following  expression  is  used, 
which  is  the  weighted  average  of  the  delays  from  low  tc  high  and  from  high  to  low. 


E[‘dn]^  Pou,M  +  ]  X  P, 


oui  .h  I 


+  P 


ip 

out  .1  h 


oul  ,!  h 


(2.3.*;) 
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Example  2:  This  example  illustrates  the  accuracy  of  Eq.  (2.35)  compared  to  [1]  and  to 


exhaustive  SPICE  simulation.  Figure  2.5(b)  shows  the  graph  representation  of  the  n-part  of 
the  CMOS  gate  shown  in  Fig.  2.5(a).  All  inputs  are  in  steady-state,  except  input  A  to  which 
the  probabilistic  waveform  shown  in  Fig.  2.5(c)  is  applied.  There  are  two  numbers  associated 
with  each  edge  of  the  graph.  The  top  one  is  the  probability  that  the  gate  node  of  the 
corresponding  transistor  is  high  and  the  bottom  one  represents  its  on  conductance  in  m  C'  ■ 
The  total  load  capacitance  =  Cp  +  used  is  lOpF . 

The  average  delay  calculated  by  Eq.  (2.35)  is  found  to  be  3.64ns  compared  to  3.60ns  by 
exhaustive  SPICE  simulation,  while  the  method  used  in  [1]  yields  0.66ns ,  which  shov/s  that 
our  result  is  much  closer  to  SPICE.  The  improvement  in  accuracy  is  the  result  of  the  improve¬ 
ment  in  the  probabilistic  graph  elimination  algorithm,  the  elimination  during  switching  algo¬ 
rithm,  and  the  delay  calculation  from  Eqs.  (2.33)-(2.35). 


Figure  2.5  (a)  A  CMOS  gate;  (b)  Graph  of  the  n-part;  (c)  Input  waveform  at  node  A 
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2.5  Output  Event  Merging 


So  far,  in  computing  the  cunent  waveform  drawn  by  a  gate,  it  has  been  assumed  that 
every  event  at  an  input  to  the  gate  produces  one  distinct  event  at  its  output.  However,  when 
the  focus  is  shifted  to  the  output  voltage,  depending  on  the  delay  of  the  gate,  and  the  time 
difference  between  the  input  events,  two  input  events  may  produce: 

i.  One  output  event,  as  in  Fig.  2.4(b); 

ii.  one  output  event  with  a  small  glitch  that  does  not  affect  the  operation  of  any  of  the 
fanout  gates,  as  Fig.  2.4(c);  and, 

iii.  two  distinct  events,  as  in  Fig.  2.4(d). 

After  performing  a  number  of  detailed  circuit  level  simulations  on  typical  ga*e  structures,  we 
have  deterrmned  that  if  two  output  events  occur  within  a  time  interval  which  is  less  than  the 
propagation  delay,  of  the  gate,  then  the  two  events  should  be  merged  into  one  according 
to  Eqs.  (2.36)-(2.39)  below.  If  more  than  two  events  occur,  merging  is  done  by  taking  two 
events  at  a  time,  starting  with  the  first  two  events  in  the  queue. 

If  eventl  and  event2  are  to  be  merged,  then  the  corresponding  probabilities  of  the 
equivalent  event,  which  will  be  referred  to  as  event  e,  are  given  by: 

P  =  P  xh(r)  (2.36) 

=  (2.37) 

Pcih  ~  Plh\P hhZ  +  Pll\P lh2  (2,38) 

where  the  events  arc  assumed  indeciendent,  is  the  probability  of  signal  staying 

high  before  and  after  switching  and  P ^  \  ~  is  its  probability  staying  low  before  and 

after  switching,  .Since  f  is  the  combination  of  the  two  events,  the  new  event  is  scheduled  for 
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time 


I  -  +  ^hlO  ^  h  •*'  (^//i2  ^hl2^  ^  h 

^  ^  lh\  ^hix  P  Ihl'^  ^hll 

which  is  the  weighted  average  of  the  times  of  occurrence  of  the  two  events,  a  result  suggested 
by  simulation  results. 

The  NAND  gate  example  depicted  in  Fig.  2.6  further  illustrates  the  merging  process.  In 
this  example  we  have  a  two-input  NAND  gate,  which  is  driven  by  three  events  on  its  inputs, 
two  of  which  are  on  the  same  input  node.  Assume  that  events  1  and  3  are  closer  than  the 
defined  by  the  driver  gate  to  the  NAND  gate.  Thus,  the  two  events  are  merged  into  the 
equivalent  event  la.  Events  la  and  2  produce  two  distinct  events  at  the  output,  la’  and  2’ 
respectively.  But  again  these  two  events  are  closer  than  for  the  NAND  gate  and  are 
merged  into  the  equivalent  event  4. 

Gate  output  event  merging  is  an  integral  part  of  our  probabilistic  simulation  approach. 
Its  aim  is  to  suppress  small  glitches  at  the  output  of  the  gate  and  to  identify  large  glitches  that 


Figure  2.6  Example  of  merging  at  the  input  and  output  nodes  of  a  gate. 
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may  affect  fanout  gates.  It  also  dramatically  reduces  the  number  of  events  that  need  to  be 
simulated,  thus  reducing  the  overall  simulation  time  significantly  as  shown  in  Table  2.1.  In 
generating  the  results  in  Table  2.1,  signal  correlation  has  been  ignored.  The  effects  of  signal 
correlation  on  the  results  are  considered  in  the  next  .section. 


I  Table  2.1:  Simulation  Results 

Primary 

inputs 

Time  with  no 
mergingt 

Time  with 
mergingt 

Memory  usage 
(Kbytes) 

c432 

37 

388.1 

1.5 

456 

c880 

61 

31.2 

2.1 

680 

C1355 

42 

7S56.2 

4.1 

804 

cl  908 

34 

983.4 

7.2 

1128 

c2670 

234 

978.0 

7.6 

1512 

c3540 

51 

59201.9 

20.9 

c5315 

179 

4126.5 

19.4 

c6288 

33 

>24h 

95.6 

3796 

c7552 

208 

1809.9 

32.2 

3620 

t  In 

CPU  seconds  on  an  HP  90tX)/730 

2.6  Signal  Correlation 

So  far  we  have  assumeu  that  all  inputs  to  a  subcircuit  are  independent.  In  this  section 
we  consider  the  effects  of  signal  correlation.  Correlation  is  a  measure  of  the  interdependency 
of  two  or  more  random  variables  Correlation  can  be  introduced  between  two  or  more  node 
variables  due  to  the  presence  of  reconvergent  fanout  and/or  feedback  in  a  circuit,  as  in  F'g. 
2.7. 


There  are  two  issues  that  have  to  be  considered,  One  is  the  computation  of  signal  corre¬ 
lations  between  the  circuit  nodes,  and  the  second  is  the  calculation  of  the  effects  of  signal 
correlations  at  the  inputs  of  a  gate  on  its  probabilistic  analysis.  To  deal  with  the  first  issue. 


namely  to  obtain  signal  probabilities  and  correlation  coefficients  between  the  nodes,  we  use 
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Figure  2.7  A  16-transistor  implementation  of  an  XOR  gate, 
the  approach  described  in  [5],  with  some  modifications  needed  for  probabilistic  simulation. 
The  approach  in  [5]  estimates  the  node  probabilities  based  on  an  estimate  of  the  first  order 
correlations. 

The  second  issue,  namely  the  analysis  of  a  gate,  given  the  signal  probabilities  and  signal 
correlations  at  its  inputs,  is  explained  in  detail  in  [7]. 

Example  3:  Referring  again  to  the  i6-transistor  XOR  gate  shown  in  Fig.  2.7,  this  circuit 
is  selected  becau.se  of  the  very  strong  reconvergence  observed  at  the  inputs  of  three  of  the 
four  total  gates  comprising  this  example  (namely  gates  2,  3,  and  4).  The  total  average  current 


as  estimated  by  our  approach  with  and  without  correlation  considerations  is  shown  in  Fig.  2.8 
along  with  the  results  of  exhaustive  simulation  as  reference.  We  observe  that  the  waveform 
produced  without  coixelation  considerations  demonstrates  high  error  at  time  i  .4n^ ,  This  is  due 
to  the  fact  that  the  inputs  to  gates  2  and  3  are  highly  correlated,  and,  thus,  the  current  drawn 
by  those  two  gates  is  overestimated.  .Again  the  same  w-avefortAi  shows  significant  error  at 
time  point  1.35ni'  where  the  current  drawn  by  gate  4  (who.se  inputs  are  also  correlated)  is 
underestunafed  From  the  same  figure  we  can  observe  that  the  implementation  with  correlation 
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accurately  follows  the  exhaustive  current  until  time,  point  LSSns.  The  error  that  arises  there  is 
due  to  the  inability  of  the  correlation  coefficient  method  [5]  to  accurately  calculate  the  correla¬ 
tion  coefficient  between  the  input  nodes  of  gate  4.  Another  measure  of  the  accuracy  of  the 
methods  proposed  in  this  paper,  which  is  directly  related  to  the  applications  of  probabilistic 
simulation  (i.e.,  average  power  estimation  and  electromigration  degradation)  is  to  calculate  the 
integral  of  the  expected  current.  For  this  example  the  following  results  are  obtained  (5V  sup¬ 
ply  voltage): 

•  Exhaustive:  1.335  mW 

•  Without  correlation:  1.419  mW  (6.0  error) 

®  With  correlation:  1.387  mW  (3.9  enor) 

As  we  can  see  the  errors  a^e  small,  even  though  the  solution  with  correlation  produces  less 
error. 

Example  4:  In  this  example  we  consider  the  smallest  of  the  ISCAS85  benchmark  cir¬ 
cuits.  c432  [6].  To  be  able  to  perform  exhaustive  simulation  on  the  circuit  within  reasonable 
time  all  the  primary  inputs  are  assigned  to  logical  "1"  except  two.  Figure  2.9  shows  the  results 
from  exhaustive  simulation  and  from  probabilistic  simulation  with  and  without  correlation 
considerations.  As  we  can  see  the  estimated  waveforms  generally  foliow  the  exhaustive  one 
and  the  integnds  under  the  curves  are  about  the  same  (less  than  10%  error  for  both  implemen¬ 
tations  with  respect  to  exhau.stive). 

However,  if  we  look  at  the  current  drawn  by  an  individual  gate,  within  the  design,  that  is 
close  to  the  output  of  the  circuit  whose  input  signals  are  correhued,  it  c.an  be  observed  that  the 
waveforms  are  different  and  the  integrals  differ  by  a  factor  of  2  when  correlations  are  ignored. 
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while  the  error  with  correlation  considerations  is  only  25%, 

When  the  ISCAS85  benchmark  circuit  c432  is  analyzed  with  all  the  primary  inputs 
switching,  the  two  waveforms  produced  with  and  without  correlation  considerations  are  margi¬ 
nally  different,  as  shown  in  Fig.  2.9.  On  all  the  circuits  we  tested,  the  effect  of  correlation 
on  the  total  current  is  small,  a  result  consistent  with  previous  published  work  [1],  while  the 
effect  on  individual  gate  currents  teuds  to  be  significant. 

Table  2.2  shows  the  impact  of  correlation  calculations  on  the  overall  speed  of  the 
approach  for  the  ISCAS85  benchmark  circuits.  The  inclusion  of  correlation  estimation 
increases  the  computation  time  up  to  threefold  for  the  examples  shown  here,  and  could  be 
even  greater  for  larger  examples,  but  the  speed  remains  quite  reasonable,  and  certainly  better 
than  doing  an  exhaustive  analysis.  Thus,  there  is  a  tradeoff  between  accuracy  and  speed. 


Table  2.2:  Simulation  Results 

Time  with  no 
mergingt 

Time  with 
mergingt 

Time  with 
correlationt 

c432 

388.1 

1.5 

1.6 

2.1 

2.8 

cl  355 

7856.2 

5.0 

ci908 

983.4 

9.0 

c2670 

978.0 

14.7 

c3540 

59201.9 

20.9 

40.3 

c53l5 

4126.5 

19.4 

46.2 

c6288 

>24h 

95.6 

110.0 

c7552 

1809.9 

32.2 

140.2 

t  In  CPU  seconds  on  an  HP  90 

00/730 

In  our  program  two  methods  of  calculating  the  expected  current  are  implemented.  The 
first  is  described  in  Sections  2,2  through  2.5  where  correlation  is  ignored.  It  is  fast,  with  ,s.mall 
memory  requirements  and,  as  our  examples  up  to  now  show,  a«';te  accurate.  It  can  be  used  for 
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exhaustive 
without  correlation 


a  quick  yet  reliable  estimate  of  the  expected  cunent  waveform  of  a  given  circuit.  The  exten¬ 
sion  to  this  method  described  in  Section  2.6  takes  into  account  second  order  correlation,  that 
is  correlation  between  two  nodes  at  a  time.  This  is  by  construction  more  accurate  than  the 
previous  one,  when  considering  individual  gates,  since  it  takes  into  account  correlation,  but  is 
slower  (up  to  three  times  slower  in  the  examples  that  were  presented  here)  and  requires  more 
memory  space.  However,  it  does  not  seem  to  provide  enhanced  accuracy  when  the  total 
current  is  being  calculated.  Thus,  there  is  a  tradeoff  between  accuracy  and  speed  that  has  to 
be  taken  into  account  when  using  the  new  approach. 

Possible  explanations  to  the  minimal  impact  of  correlation  on  the  total  current  is  that  the 
uncorrelated  gate  inputs  in  the  studied  circuits  are  more  than  the  correlated  ones  and  mask 
their  effect  and/or  that  ignoring  correlation  introduces  random  error.  This  error,  when  the 
estimated  currents  are  superimposed  (added),  tends  to  zero.  But,  when  the  current  drawn  by 
each  individual  gate  is  estimated,  as  some  applications  require,  large  discrepancies  between 
estimated  and  actual  current  could  occur.  In  thv^se  applications  inclusion  of  correlation  in  the 
computations  is  imperative. 


2» 


3.  CURRE  VT  ESTIMATION  F^R  HOT-CARRIER  ETFECTS 


3.1  Introduction 

In  Chapter  2,  we  explained  the  application  of  probabilistic  simulation  to  estimate  the 
average  and  variance  curreitt  waveforms  drawn  by  the  gates  at  contact  points  to  the  power 
bus.  In  this  chapter,  we  report  on  the  work  we  have  done  on  the  application  of  probabilistic 
simulation  for  estimating  the  average  current  waveforms  flowing  in  the  devices  within  the 
gates  of  a  given  design.  These  currents  are  then  used  to  estimate  the  expected  hot-carrier 
(HCE)  induced  degradation  in  each  transistor.  We  then  propose  some  design  strategies  to 
reduce  such  degradation  and  its  impact  on  long-term  circuit  performance. 

Existing  work  on  HCE  mainly  focuses  on  developing  degradation  models  of  MGS 
transistors.  These  models  are  then  linked  to  general  purpose  circuit  simulators  such  as  SPICE 
to  simulate  circuit  performance  over  time  [8J,  [9],  [10],  [11].  This  type  of  approach  is  applica¬ 
ble  when  the  circuit  under  simulation  contains  only  a  few  transistors,  and  is  used  to  verify  the 
accuracy  of  the  degradation  models  For  VLSI  circuits  with  hundreds  of  thousands  of  transis¬ 
tors,  however,  the  use  of  a  circuit  simulator  such  as  SPICE,  even  with  a  simple  (level  one) 
transistor  model,  is  not  feasible,  especially  since  HCE  is  a  cumulative  effect,  and  the  com¬ 
bined  effects  of  all  possible  inputs  have  to  be  evaluated. 

Recently  timing  simulation  has  been  considered  as  a  tool  for  HCE  estimation  [12],  [1,3]. 
The  approach  in  [12]  employs  fast-timing  simulation  to  simulate  s  jbcircuits  with  no  damage, 
and  more  detailed  circuit  simulation  to  estimate  HCE  on  the  dama';;ed  subcircuits.  In  [13],  a 
macromodel  of  channel-connected  MOS  blocks  is  constructed  anc  fast  timing  simulation  is 
used  to  predict  the  degradation  effects  over  time.  This  approach,  however,  estimates  the 
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degradation  of  the  simplified  macromodels  instead  of  each  MOS  transistor  inside  the  original 
circuit.  In  both  these  approaches  deterministic  simulation  is  used,  and  thus  the  most  damaged 
subcircuit/transistor  is  obtained  assuming  that  a  user  specified  input  waveform  is  continuously 
applied  to  the  circuit.  Since  HCE  is  a  long  term  process,  the  combined  effects  of  a  large 
number  input  waveforms  need  to  be  applied  to  find  out  which  transistors  are  expected  to  accu¬ 
mulate  sufficient  degradation  so  as  to  affect  the  performance  of  the  entire  circuit.  To  deter¬ 
mine  the  most  d  imaged  transistor  with  respect  to  all  possible  input  waveforms  using  conven¬ 
tional  determini  St  c  timing  simulators,  exhaustive  simulation  or  at  least  a  ven/  large  number  of 
simulations  have  to  be  dore. 

In  this  work,  instead  of  performing  exhaustive  or  Monte  Carlo  simulation,  we  apply  pro¬ 
babilistic  simulation  techniques  to  estimate  HCE.  In  this  application,  we  are  interested  in 
computing  the  average  current  flowing  in  individual  transistors  within  subcircuits  during 
switching  time.  Thus  more  accurate  and  detailed  subcircuit  analysis  is  required.  To  attain 
this  level  of  accuracy,  the  probabilistic  voltage  waveform  representation  is  extended  to  include 
finite  t'-ansition  times  during  switching,  and  each  subcircuit  is  analyzed  in  more  detail,  as  will 
be  explained  in  Section  3,3.  With  this  additional  information  we  are  able  to  estimate  the  level 
of  expected  degradation  individual  transistors  might  experience,  relative  to  other  transistors  in 
the  circuit,  as  well  as  predict  the  effects  of  such  degradation  on  the  circuit  performance.  Our 
aim  is  not  to  estimate  how  long  it  takes  a  transistor  to  reaci^  a  certain  level  of  degradation,  but 
rather  to  suggest  design  changes  so  as  to  reduce  circuit  ui  gradation  that  may  be  caused  by 
HCE. 

In  Section  3.2,  we  de.scribe  the  measures  we  use  to  estimate  for  hot-carrier  effects.  In 
Section  3.3,  we  briefly  review  the  extension  of  the  input  irobability  waveform  to  include 
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transition  times,  and  explain  the  analysis  of  a  siniple  gate.  Voltage  signal  overlap  handling  is 
given  in  Section  3.4.  Implementation  and  simulation  results  are  presented  in  Section  3.5.  In 
Section  3.6,  we  propose  a  rede.sign  technique  derived  as  a  result  of  the  simulation,  for  reduc¬ 
ing  HCE  in  a  circuit  and  thus  improving  the  long  term  design  reliability  without  increasing 
the  design  area. 

3.2  Hot-Cs  rrier  Effects 

Hot-carrier  effects  in  MOS  transistors  are  caused  by  the  injection  of  high-energy  carriers 
into  the  gate  oxide  region  of  the  transistors  These  carriers  manifest  themselves  in  the  form 
of  oxide  charg  e  trapping  and  interface  trap  generation  which  accumulate  over  time,  depending 
on  the  operatic  n  of  the  transistor  within  a  circuit,  and  cause  a  permanent  change  in  the  oxide 
and  interface  charge  distribution  [14],  [15],  This  change  in  the  charge  distribution  causes 
degradation  in  {he  transistor  characteristics,  such  as  shift  in  the  threshold  voltage  V-j-,  and 
decrease  in  the  transconductance  gm  and  electron  mobility  |x„  in  the  channel  [14],  This 
degradation,  in  turn,  affects  the  performance  of  the  circuit  in  which  the  transistor  i.s  embed¬ 
ded,  such  as  increase  in  the  circuit  delay. 

With  current  semiconductor  process  technologies,  with  effective  channel  length  longer 
than  0.5  micron,  hot-camer-induced  degradation  is  more  severe  in  nMOS  transistors  than  in 
pMOS  [16].  Therefore,  in  this  study  we  focus  our  attention  on  n-channel  transistors.  Our 
approach,  however,  can  be  applied  as  well  to  estimate  HCE  in  pMOS  transistors  by  using  the 
appropriate  pMOS  degradation  model.  It  has  been  shown  that  HCE  in  nMOS  transistors  in 
digital  circuits  is  dominated  by  interface  trap  generation  which  occurs  mostly  when  the 
transistor  is  operating  in  or  near  the  saturation  region  [14],  [17],  !8],  and  thus  the  extent  of 
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the  damage  depends  on  the  operating  conditions  of  the  transistor  within  a  circuit  design.  It 
lias  also  been  shown  that  the  damage  can  be  measured  using  the  interface  trap  generation 
model  [19]: 

where  W  is  the  transistor  width  and  H  is  a  technology  dependent  parameter;  m  is  approxi¬ 
mately  3.  Thus,  to  estimate  HCE  induced  damage,  one  has  to  compute  the  drain-to-source 
current  and  substrate  current  waveforms,  losiO  and  /j„^(r),  for  the  lifetime  duration  of  the 
transistor  under  all  possible  inputs,  with  the  transistor  model  degraded  in  time.  In  our 
approach,  however,  we  are  interested  i*.  'Jentifying  those  transistors  that  are  most  susceptible 
to  HCE  induced  damage  for  redesign  purposes,  rather  than  estimating  the  lifetime  of  the 
design.  We  will  consider  the  average  value  of  D„  in  (3.1),  where  the  average  is  taken  over  all 
possible  inputs  to  the  circuit,  as  a  measure  of  relative  transistor  damage.  Alternatively,  we 
could  use  the  average  integral  of  I )  as  a  measure  of  damage.  As  we  will  show  below 
both  of  these  measures  agree  as  fai‘  as  detecting  relative  damage.  We  compute  these  averages 
using  the  undamaged  transistor  model. 

It  has  been  established  that  )  flows  when  the  transistor  is  operating  in  the  saturation 
region  and  is  a  function  of  the  drain-to-source  current,  /q.;,  drain-to-source  voltage,  with 
exponential  dependency  on  and  process  paraineters,  including  the  channel 

length  [14].  By  performing  exhaustive  and  detailed  circuit  simulation  of  typical  complemen¬ 
tary  CMOS  gates,  we  have  identified  a  number  of  factors  that  affect  the  time  interval  a 
transistor  stays  in  saturation  during  switching  and  consequently  the  current  waveform 

in  the  transistor.  These  factors  are;  (1)  the  type  of  gate  in  which  the  transistor  is  embedded; 
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(2)  the  frequency  of  transistor  switching  that  causes  the  gate  output  to  switch;  (3)  the  input 

slew  rate  and  the  gate  output  loading;  and  (4)  the  position  of  the  transistor  within  the  gate. 

We  will  now  justify  each  of  these  factors  in  more  detail. 

1.  Since  I exponentially  proportional  to  (V/35  -  if  transistor  is  prevented 

from  entering  saturation  region,  or  at  least  the  voltage  across  the  channel  is  reduced,  the 
substrate  current  can  be  reduced  dramatically.  As  a  consequence,  the  nMOS  transistors  in 
a  NAND  gate,  for  example,  suffer  less  than  those  in  a  NOR  gate  configuration  because 
the  serial  connection  of  nMOS  transistors  reduces  ^DS  of  each  nMOS,  and  hence  hot- 
carrier  generation  is  reduced. 

2.  Since  HCE  occurs  when  the  transistor  is  in  the  saturation  region,  which  in  CMOS  cir¬ 
cuits  happens  only  during  gate  output  transition,  it  follows  that  the  higher  the  switching 
frequency  of  a  transistor  in  a  gate  that  causes  the  gate  output  to  switch,  the  more  damage 
the  transistor  will  experience  over  a  period  of  time. 

.3.  The  input  slew  rate  and  the  gate  output  capacitance  also  affect  the  duration  in  which  a 
transistor  stays  in  saturation;  slower  slew  rate  causes  wider  I time-span,  and  hence 
more  degradation.  Similarly,  larger  output  capacitance  increases  output  transition  time 
and  thus  the  duration  the  transistor  stays  in  saturation.  F  gure  3.1  shows  the  relationship 
between  HCE  induced  damage  (D„  and  integral  of  I and  the  input  slew  rate  of  a 
CMOS  inverter.  Note  that  if  the  input  is  considered  to  be  a  step  function,  while  in  actu¬ 
ality  it  has  a  slew  rate  of  2  ns,  the  HCE  induced  damage  D„  could  be  underestimated  by 
twenty  times!  Figure  3.1  also  shows  the  relationship  between  output  capacitance  and 
relative  damage.  We  note  here  that  the  dependence  of  HCE  on  input  slew  rate  and  out- 
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put  capacitance  has  also  been  reported  in  [20].  Note  that  both  D„  and  show  the 
same  trend.  They  are  both  sublinear  functions  of  the  slew  rate.  This  fact  justifies  using 
average  slew  rate  valued  in  the  probabilistic  simulation  method  described  in  Section  3.3. 

4.  The  position  of  the  transistor  within  the  gate  also  affects  its  susceptibility  to  HCE.  Con¬ 
sider,  for  example,  the  NAND  gate  shown  in  Fig.  3.2.  Since  hot-carrier  induced  degrada¬ 
tion  in  the  nMOS  transistors  occurs  mainly  only  when  switches  from  high  to  low, 
as  will  be  shown  below,  three  cases  may  occur:  (i)  Ml  and  M2  turn  on  simultaneously, 
(ii)  M2  turns  on  last,  and  (iii)  Ml  turns  on  last.  Table  3.1  shows  and  sub  of  Ml 
and  M2,  with  different  capacitance  values  for  Cj,  after  a  simulated  three  months  continu¬ 
ous  operation.  It  is  clear  from  the  table  that  in  cases  (i)  and  (ii),  either  Ml  or  M2  experi¬ 
ences  light  degradation,  while  the  other  transistor  experiences  minimal  degradation,  This 
degradation  may  not  be  severe  enough  to  cause  circuit  malfunction.  On  the  other  hand, 
in  case  (iii)  Ml,  which  is  connected  to  the  output  node  and  switches  last  and  thus  is  the 
one  that  causes  the  output  to  switch  from  high  to  low,  suffers  extensive  degradation 
many  order  of  magnitude  compared  to  M2,  and  compared  to  both  M 1  and  M2  in  cases  1 
and  2.  These  results  can  l)e  explained  as  follows.  In  case  (i)  both  Ml  and  M2  are 
switching  at  the  same  time,  therefore,  the  drain-to-source  voltage  transistors 

IS  about  half  of  Vout-  mentioned  above,  smaller  due  to  serial  transistor  connec¬ 
tion  results  in  light  damage  to  both  Ml  and  M2.  In  cases  (ii)  and  (iii)  either  transistor 
causes  an  output  high-io-low  transition,  the  initial  value  of  across  M  i  is  while 
that  across  M2  is  Vj,/;  -  V-/ .  This  weak  high  stale  of  M2  is  the  main  factor  protecting  it 
from  .severe  degradation.  As  a  result.  Ml  stays  in  the  saturation  region  longer,  has  higher 
peak  substrate  current,  and  hence  experiences  more  degradation.  This  is  also  true  when 
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more  than  two  transistors  are  in  series;  a  transistor  connected  to  the  gate  output  stays  in 
saturation  during  switching  for  a  longer  period  of  time  compared  to  when  it  is  not  con¬ 
nected  to  the  output  node.  Thus,  we  conclude  the  top  nMOS  transistors  that  are  directly 
connected  to  the  output  node  have  the  potential  of  experiencing  severe  damage.  This 
observation  has  also  been  reported  in  [21]. 


Figure  3.3(a)  illustrates  the  relation  between  V^,  and  of  the  nMOS 


i, 
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transistor  of  a  typical  CMOS  inverter  obtained  using  SPICE.  We  have  found  that,  by 
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Figure  3.2  A  CMOS  NAND  gate  circuit  and  three  input  combinations  leading 

to  hot-carrier  induced  damage. 


(b)  Time 


Figure  3.3  (a)  The  waveforms  of  I  sub 

of  the  nMOS  transistor  of  a  CMOS  inverter  obtained  using  SPICE. 

(b)  The  waveforms  after  replacirig  with  a  ramp.  Th'*  ramp 
is  obtained  by  connecting  the  two  points  (r  1)  =  4)  and  =  1). 

approximating  in  Fig.  3.3(a)  with  a  ramp,  as  drawn  in  Fig.  3.3(b),  the  error  of  7^^;,  peak 

value  is  within  2%.  Compared  with  possible  modeling  errors,  this  result  is  accurate  enough 

especially  that  our  goal  is  to  obtain  the  relative  degradation  levels  of  the  transistors  in  CMOS 

circuits  rather  than  the  absolute  degradation  of  each  transistor.  We  will  thus  use  ramp 
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approximation  to  the  voltage  waveforms  during  transition.  Note  also  from  Fig.  3.3  that  HCE 
degradation  in  the  nMOS  transistor  during  output  low-to~high  transition  is  negligible,  and 
hence  we  estimate  HCE  degradation  only  during  output  high-to-Iow  transitions. 

33  Gate  Probabilistic  Simulation 
33.1  Probabilistic  Waveform  Description 

As  described  in  [1],  in  probabilistic  simulation,  a  signal  waveform  is  specified  in  terms 
of  its  probability  of  being  high  during  a  time  interval  and  the  probability  of  its  switching  from 
low  to  high  at  specified  time  points.  In  this  section  we  extend  this  description  by  including 
the  signal  transition  times  during  switt  ang.  These  transition  times  are  needed  to  estimate 
HCE  more  accurately,  as  shown  in  the  previous  section.  Figure  3.4  shows  a  typical  proba¬ 
bilistic  waveform,  description.  For  an  event  occurring  at  time  Iq,  five  values  ar  specified: 
Ph(t^),  PhitQ),  Plh{tQ).  tufito)  and  tHiAh)'  where  Ph(tQ  )  is  the  probability  of  the  signal 


Normalized 

transistor 

case  (i) 

case  (ii) 

case  (iii) 

Cl  =  0.2pf 

1.28e-3 

8.95e-8 

0.84 

Cl  --  0.02pf 

1.30e-12 

3.2Se-5 

1.30e-12 

Cl  =  0.2pf 

3.38e-8 

1.40e-10 

1.223 

Cl  =  0.2pf 

1.47e-12 

6.29e-t) 

1.47e-12 

Cl  =  0.2pf 

BS 

9.99e-ll 

1.48e-12 

1.42 

Cl  =  O.Spf 

1.50e-12 

4.84e-7 

1.50e-12 

I  I jub 

transistor 

case  (i) 

case  (ii) 

case  (iii) 

Cl  =  0.2pf 

Ml 

2.104e-2 

1.412e-03 

1.291e-01 

Cl  =  0.02pf 

M2 

1.260e-06 

1.2768-03 

1.260e-06 

Cl  =  0.2pf 

l.Olle-03 

1.758e-04 

1.9778-01 

Cl  =  0.2pf 

■91 

2.848e-06 

1.985e-03 

2.848e-06 

Cl  0.2pf 

Ml 

1.192e-04 

1.123e-05 

2.420e-01 

Cl  =  O.Spf 

M2 

3.784e-06 

1.135e-03 

3.781e-06 

Table  3.1  The  damage  of  nMOS  transistors  in  an  NAND  gate 
with  various  internal  capacitances. 


37 


l)eing  high  right  before  Iq,  Phi:^)  the  probability  of  the  signal  being  high  after  Iq  +• 


I 

i 

I  Plh(tQ)  the  probability  of  the  signal  switching  from  low  to  high,  and  (^q)  the  aver 


age  low-tO“high  and  high-to-low  transition  times  of  the  signal.  If  there  is  a  complete  transi¬ 
tion,  then  the  voltage  waveform  (a  ramp)  can  be  reconstructed  as 


V(0=1 


^DD  * 


t~to 


if  Vfto  )  “  0,  /jj  <  r  ^  to 

if  V(to  )  =  V^p,  to  S  r  fo  +  ^//£,(to) 


(3.2) 


Note  that  Phlito),  the  probability  of  switching  from  high  to  low,  can  be  derived  from  [1], 


PhiitQ)  =  Ph{to)  +  Plhito)  -  Ph{t^)  (3.3) 

/Jso  note  that  Iq  is  an  average  reference  time  point  and  not  the  actual  time  when  a 
transistor  switches.  The  timing  simulation  algorithm  descri'oed  below  computes  the  average 
transistor  switching  time  relative  to  tQ  based  on  its  value. 


3.3.2  Subcircuit  Probabitistk  Simulation  Including  Transition  Times 

We  consider  for  the  moment  static  CMOS  circuits  that  consist  of  serial-parallel 
configurations  with  no  pass  transistors.  We  define  a  gate  to  be  a  channel-connected  subcircuit. 
The  simulation  of  a  single  gate  consists  of  two  parts:  One  is  steady-stare  analysis,  which 

tLjj(tO) 

Plh(to) 


Ph(to“ ) 

PhtV) 

to 

Figure  3.4  A  typical  input  waveform  for  probabilistic  timing  simulation. 


involves  the  coirjputation  of  Ph{t)  at  the  gate  output  node  in  the  time  intervals  before  it  starts 
and  after  it  completes  switching.  The  other  is  switching  analysi.s  which  involves  the  computa¬ 
tion  of  the  average  delay  A/.  Plh{t^),  Phl{t^),  ryjy(t,).  when  a  subcircuit  output 

switches,  where  tj  =  tg  -f  A/. 

Calculating  Ph{t)  of  the  output  node  is  equivalent  to  finding  the  probability  of  haviitg  a 
conducting  path  from  the  output  node  to  Vqj)  ;  i.e. 

Ph{i\^  -  -P„„(f.e,-(oucV£,c))  -  1  -  P^„(r.e2  =  (out,G^o))  (3-4) 

Where  P„„(t.e)  represents  the  probability  of  having  a  conducting  path  through  equivalent 
edge  e.  If  the  edge  represents  a  transistor,  then  P^„{t ,e)  -  Ph(t)  (or  Pl(i))  of  the  gate 
node  of  the  nMOS  (pMOS)  transistor.  For  serial-parallel  CMOS  circuits,  the  graph  reduction 
orocedure  consists  of  serial  and  parallel  combinations  of  edges  as  described  in  [1]. 

In  switching  analysis  we  follow  the  approach  described  in  Chapter  2,  where  signals  arriv¬ 
ing  at  the  inputs  of  a  gate  are  considered  to  switch  at  separate  time  instances,  even  if  they 


Vao 


Figure  3.3  A  fully  complementary  CMOS  circuit. 
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appear  to  arrive  closely  together.  Events  arriving  at  a  gate  tire  propagated  individually  and 
merged  at  the  output  if  they  occur  there  within  a  certain  time  interval;  and  when  an  nMOS 
transistor  in  a  CMOS  gate  switches  and  its  switching  causes  a  transition  at  the  gate’s  out¬ 
put,  then  at  least  one  conducting  path  is  formed  between  the  output  and  ground,  and  all 
trairsistors  that  form  parallel  paths  with  must  be  off.  The  same  observation  holds  for  the 
pMOS  transistoR  with  ground  replaced  by  • 

Consider  the  nJMOS  part  of  a  CMOS  logic  gate  thown  in  Fig.  3.5,  where  the  capacitance 
C  is  lumped  at  the  output  norle.  Depending  on  the  position  of  the  switening  transistor,  the 
nMOS  part  is  reduced  to  one  of  the  three  macromodels  as  shown  in  Fig.  3.6.  In  Fig.  3.6(a) 
the  switching  transistor  is  connected  to  the  output  node;  in  Fig.  3.6(b)  it  is  connected  to 
ground;  and  in  Fig.  3.6(c)  it  is  connected  to  other  transistors  and  not  to  or  ground  All  of 
the  nMOS  transistors,  other  than  the  switching  transistor,  are  combined  into  equivalent  transis- 
toR.  The  transconductance  parameters  p,  and  P2  of  ihe  equivalent  transistors  are  random 
variables,  with  expected  values  derived  using  series/parallel  probabilistic  reduction  techniques 
given  in  Chapter  2.  The  same  method  is  applied  to  pMOS  part  of  the  logic  circuit.  Finally, 


V«.» 

Vo.« 

V<»i 

v»o  P’ 

i  •  , 

vdo  H  1  pi 

n 

V(«  hl 

Ht 
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Vto  — jn 

vbd  —j  jr 

_ 11” 

J  switching 

4.- 

4. 

4* 
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(a) 

(b) 

(C) 

Figure  3,6  Three  mad'ornodeis  for  the  nMOS  pait  of  a  CMOS  logic  circuit, 
depending  on  the  position  of  the  switching  transistor. 
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the  graph  reducticn  is  continued  for  both  the  n-  and  p-parts  to  obtain  an  equivalent  inverter 
macromodel  shown  in  Fig.  3.7. 


3.3.3  Propagation  Delay  and  Output  Slew  Rate 

The  approach  we  follow  is  to  first  calculate  the  output  voltage  waveform  and  the  current 
using  the  equivalent  ciicuit  obtained  in  Fig.  3.7.  The  output  waveform  is  then  applied  to 
the  appropriate  equivalent  circuit  shown  in  Fig.  3.6(a),  (b),  or  (c),  to  estimate  the  substrate 
current  waveform  in  the  switching  transistors.  The  same  method  can  be  applied  to  the  pMOS 
part  as  well,  if  necessary.  The  equation  describing  the  equivalent  inverter  circuit  in  Fig.  3.7 
can  be  written  in  the  form: 


(3.5) 


dt  C 

Equation  (3.5)  is  solved  with  Backward  Euler  integration  to  obtain  the  output  waveform  Vg^^, 
and  the  current  waveforms  ip  and  Other  faster,  but  less  accurate  methods  such  as  [13], 

[22],  or  [23]  can  be  applied  as  well.  During  the  simulation,  two  voltage-time  points  on  the 
output  waveform,  (V^o  ^a)  and  (Vy^,  r^),  are  selected  in  order  to  calculate  the  delay 

and  output  slew  rate. 


Vdu 


Figure  3.7  The  equivalent  inverter  circuit  of  the  CMOS  gate  in  Fig.  3.5. 


41 


Based  on  our  ramp  waveform  definition,  for  output  high-lo-low  transition, 


^HL  -  (h~^a)  * 


DD 


^DD  -Vj-  ~\  Vj 


\Vt 

‘  D 


a  ~^HL 


'DD 


-]  -  /o 


(3.6) 


(3.7) 


Similai-  equations  are  applied  to  output  low-to-high  transition.  In  general  the  propagation  delay 
of  low-to-high  and  high-to-low  transitions  are  different.  In  order  to  simplify  the  problem,  we 
take  the  weighted  average  of  the  up  and  down  transition  propagation  delays  using  their  "pro¬ 
bability  to  occur"  as  weights  [1];  i.e.. 


At. 


average 


Ati,Plh(t)^^,+At,iPhlit)^^, 

PiHt)ou,+Phl(t)ou. 


(3.8) 


In  this  way,  the  up  and  down  transitions  of  each  event  are  always  bound  together  a^ter 
the  signal  is  propagated  through  the  circuit.  This  binding,  however,  is  not  necessary  for  the 
approach.  The  cost  is  added  storage  and  computation.  The  output  waveform  (defined  by 
Ph{t{),  Phitp),  PlhitO,  where  h  =  h ^t^verage  and 

r =  r ,  +  maxitf^i^it ^),  ,)))  now  becomes  an  input  to  the  fanout  subcircuits. 


3.3.4  Substrate  Current  Estimation 

After  the  computation  of  the  output  node  voltage  waveform  and  the  current  waveform 
through  the  equivalent  inverter  circuit,  at  each  time  point  and  D„  in  the  switching  transis¬ 
tor  are  computed  using  one  of  the  macromodels  in  Fig.  3.6.  For  example,  in  Fig.  3.6(a), 
applying  KCL  at  the  internal  node  1,  the  following  equation  applies: 

(  ^oul ' ^in  » ^  I )  ~  ( ^  1  >  ^^DD  '^ND  )  (3.9) 
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In  this  equation,  the  values  of  i„ ,  ,  and  are  known  at  each  time  point.  The  only 

unknown  variable  is  V].  Note  also  that  the  bottom  equivalent  transistor  in  Fig.  3.6(a)  never 
enters  into  the  saturation  region.  The  procedure  of  estimating  and  D„  in  Fig.  3.6(a)  is 
listed  as  follows: 

(i)  At  each  lime  step,  Eq.  (3.5)  is  solved  with  Backward  Eiuler  integration  to  obtain  the 
values  of  and  i„ . 

(ii)  The  value  of  i„  is  used  to  find  the  voltage  Vj  across  the  equivalent  transistor  p2> 
the  linear  region  current  equation  of  the  equivalent  transistor. 

(iii)  The  voltage  ^^5.  =  across  the  switching  transistor  is  used  to  check  if  the 

transistor  is  in  saturation  region.  If  so,  and  are  applied  to  estimate  the 

substrate  current  of  the  switching  transistor  at  the  particular  time  point.  The  trapezoidal 
method  is  employed  to  approximate  D„  and  the  integral  of 

(iv)  Repeat  (i)  to  (iii)  until  the  switching  transistor  no  longer  stays  in  saturation  region. 

Similar  .ocedures  are  applied  to  the  macromodels  in  Fig.  3.6(b)  and  (c).  The  output 
voltage  wavefomi  and  current  waveform  are  used  to  deteiinine  the  voltage  across  the 
switching  transistor,  check  its  operating  condition,  and  calculate  its  substrate  current  waveform 
as  well  as  £)„ .  It  follows  that  the  average  damage  the  switching  transistor  experiences  with 
respect  to  one  incoming  event  is 


=^''' j  l„s(~rd,  (3.10) 

Yvn  ^  DS 


Note  that,  as  mentioned  in  previous  section,  the  HCE  degradation  during  low-to-high 
transitions  can  be  ignored.  Also  note  that  during  the  operation  of  the  circuit  one  transistor 
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may  be  driven  by  several  incoming  events.  'Fherefore,  the  average  total  damage  the  switching 
transistor  experiences  during  the  entire  period  of  operation  becomes 

or  (£>„)„„,  (3.11) 

DS 

3.4  Signal  Overlap  Handling 

When  voltage  waveform  slew  rate  is  included  in  the  waveform  specification,  overlap  dur¬ 
ing  transition  could  occur  and  needs  to  be  taken  into  account.  Figure  3.8  shows  the  relation 
between  two  consecutive  events  at  a  gate  output  caused  by  two  different  switching  transistors. 
In  Fig.  3.8(a),  the  two  events  do  not  overlap,  so  they  can  be  handled  independently.  In  Fig. 
3.8(b),  the  two  events  are  close  to  each  other,  and  thus  tliere  is  certain  possibility  that  one 
transition  happens  before  the  other  transition  is  completed.  :\s  shown  in  this  figure,  the  high- 
to-low  (low-to-high)  transition  of  the  second  event  may  occur  before  the  low-to-high  (high- 
to-low)transition  of  the  first  event  reaches  steady  state.  Assuming  these  two  events  at  the  out¬ 
put  node  are  independent,  the  corresponding  probability  of  having  overlap  is  as  follows. 


Figure  3.8  (a)  No  overlap  between  two  events;  (b)  two  events  overlap  each  other. 


r 

JO  //  ^2-^!  +0^ 

ProH^-,ype  ^  if 

I 


Proh{V—type  overlap)—'^ 


0 


If  r2>^i  +  ^//L 


Ph^i \)event\  *  P  ^  ^(h)  event  2^  P^^l  ^eventl  +^WL 


(3.13) 


However,  the  two  output  events  are  in  general  dependent.  Depending  on  the  relative 
position  of  the  switchii  g  transistors,  only  one  of  the  two  overlap  types  would  occur.  If  the 
switching  transistors  ;ire  in  series  in  n-part  of  a  CMOS  gate  (e.g.,  a  NAND  gate),  then 
V-type  overlap  would  occur  with  probability  given  in  (3.13),  but  A— type  overlap  will  not 
occur;  if  the  switching  transistors  are  in  series  in  p-part  of  a  CMOS  gate  (e.g.,  a  NOR  gate), 
then  A~type  overlap  would  occur  with  probability  given  in  (3.12),  but  V-type  overlap  will 
not. 


The  reason  is  stated  as  follows;  Consider  the  case  when  the  switching  transistors  are  in 
series  in  n-pan  of  a  CMOS  gate.  In  order  to  have  a  A-type  overlap,  the  first  switching 
transistor  must  switch  on  — >  off  (the  output  node  switches  from  low-to-high),  and  the 
second  switching  transistor  needs  to  switch  off  on .  However,  when  the  first  switching 
transistor  is  causing  an  output  transition,  the  second  switching  transistor  needs  to  be  conduct¬ 
ing  since  it  is  in  series  with  the  first  one.  Therefore,  the  second  switching  transistor  could  only 
switch  on  off  after  the  first  event  is  propagated  through,  which  is  contradictory  to  the 
condition  of  having  a  A-type  overlap,  and  hence  A-type  overlap  could  not  occur  when  the 
switching  transistors  are  in  series  in  n-part  of  a  CMOS  gate.  Similar  reasoning  can  be  applied 
when  the  switching  transistors  are  in  series  in  p-parl  of  a  CMOS  gate. 
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We  now  use  an  example  to  explain  overlap  in  more  detail.  Consider  the  NAND  gate 
shown  in  Fig.  3.9(a).  All  possible  waveforms  incident  at  input  nodes  A  and  B  are  displayed 
within  the  shaded  areas  in  Fig.  3.9(a).  Four  waveforms  with  equal  probabilities  are  possible 
at  each  A  and  B.  The  expected  values  of  the  switching  times  at  A  and  B  are  and  t^, 
respectively.  The  corresponding  probability  v/aveforms  are  shown  directly  above  the  shaded 
areas.  The  slew  rate  information  has  been  omitted  from  the  probability  waveform  figures  for 
clarity.  Figure  3.9(b)  shows  within  the  shaded  area  all  the  possible  output  waveforms  and 
their  frequency  of  occurrence  out  of  the  sixteen  possible  scenarios  obtained  using  exhaustive 
simulation.  Output  switching  time  1 1  is  equal  to  plus  the  expected  value  of  the  gate  delay 
when  A  is  switching,  and  t2  is  equal  to  plus  the  expected  value  of  the  gate  delay  when  B 
is  switching.  In  general,  these  delay  values  could  be  different,  depending  if  A  or  B  switches. 

The  results  show  that  the  output  signal  at  C  will  switch  at  t  ]  twice  from  low  to  high  and 
once  from  high  to  low.  It  will  switch  at  ^2  ottce  from  low  to  high  and  twice  from  high  to 
low.  In  one  case,  it  will  switch  from  high  to  low  at  tj  and  back  to  high  at  time  ^2-  Itt  eight 
cases  out  of  sixteen,  it  will  stay  high  md  in  one  case  it  will  stay  low.  The  probabilistic 
waveform  above  the  shaded  area  in  Fig.  3.9(b),  which  is  obtained  directly  by  performing  pro¬ 
babilistic  simulation  on  the  gate,  matches  the  statistics  obtained  from  exhaustive  simulation. 
Figure  3.9(c)  show  the  results  when  is  closer  to  so  that  ^2  “  less  than  the  gate 
delay.  In  this  case,  an  incomplete  glitch  occurs  at  C,  as  shown  in  waveform  number  5  in  Fig. 
3.9(b).  Incomplete  glitches,  or  overlapped  output  waveforms  are  detected  in  probabilistic 
simulation  by  monitoring  *2  slew  rate  at  t If  ^2  -  is  les,s  than  the 

time  it  takes  a  waveform  starting  at  to  reach  steady-state,  then  an  incomplete  glitch  might 
occur.  This  case  is  shown  in  Fig.  3.9(c). 


no  overlap 


overlap 
0.125  0.125 

r  0-75  4o.7^  0.75  ►Merge 
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Figure  3.9  (a)  A  NAND  gate  and  the  input  waveforms; 
(b)  the  output  waveforms  when  tnere  is  no  overlap; 
(c)  the  output  waveforms  when  there  is  overlap. 


Overlap  can  cause  incomplete  transitions  and  hence  two  issues  need  to  be  addressed:  ( 1 ) 


the  overlapped  waveform  may  not  propagate  beyond  the  fanout  gates;  and  (2)  the  overlapped 


waveform  may  affect  the  degradation  within  the  subcircuit.  Whether  an  overlapped  waveform 


will  propagate  depends  on  the  voltage  of  the  intersection  point  of  the  overlapped  events,  as 
shown  in  Fig.  3.10.  If  the  intersection  point  is  close  to  for  A-type  waveform  (or  G/yo  for 
V-type),  then  the  two  events  can  be  considered  as  two  complete  transitions  by  the  fanout 
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gates.  On  the  other  hand,  if  the  intersection  point  i.s  too  low  (or  too  high  for  V-type),  then  it 
is  only  a  glitch  for  the  fanout  gates  and  the  overlapped  wavefoim  will  not  be  propagated. 

Based  on  circuit  simulation  results  on  typicaJ  designs,  a  threshold  of  V'po/2  is  chosen  to 
filter  out  some  overlapped  waveforms;  that  is,  if  the  intersection  point  of  a  A-type  (V-type) 
waveform  is  VddI2  or  higher  (lower),  then  the  overlapped  waveforms  are  propagated  as  if 
there  is  no  overlap;  otherwise  the  two  events  are  considered  overlapped  and  merged  into  one 
before  being  propagated  to  the  fanout  gates.  Note  that  only  events  on  the  same  signal  line 
would  be  merged.  Figure  3. 1 1  shows  the  merging  of  two  adjacent  events  with  the  following 
equations; 

Plh(a)  Plh(l)  +  Plh{2)  ~  Prob(A  ~  type  overlap)  (3.14) 

if  events  I  and  2  are  caused  by  two  transistors  in  series  in  the  p-part  of  the  gate.  Or 

P lh(a)  =  P !h(\)  +  P I h(2)  -  Prob(V  -  type  overlap)  (3.15) 

if  events  1  and  2  are  caused  by  two  transistors  in  series  in  the  n-part  of  the  gate. 


Figure  3.10  Two  types  of  overlapped  waveforms  and  the  intersection  points. 
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Figure  3.1 1  The  merging  of  events 
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{p/;i(i)+p/i/(i))x?i+(p/;i(2)+/>;z/(2))xr2 

P//i(l)  +  P/i/(r)  +  P/A(2)  +  P/z/(2) 


In  addition. 


Phl{\)tni^{{\)+PhH2)t„^X2)) 

Phl{\)+Phl{l) 

Pf;>(i)f^((i)+p/;i(2)f^^(2)) 

Plh{\)-^Plh{2) 


(3,16) 


(3.17) 


(3.181 


The  effect  of  overlap  on  damage  estimation  is  that  an  overlapped  waveform  does  not 
draw  the  same  amount  of  current  as  two  complete  transitions  do,  resulting  in  less  degradation 
than  two  complete  transitions.  Based  on  circuit  simulation  exiTcriments,  we  have  found  tha 
the  degradation  in  nMOS  transistor  caused  during  the  falling  edge  of  an  incomplete  A-type 
glitch  is  reduced,  according  to  the  following  equation. 


(^n  )ave  ={Phl  - kProh (A  -  type  overlap ))  —■  J  (-~T  dt  (3.19) 

WH  Ids 

where  Ph  I  is  the  output  high-to-  low  probability  caused  by  the  second  event  before  rnergin:', 
and  where  ^  is  a  function  of  the  intersection  point  voltage  V , 


k  =  J 


2  - 


2V 


'DO 


Vdo/2  <V  <  Vdo 


^  -  ^00^2. 


(3.2,0) 


The  modification  applies  to  the  calculation  of  (j  hubXive  as  well.  Note  that  damage  occurs  in 
the  nMOS  transistor  during  the  falling  edge  of  a  wavefonn  and  none  during  the  rising  edge. 


3.5  Implementation  and  Simulation  Results 


We  have  implemented  the  above  approach  in  program  iProbe-d.  The  program  consists  of 
four  main  parts:  (1)  Input  circuit  reading  subroutines,  (2)  partitioning  and  scheduling,  (3)  cir¬ 
cuit  simuhition  based  on  waveform  relaxation,  and  (4)  probabilistic  and  timing  sirniilalion. 
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iProbe-d  reads  in  SPICE-lypQ  circuit  description  and  creates  a  circuit  database,  the  input 
circuit  can  be  defined  hierarchicaliy.  Instead  of  flattening  the  input  circuit,  the  program 
preserves  the  circuit  hierarchy  to  save  memory  space. 

The  partitioning  and  scheduling  parts  of  the  program  are  flexible  in  the  sense  that  the 
user  may  choose  either  to  partition  the  circuit  or  keep  die  user  specified  partitions.  If  partition¬ 
ing  is  activated,  the  circuit  at  the  lowest  level  in  the  hierarchy  is  partitioned  into  channel- 
connected  subcircuits.  The  scheduling  part  orders  subcircuits  for  simulation  and  detects  feed¬ 
back  loops. 

The  circuit  simulation  part  is  implemented  for  the  situation  when  more  accurate  simula¬ 
tion  is  necessary  and  for  verifying  the  accuracy  of  the  results  during  the  program  development 
stage.  The  probabilistic  and  timing  simulations  are  combined  together.  In  fact,  timing  simula¬ 
tion  is  a  special  case  of  probabilistic  timing  simulation  by  specifying  the  values  of  Ph,  PI, 
Plh ,  Phi  to  be  f  or  I . 


3t 


Figure  J,12  A  test  CMOS  combinational  circuit  and  the  input  statistical 
descriptions,  and  are  both  1 .0  ns  for 
each  input  event.  Simulation  period  T  =  20  ns. 
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iProbe-d  has  been  used  to  simulate  several  simple  and  benchmark  circuits  [6],  Figure 
3.12  shows  a  CMOS  circuit  with  two  NAND  and  one  NOR  gates.  The  statistical  descriptions 
of  the  four  input  nodes  are  also  shown  in  the  figure.  Table  3.2  illustrates  the  degradation 
(D„)  of  the  six  nMOS  transistors  using  iProbe-d,  SPICE  exhaustive  simulation,  and  SPICE 
one  run,  respectively.  The  inputs  are  assumed  to  be  uncorrelated.  The  circuit  is  simulated  for 
a  time  period  T  =  20  ns;  the  input  events  are  then  assumed  to  repeat  periodically  and  the 
damage  is  extrapolated  to  three  months.  For  each  primary  input  node,  there  are  four  possible 
waveforms  (staying  high,  staying  low,  switching  from  low  to  high,  from  high  to  lov/);  hence, 
there  ure  4^^  256  possible  input  combinations.  Therefore,  for  the  SPICE  exhaustive  simula¬ 

tion,  256  runs  are  needed.  On  the  other  hand,  only  one  run  is  necessary  when  iProhe-d  is 
employed.  The  input  voltage  waveforms  of  the  third  row  in  Table  3.2,  SPICE  one  run,  are 
A  -  B  -  D  =  Vjjf-)  and  C  -  a  periodical  square  wave.  From  the  result,  iProhe-d  and  SPICE 
exhaustive  simulation  show  the  same  tendency,  while  running  SPICE  with  one  particular  input 
waveform  points  to  the  wrong  critical  transistor.  We  have  also  compared  the  waveform 
obtained  from  iProbe-d  and  SPICE  exhaustive  simulation  for  this  circuit.  The  slew  rate  and 
propagation  delay  from  iProbe-d  are  within  5%  error  compared  to  the  average  '/alue^  from 
SPICE  exhaustive  simulation  on  all  the  internal  nodes. 

Table  3.2  demonstrates  the  effect  of  overlap  as  wed.  For  the  two  NAND  gates,  as 
pointed  out  in  previous  section,  there  exists  only  V-type  overlap  and  hence  the  degradation  in 
the  top  nMOS  transistors  in  these  two  gales  would  not  be  affected.  Both  iProhe-d  and  SPICE 
exhaustive  simulation  prove  the  prediction.  At  the  NOR  gate,  there  exists  A  -  type  overlap 
and  the  right  transistor,  which  switches  last,  ha.s  less  degradation  than  the  left  one  w'hen  there 
is  overlap.  The  decreasing  degradation  in  the  left  transistor  in  the  NOR  gate  is  the  result  of 


overlapj  in  the  NOR  gate  as  well  as  the  event-merging  at  node.  E.  A  number  of  ISCAS85 
benchmark  circuits  [6]  have  also  been  analyzed.  These  benchmark  circuits  have  been  con¬ 
verted  to  SPICE  input  file  with  complementary  CMOS  transistor  circuit  implementation  of 
each  gate.  Table  3.3  shows  the  computation  time  on  a  SPARC  station  )0,  together  with  the 
damage  estimation  after  a  simulated  three  months  continuous  oiTeration.  Here  we  assign  otte 
event  to  each  primary  input  node,  with  50ns  as  the  simulation  period,  and  assume  the  primary 
input  signals  are  uncorrelated.  We  have  also  simulated  C432  with  the  circuit  simulation  mode 
of  iProbe-d,  which  in  this  case  is  faster  than  SPICE^  since  the  circuit  is  combinational  with 
no  feedback,  the  simulation  time  for  one  input  vector  specification  on  a  SPARC  station  10  is 
[  1  seconds.  This  circuit  has  36  inputs  and  exhaustive  simulation  would  require  4^^  x  1 1 
seconds  =  1  65  x  lO’^  years!  The  results  show  that  iProbe-d  can  handle  large  circuits  with 
reasonable  computation  time.  The  results  shown  in  Table  3.3  ignore  signal  correlations  in 
analyzing  the  gates  within  the  circuit,  To  take  into  account  signal  coiTelation  effects,  the 
method  proposed  in  [5]  and  modified  in  [4]  to  nulude  delay  information  is  being  imple¬ 
mented. 

3.6  Transistor  Degradation  Model  and  Effects  on  Circuit  Performance 

So  far,  we  have  estimated  relative  HCE  degradation  in  a  MOS  transistor  by  estimating  its 
substrate  current  which  i.*:’  a  function  of  drain  curreni  drain  to  source  voltage 

Vqv.  The  actual  physical  damage  due.  to  HCE  ir*  nMOS  transisters  is  due  in  great  part  to 
interface  charge  generation,  which  can  be  represented  by  the  inferface  state  density  A',..  The 
relation  between  /  and  is  grve:i  b’'  (14],  [24]: 

'If  (he  circuit  has  lecdbcck,  the  Jctemirmsiic  ciixuit  simulation  niadtr  of  i Probe  d  Lcrs  waveform  reiaxa'io'';  i  .  this  case,  iProbc-d  may 
or  may  not  be  faster  lh.in  SPICE  iii  ;uialyr:irig  the  circuit  for  one  input  waveform  (depending  an  the  rale  of  cr'nve'i’eoce) 
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NANDI 

top 

NANDI 

bottom 

NAND2 

top 

NAND2 

bottom 

NOR 

left 

NOR 

right 

SPICE  exhaustive 
r  =  5ns 

2,969 

7.56e-5 

1 

2.975 

7.80e-5 

2.285 

2.280 

T  —  1.5ns 

2.966 

5.63e-5 

2.975 

5.65e-5 

1.718 

0.911 

T  —  O.Sns 

2.976 

4.63e-5 

2.978 

4.73e-5 

1.648 

0.397 

iProbed 
r  =  Sns 

6.28e-5 

6.28e-5 

2.045 

2.045 

r  —  1.5n3 

S.28e-5 

6.28e-5 

1.713 

1.111 

r  =  0.5n5 

2.870 

6.28e-5 

2.870 

6.28e-5 

1.535 

0.384 

0 

0 

6.43 

0 

0 

9.62 

Table  3.2  The  relative  damage  of  all  the  nMOS  transistors  in  the  circuit 
shown  in  Fig.  12.  '''A  =  B  =  D  =  high(V£,£j); 

C  =  a  periodical  square  v/ave. 


ISCAS85 

Circuits 

I  #  of  nMOS  transistors  in  each  damage  level  (  Normalized  Z?n)  1 

CPU 

time  (sec) 

<  1.0 

1.0  -  2.0 

2.0  -  3.0 

3.0  -  4.0 

4.0  -  5.0 

>  5.0 

C432 

232 

37 

16 

65 

9 

53 

1.00 

C499 

576 

0 

8C 

0 

0 

226 

1.65 

C880 

443 

62 

40 

39 

52 

265 

3.37 

C1355 

776 

176 

56 

32 

0 

114 

4.43 

C1908 

653 

30 

29 

65 

86 

780 

7.60 

C2670 

971 

209 

151 

80 

125 

1146 

3.49 

C3540 

1450 

340 

225 

171 

143 

1423 

17.12 

C5315 

2147 

274 

374 

261 

192 

238S 

18.37 

C6288 

512 

92 

495 

417 

156 

3384 

92.08 

C7552 

2767 

332 

861 

262 

191 

3285  j 

26.84  j 

Table  3.3  The  conjputation  time  of  iProbe-d  on  several  benchm^ark 
circuits,  and  predicted  average  damage  after  three  months  continuous 
operation.  The  simulation  is  done  on  a  SUN4  SPARC  station  10. 
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Where  lO'^V/cm;  Du  and  Xu  repiesen'"  diffusion  constants;  K  and  fi  represent 


process-dependent  coefficients;  where  is  the  oxide  thickness  and  Xj  is 

the  junction  depth. 

Since  HCE  can  cause  a  shift  in  Vj  and  degradation  in  gm  and  }i„ ,  it  may  result  in  indi¬ 
vidual  device  failure  as  well  as  in  overall  circuit  performance  degradation.  In  digital  circuits, 
tlie  latter  is  usually  an  increase  of  circuit  delay.  We  thus  set  two  criteria  to  evaluate  HCE 
degradation.  The  first  criterion  is  the  local  damage,  represented  by  N^, ,  in  each  nMOS  transis¬ 
tor;  and  the  second  criterion  is  the  increase  in  total  delay  in  the  circuit,  which  can  be  deter¬ 
mined  by  searching  the  longest  true  path  before  and  after  degradation.  If  either  criterion 
exceeds  a  user  specified  limit,  it  is  then  necessary  to  redesign  the  circuit. 

In  estimating  the  circuit  delay  we  employ  a  linearized  RC  model.  In  each  subcircuit,  the 
transistors  are  replaced  by  equivalent  resistors.  The  subcircuit  capacitors  are  lumped  together 


VDD  .  VDD 


(a)  (b) 

Figure  3.13  (a)  A  CMOS  circuit  raid  (b)  its  equivalent  RC  network  and  worst-case  delay  paths; 
the  an-ows  represent  longest  resistive  paths  in  the  n-patl  and  p-part  of  the  subcircuit. 
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at  the  output  node.  Figure  3.13  sh(!W.s  a  CMOS  circuit  and  it.s  equivalent  RC  network. 
Worst-case  delay  estimation  is  assumed.  In  this  case,  the  worst-case  rise  and  fall  delay  of  the 
circuit  are; 


'^Drise  ~  (^/l,  *  ^load  (3.24) 

^Dfall  =  (^8,  *  ^load  (3.25) 


assuming  that 


and 


^/i„  -  ^s„  “ 


=  Rr  =R 


(3.26) 


•eq 


^\^nC„(.VDD-VT.) 


for  nMOS 
for  pMOS 


(3.27) 


The  equations  above  show  that  HCE  may  affect  the  nMOS  equivalent  resistance,  which 

depends  on  both  Vj  and  u.„ ,  but  the  capacitance  is  not  affected.  Earlier  work  on  HCE 

assumed  that  the  oxide-interface  charge  (Q„  =  *  ^ )  is  uniformly  distributed  over  the 

entire  channel  region  and  hence  the  change  of  threshold  voltage  is  linearly  proportional  to  the 

generated  interface  charge  [25].  However,  experiments  have  proved  that  the  damage  causod  by 

HCE  is  localized  in  the  channel  near  the  drain  node  (for  nMOS)  and  a  unifoun  model  cannot 

correctly  characterize  the  damaged  device.  A  one-dimensional  model  has  been  proposed  by 

Leblebici  and  Kang  for  nMOS  transistors  with  hot-carrier  induced  oxide  damage  [26].  In  this 

model,  the  transistor  channel  is  divided  into  two  regions,  the  damaged  and  undamaged  region. 

« 

As  reported  in  [26],  based  on  experimental  data  for  l|J.m  technology,  the  damaged  region  is 
located  within  0.1  |.mi  of  the  drain  node  and  the  interface  charge  is  uniformly  distributed  in 
this  region.  The  rest  of  the  channel  emains  undamaged  with  no  generated  interface  charge. 
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For  subniicron  technology,  the  length  of  the  damaged  region  needs  to  be  deduced  from  new 
experimental  data. 

We  adopt  the  concept  of  this  model  and  use  a  two-transistor  degradation  model,  as 
shown  in  Fig.  3.14,  where  a  damaged  nMOS  transistor  is  considered  as  two  transistors  con¬ 
nected  in  series,  with  different  threshold  volcagcs  and  mobilities.  The  channel  length  of  the 
two  transistors  are  0.1  pim  and  L  ~  0.1|im  respectively,  where  L  is  the  channel  length  of  the 
original  transistor  and  the  transistor  with  the  0.1  )i;n  length  is  considered  to  be  the  damaged 
one.  The  equivalent  conductances  of  the  transistors  are  calculated  with  the  following  equa¬ 
tions: 


3  *  10'‘  T 

u„.  =  1000  *  •  *  — 

Ni,  80 

(3.28) 

1  Hoot 

(3.29) 

Vr.-yr--^ 

^OX 

(3.30) 

(3.31) 

w 

Ri-  =G2=  -Vr) 

(3.32) 

L  |im,  jin,  Vt,  R 


0.1  |im,  jini,  Vti  Ri 


L-0.1  jim,  jin,  Vt,  R2 


(a)  (b) 

Figure  3.14  An  riMOSFET  (a)  before  and  (b)  after  HCE  degradation. 
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and  the  increase  of  fall  delay  is: 


*  Cl,,,  (3,33) 

Note  that  there  is  no  change  in  the  rise  delay  since  the  damage  in  pMOSiFETs  is  con¬ 
sidered  to  be  negligible. 

If  the  circuit  under  consideration  does  not  have  a  logic  description  available  and  its  logic 
function  cannov  be  retrieved  from  its  netlist,  then  the  delay  of  its  longest  path  will  be  con¬ 
sidered  as  its  delay.  In  that  case,  a  combination  of  topological  sort  and  worst-case  delay  pro¬ 
pagation  is  capable  of  finding  the  longest  path  and  delay  [27],  When  the  logic  function  of  the 
circuit  is  known,  then  the  sensitiz,..jility  [28]  of  the  longest  path  should  be  checked. 

Several  sensitization  criteria  of  true/false  path  identification  have  been  investigated  [28]. 
In  our  work,  we  use  static  sensitization  as  in  [29],  with  some  modification.  For  a  path 
P  =  <f  Q,f  >.  where  f  qJ internal  nodes,  to  be  statically  sensitizable,  a  tran¬ 
sition  at  fi  must  cause  a  transition  at  /^^i.  This  is  equivalent  to  saying  that  all  of  the  other 
inputs  to  the  gates  on  the  path  must  have  non-controlling  stable  values  (eg,,  i  for  AND  gates 
and  0  for  OR  gates).  The  delay  of  the  longest  statically  sensitizable  path  is  defined  as  the 
delay  of  the  circuit.  With  binary  decision  diagrams  (BDD)  [30],  [31],  we  can  efficiently  cal¬ 
culate  the  sensitizability  of  a  given  path  without  backtracking.  For  example,  in  order  for  the 
path  P  =<X,bM,€,g>  in  Fig.  3.15  to  be  statically  sensitizable,  the  condition 
(Z  =  c  =  1,  /  =0)  must  be  satisfied.  Since  c  and  /  are  not  primary  inputs,  backtracking  is 
necessary’  to  check  whether  the  logic  constraint  can  be  satisfied.  However,  if  we  generate  a 
BDD  for  each  internal  node  at  the  outset,  namely,  c  =  (X'+Y'  +  Z')  and  /  ~XZ.  The  sensit¬ 
ization  problem  can  he  transformed  into  a  satisfiability  problem  [29].  In  this  ca.se, 
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Z  nc  c^f'  =  Z  n{X'  +  Y'  +  Z')n{X'+Z')  ^  0;  therefore,  path  P  is  statically  sensitizable. 


Most  previous  work  on  critical  path  analysis  makes  the  assumption  that  each  gate  has 
only  one  delay  value  associated  with  it,  and  hence  each  path  has  one  delay.  However,  when 
the  gate  delay  is  modeled  with  macromodeling  and  RC  delay  definition,  it  is  possible  that  the 
rise  delay  of  a  gate  is  different  from  its  fall  delay,  and  hence  there  are  two  delay  values  asso¬ 
ciated  with  a  given  path.  For  example,  the  path  P  in  Fig.  3.15  has  two  delay  values. 


Pq  -f"  7/)  "4“  ITp  -f*  ITf^ 

OriMM  Lrn 


(3.34) 


and 


+  To^,  +  +  Tof^  (3.35) 

For  XOR  gates  in  which  the  output  may  switch  either  way  when  the  input  has  a  transition,  we 
use  the  larger  one  of  rise  and  fall  delay  to  make  sure  the  path  delay  is  not  underestimated. 


Figure  3,15  A  logic  circuit  example. 


58 


NANDI 


NANDI 


NAND2 


NAND3 


Figure  3.16  (a)  Probabilistic  voltage  waveforms  at  the  four  input  nodes,  A,  B,  C,  and  D. 
tfii  and  tiji  are  both  02ns  for  each  input  event; 

(b)  detailed  structure  of  a  CMOS  NAND  gate; 

(c)  original  connection  of  a  three-NANTD-gate  circuit; 

(d)  a  better  connection  of  the  three-NAND-gate  circuit  for  reliability. 


NANDI 

top 

NANDI 

bottom 

NAND2 

top 

NAND2 

bottom 

NAND3 

top 

NAND3 

bottom 

connection! 

Dr. 

P^^out 

2.223e02 

0.3 

0.12 

. . . . . 

5.927e-03 

0.08 

0.12 

3.323e-02 

0.216 

0.168 

connection2 

Dr. 

T  Phlout 


8.890e-03 

0.12 


5.927e-03 


2.584e-02 


C.08  I  0.12  I  0.168  0.216 


Table  3.4  The  damage  10**  cm~^)  o  all  the  nMOS  transistorr.  in 
the  circuit  shown  in  Fig.  3.16(c)  and  (d). 

3.7  HCE  Resistant  Design 

The  factors  that  influence  HCE  degradation  in  a  transistor  within  a  design  include  its 
feature  size,  its  switching  probability,  its  position  in  the  circuit,  the  loading  capacitance,  and 
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input  slew  rate.  Several  strategies,  based  on  these  faeiors,  have  been  previously  proposed  to 
improve  circuit  reliability.  However,  these  proposed  i  nprovements  have  been  suggested  on 
small  circuit  structures  in  isolation,  without  taking  ir.to  consideration  the  operation  of  ihe 
small  circuit  within  a  large  design,  as  we  do  here.  In  the  following,  we  first  propose  a  design 
strategy  to  reduce  HCE  effects  without  incre^.sing  area,  and  then  review  and  compare  previ 
ously  proposed  techniques. 

1.  Since  the  damage  in  a  transistor  depends  on  its  swi  ching  probability  as  well  as  its  posi¬ 
tion  in  the  circuit,  the  damage  can  be  reduced  by  connecting  nMOS  transistors  that  have 
the  least  switching  probability  directly  to  the  output  node  [24].  This  is  done  by  recon 
necting  the  circuit  input  wiring  so  that  the  top  nMOS  transistor  always  has  the  lerst 
switching  probability.  For  example.  Fig.  3.16(a)  shows  the  waveform  description  of  four 
input  nodes  A,  B,  C,  and  D.  The  circuit  under  test  is  a  three-NAND-gate  circuit  with 
detailed  gate  structure  in  Fig,  3.16(b).  When  the  circuit  is  connected  as  in  Fig.  3.16(c), 
the  HCE  degradation  is  listed  in  row  1  of  Table  3.4,  as  well  as  the  switching  probability 
of  each  nMOS  transistor.  It  is  noticed  that  the  bottom  transistors  of  gate  NANDI  and 
NAND3  have  lower  switching  probability  value.  After  reconnecting  the  input  of  these 
two  gates  as  in  Fig.  3.16(d),  the  HCE  degradation  is  reduced,  as  shown  in  the  second 
row  of  Table  3.4.  The  advantage  of  this  method  is  that  the  circuit  area  is  not  increased. 
However,  this  method  of  reordering  the  transistors  cannot  be  applied  to  NOR  gates  in 
which  all  the  nMOS  transistors  are  directly  connected  to  the  output.  In  this  case,  some 
other  method,  such  as  increasing  the  area  of  the  affected  transislcr  or  applying  method  3 
below  could  tie  used,  with  some  sacrifice  in  performance.  A  better  approach  would  be  to 
eliminate  NOR  gates  with  high  switching  activity  through  functional  transformations. 


2.  By  adding  a  self-bootstrapping  nMOS  transistor,  as  shown  in  Fig.  3.17(a),  the  peak  sub¬ 
strate  current  in  transistor  NMl  can  be  reduced  by  as  much  as  3.9  times  according  to 
[32]  and  the  interface-state  generation  can  be  suppressed  by  one  to  two  order  of  magni¬ 
tude.  In  addition,  this  structure  also  has  shorter  rise  and  fall  time  compared  to  conven¬ 
tional  NAND  logic  circuit.  A  disadvantage  of  this  structure  is  an  increase  in  circuit  area; 
therefore,  it  is  preferred  to  add  a  small  self-bootstrapping  nMOS  transistor.  In  addition, 
this  structure  cannot  be  used  to  improve  NOR  gates. 


Figure  3.17  Several  HCE  resistant  circuit  design  strategies. 
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3.  Another  HCE  suppressing  circuit  redesign  technique,  shown  in  Fig.  3.17(b),  has  been 
proposed  in  [33],  where  a  normally-on  nMOS  transistor  is  added  between  the  n-part  of  a 
CMOS  circuit  and  the  output  node.  Since  the  added  transistor  is  always  on,  the  logic 
function  of  the  circuit  is  not  affected.  In  addition,  the  added  transistor  is  always  operat¬ 
ing  in  the  linear  region  and  hence  has  negligible  hot-carrier  generation.  The  rest  of  the 
nMOS  transistors  are  not  directly  connected  to  the  output  node,  and  thus  are  relatively 
safe  from  HCE  degradation.  This  technique  can  also  be  applied  to  NOR  gates,  as  shown 
in  Fig.  3.17(c).  Its  disadvantage  is  the  increase  in  layout  area;  therefore,  it  should  be 
selectively  used. 

4.  In  Clocked  CMOS  (C^MOS)  circuits  where  the  clock  signal  arrives  after  all  the  inputs 
settle,  we  recommend  that  the  control  nMOS  transistors  be  put  anywhere  but  the  top. 
With  this  restriction,  the  top  nMOS  transistor(s)  never  switches  last,  i.e.,  their  switching 
doe.s  not  cause  output  transitions,  and  hence  HCE  degradation  can  be  prevented.  Fig. 

1 7(d)  illustrates  a  HCE  resistant  C^MOS  design.  The  same  strategy  can  be  applied  to 
dynamic  CMOS  design  as  well.  However,  the  effects  of  such  a  design  strategy  on  cir¬ 
cuit  delay  and  performance  need;;  to  be  investigated. 

A  CMOS  circuit  design,  shown  in  Fig.  3.18,  has  been  proposed  for  testability  in  [.34], 
where  it  is  proposed  that,  with  the  four  combinations  (00,  01,  10,  1 1)  of  the  control  sig¬ 
nals  <|),  and  4)2,  all  stuck-open  and  stuck-on  faults  of  the  circuit  can  be  detected.  The  cir¬ 
cuit  operates  as  follows;  (1)  When  4>i  =0,4)2=  s'ircuit  operates  normally  and 

implements  the  same  function  as  the  one  without  the  transistors  and  (shorted).  (?.) 
When  4>i  =  1 .4>2  -  d,  both  and  art.  off.  and  the  output  state  does  not  change  even 
if  the  input  signals  have  transitions.  (.3)  When  (j),  =  0,4)2  ~  d,  only  P ^  is  on  and  all  the 
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Figure  3.18  A  CMOS  circuit  sthjcture  for  testability, 
stuck-open  and  stuck-on  faults  of  the  pMOSFET  network  can  be  detected.  (4)  Finally, 
when  (j)j  =  1,<|)2  =  1,  only  is  on  and  all  the  stuck-open  and  stuck-on  faults  of  the 
nMOSFET  network  can  be  detected.  We  found  that  this  structure  can  be  employed  for 
improving  reliability  as  well.  There  are  two  situations;  (1)  If  fbe  circuit  is  in  normal 
operation  most  of  the  time,  i.e.,  the  control  signals  rarely  switch,  then  this  design  is  simi¬ 
lar  to  the  structure  in  Fig.  3.17(b).  (2)  If  the  system  clock  is  used  as  the  control  signal, 
then  we  need  to  modify  the  design  by  moving  to  the  bottom  of  the  nMOSFET  net¬ 
work,  similar  to  the  c  isc  of  C^MOS  design. 

The  algorithms  and  strategies  above  have  been  included  in  iProbe  d  and  used  to  simulate 
and  redesign  some  benchmark  circuits.  The  results  of  simulating  the  ISCAS85  benchmark 
circuits  are  shown  in  Table  3.3.  After  obtaining  the  damage  and  switching  probability  of  each 
transistor,  the  critical  path  analysis  algorithm  is  employed  to  check  the  global  degradation  cri¬ 
terion.  Table  3.5  shows  the  circuit  delay  before  and  after  HCE  dv-tgradation.  The  delay 
increase  mostly  comes  from  INVERTERS  and  NOR-type  gates;  all  of  the  NAND-type  gates 
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with  more  than  two  inpuUi  have  less  than  0.05ns  delay  increase.  TTiis  observation  matche,s  our 
suggestion  that  NAND-type  gates  are  more  HCE  resistant  than  NOR-tyiie  gates.  Table  3.6 
shows  delay  degradation  for  the  same  circuits  when  Vqd  -  5.5v .  Dote  that  by  increasing 
the  cimuit  becomes  faster,  but  the  dai..age  is  increased. 


Table  3.5  The  circuit  delay  (icngest  path)  of  several  benchmaik  circuits 
before  and  after  HCE  degradation.  ^DD  -  •  Delay  in  nsec. 


Circuit 

No. 

No 

Delay 

Delay 

Delay 

Gates 

Nodes 

before  HCE 

after  HCE 

increase 

432 

160 

198 

20.88 

20.95 

0.34% 

C499 

202 

245 

16.46 

16.50 

0.24% 

C880 

383 

445 

19.02 

19.05 

0.16% 

C1355 

546 

589 

20.11 

20.15 

0.20% 

Cl  908 

880 

915 

27.71 

27.98 

0.97% 

C2670 

1193 

1428 

31.38 

31.57 

0,61% 

C3540 

1669 

1721 

37.87 

38.13 

0.69% 

C5315 

2307 

2487 

34.69 

34.98 

0.84% 

C6288 

2H0 

2450 

80.04 

80.93 

1.09% 

C7552 

3512 

3721 

28.78 

29.01 

0.80% 

1  Table  3.6  The  circuit  delay  (longest  path)  of  several  benchmark  circuits  I 

before  and  after  HCE  degradation  when  Vqq 

=  5.5v.  Delay 

in  nsec. 

Circuit 

No. 

No 

Delay 

Delay 

Delay 

Gates 

Nodes 

before  HCE 

after  HCE 

increase 

432 

160 

198 

1798 

18.21 

1.28% 

C499 

202 

245 

14.18 

14.30 

0.85% 

C880 

383 

445 

16.40 

16.61 

1.28% 

C1355 

546 

589 

17.32 

17.49 

0.98% 

Cl  908 

880 

915 

23.90 

25.54 

6.86% 

C2670 

1193 

1428 

27.07 

27.84 

2.84% 

C3540 

1669 

1721 

32.67 

33.89 

3.73% 

C53i5 

2307 

2487 

29.93 

31.13 

4.01% 

C6288 

2416 

2450 

69.19 

77,81 

12.46% 

C7552 

3512 

3721 

24.83 

25.69 

3.46% 
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Preliminary  results  of  redesign  improvements  are  shown  in  Table  3.7.  'fhis  table  shows 
the  impact  of  three  strategies,  namely,  (1)  increasing  transistor  width,  (2)  adding  normally-on 
transistors,  and  (3)  interchanging  gate  input  lines  on  benchmark  circuit  C3540  with 
We  set  the  criteria  of  redesign  to  be  normalized  equal  to  5.0  for  individual 
transistor  damage  and  0.05ns  for  subcircuit  delay  increase.  We  also  assume  the  equivalent 
resistance  of  the  added  normaily-on  transistor  to  be  one  fourth  that  of  the  other  nMOSFETs. 
As  can  be  seen  from  this  table,  increasing  transistor  widths  can  improve  circuit  delay  but  have 
less  improvement  on  transistor  damage.  Adding  normally-on  transistors,  on  the  other  hand, 
does  eliminate  HCE  degradation  in  all  the  critical  transistors,  while  the  circuit  delay  is  slightly 
improved  compared  to  the  damaged  circuit.  Both  strategies  increase  the  total  circuit  area  by 
20%,  assuming  the  normally-on  transistor  occupies  same  area  as  one  regular  nMOSFET.  This 
increase  may  not  be  acceptable,  and  other  redesign  strategies  that  do  not  impact  the  area  need 
to  be  considered.  The  impact  of  interchanging  input  lines  on  HCE  is  shown  in  the  third  row 
in  the  table.  Note  that  for  this  example  the  improvements  in  damage  and  delay  resulting  from 
interchanging  input  lines  at  critical  NAND  gates  along  longest  paths  are  minimal.  The  reason 
is  that  in  this  example  the  switching  probabilities  at  the  inputs  of  critical  NAND  gates  are  not 
too  different,  and  interchanging  the  connections  dees  not  improve  reliability  in  a  noticeable 
way.  This  study  indicates  that  design  for  reliability  should  be  done  as  part  of  the  design  syn¬ 
thesis  process  rather  than  afterwards. 
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Table  3.7  (1)  Increasing  the  width  of  critical  transistors  (A/;.  >  5.0)  by  75%; 
(2)  Adding  a  normally-on  transistor  on  top  of  n-logic  of  critical  subcircjits; 
(3)  Interchanging  gate  input  lines.  ~  5.5V. 

Circuit  delay  before  KCE  was  32.67  nsec. 


C3540 

No 

<1.0 

of  nMOS 
1.0- 2.0 

S  transisto 
2.0-3.0 

rs  in  each 
3.0-4.() 

damage  le 
4.0-5. 0 

VC  I 
>5.0 

Circuit  delay 
after  HCE” 

initial 

950 

177 

_ 1 

142 

is 

~2i67~ 

33.89 

(1) 

950 

magm 

142 

177 

1846 

32.84 

msm 

1^1 

13S 

■on 

Bi 

0 

33.26 

(3) 

967 

199 

150 

112 

215 

2109 

33.87 
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4.  MAXIMUM  CURRENT  AND  VOLTAliJE  DROP  ESTIMAIlTON 


4.1  InlrtMlucdon 

Voltage  drop  is  another  imporrant.  reliability  issue  which  is  related  to  the  current  flow  in 
the  power  bus.  Excessive  supply  cuiTcnts  can  scvtrelj  affect  both  circuit  reliability  and  per¬ 
formance  by  causing  excessive  voltage  drops  in  the  Power  ana  Ground  (P&G)  lines.  Exces¬ 
sive  voltage  drops  manifest  themselves  as  glitches  on  tlie  P&G  lines,  and  cause  eiToneous 
logic  signals  (soft  errors)  and  degradation  in  switching  speeds.  Severity  of  the  ‘'oltage  drop 
problems  intensify  with  the  continuing  push  for  denser  chips  and  finer  technologies.  With 
higher  currents  flowing  in  narrower  lines,  the  voltage  drops  in  the  PG  lines  go  up  and  quickly 
become  a  limiting  factor  in  the  de.sign  of  VLSI  chips.  Furthermore,  a  lower  supply  voltage 
means  that  the  noise  margins  for  the  correct  operation  of  the  transistors  on  the  chip  decrease. 
In  order  to  avoid  logic  errors,  the  circuit  needs  to  be  appropriately  designed  to  take  care  of 
increased  voltage  drops  and  reduced  noise  margins.  Worst  case  voltage  drops  are  determined 
by  maximum  cument  Tows  ratler  than  averages.  Thus  this  research  subtask  is  focused  on  the 
problem  of  estimating  maximum  currents  in  the  P&G  lines, 

Current  drawn  by  a  CMC  S  circuit  depends  upon  the  specific  input  pattern  applied  at  its 
inputs,  .^n  inpiii  pattern  for  a  circuit  with  n  inputs  is  defined  as  a  vector  of  n  excitations, 
where  each  excitation  could  oe  any  one  of  four  possibilities  i.e.,  low.  high,  high  low  or 
low  to  high  .  For  different  input  patterns,  different  transient  current  waveforms  are  drawn  at 
the  contact  poirts.  Therefore,  m  the  presence  of  such  input  dependent  and  transient  current 
waveforms,  vve  need  to  define  'vhat  we  mean  by  the  maximum  current  waveform  ai  a  contact 
point. 
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Accurate  estimation  of  the  maximum  current  waveform  at  every  contact  point  is 
extremely  difficult  since  for  that  we  need  to  determine  current  waveforms  corresponding  to  all 
possible  input  patterns.  If  the  circuit  has  n  primaiy  inputs  then  we  need  to  explore  the  set  of 
2^"  input  patterns  to  calculate  the  maxiinum  current  waveforms.  This  makes  the  problem 
practically  impossible  to  handle  by  any  of  the  known  search  procedures  for  large  circuits. 

Several  papers  (such  as  [36]-[39])  have  appeared  in  literature  on  the  estimation  of  P&G 
currents  from  deterministic  input  patterns  These  methods  offer  significant  improvement  in 
execution  times  compared  to  SPICE,  while  providing  acceptable  accuracy  in  the  current 
wavefonns.  These  methods  can  be  used  for  finding  maximum  currents  for  small  circuits  hav¬ 
ing  a  few  inputs,  by  calculating  the  current  v/aveforms  corresponding  to  all  possible  input  pat¬ 
terns.  However  they  are  not  much  helpful  for  large  circuits,  as  they  do  not  guide  us  in  select¬ 
ing  an  inpL  t  pattern  that  leads  to  the  maximum  currents. 

In  [35],  Chowdhury  et  ai.,  have  addressed  the  problem  of  maxiinum  current  estimation 
for  iarge  circuits.  The  circuit  is  first  divided  into  a  set  of  macros,  where  each  macro  consists 
of  a  combinational  interconnection  of  logic  gates.  Considering  each  macro  separately,  either 
an  exact  search  technique  (namely,  branch  and  bound)  or  a  heuristic  technique  is  used  to  find 
the  maximum  of  its  transient  currents,  assuming  that  the  macro  has  only  one  contact  point  and 
its  inputs  switch  simultaneously  In  the  analysis  of  the  bus,  to  calculate  maximum  voltage 
drops,  the  authors  assume  that  all  the  macros  draw  their  maximum  currents  simultaneou.sly. 
Because  of  these  assumptions,  their  methodology  overestimates  the  worst  case  cunents  and 
voltage  drops.  Secondly,  due  to  the  huge  size  of  the  input  space,  theii  brancti  and  hound 
se;irch  technique  is  slow  on  iarge  circuits.  Furthermore,  their  heuristic  approach  does  not 
guarantee  an  upper  bound  on  the  maximum  currents. 
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In  [40],  Devadas  et  al.,  have  addressed  a  similar  problem.  They  consider  the  estimation 
of  worst  case  power  dissipation  in  CMOS  combinational  circuits.  They  reduce  this  problem 
to  a  weighted  max -satisfiability  problem,  on  a  set  of  multi-output  Boolean  functions  obtained 
from  the  circuit  logic  description.  These  functions  are  appropriately  weighted  to  account  for 
different  load  capacitances.  They  then  use  either  a  disjoint  cover  enumeration  algorithm  or 
the  branch  and  bound  algorithm  to  solve  the  (/^P -complete,  max-satisfiability  problem.  They 
are  able  to  account  for  the  glitches  at  various  internal  nodes.  However,  for  a  multilevel  logic 
circuit,  even  under  a  unit  gate  delay  assumption,  the  functions  generated  by  their  algorithm 
are  fairly  complex.  Consequently,  even  for  small  circuits,  their  analysis  is  slow.  Analysis  of 
multi-level  circuits  under  a  general  delay  model  was  not  attempted. 

In  this  chapter,  we  propose  a  better  measure  of  the  maximum  current  waveform  at  a  con¬ 
tact  point  called  maximum  envelope  current  (MEC)  waveform.  We  then  propose  a  pattern 
independent,  linear  time  (in  the  number  of  gates)  algorithm  (iMax)  that  provides  tight  upper 
bounds  on  the  MEC  waveforms.  The  proposed  approach  represents  a  trade-off  between  exe¬ 
cution  speed  and  tightness  of  these  bounds.  In  order  to  maintain  reasonable  execution  times, 
the  iMax  algorithm  neglects  various  signal  correlations  that  exist  inside  a  circuit.  As  will  be 
shown  later,  while  in  most  cases  iMax  produces  good  upper  bound  waveforms,  in  some  cases 
the  loss  due  to  signal  correlatioi.a  can  be  significant.  We  then  propose  a  new  partial  input 
enumeration  (PIE)  algorithm  that  efficiently  resolves  these  correlations  and  leads  to  significant 
improvement  in  the  upper  bound  waveforms.  The  PIE  algorithm  is  based  on  (1)  intelligently 
.selecting  a  few  critical  inputs  and  (2)  enumerating  a  limited  number  of  cases  at  these  inputs 
to  produce  an  overall  improvement  in  the  upper  bound  waveform  at  every  contact  point.  It 
turns  out  that  the  choice  of  these  c  itical  inputs  is  the  key  to  improving  the  upper  bound.  Wc 
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pre5v-nt  two  heuristics  for  automatically  selecting  the  critical  inputs,  that  have  shown  good 
results  in  practice.  While  the  PIE  algorithm  is  slower  than  the  simple  iMax  algorithm,  we 
demonstrate  good  speed  and  accuracy  performance  results  on  circuits  with  over  twenty 
thousand  gates.  Furthermore,  the  algorithm  has  the  attractive  property  that  it  does  an  iterative 
improvement,  so  that  one  can  stop  the  algorithm  at  any  time  and  obtain  a  better  upper  bound 
than  the  simple  iMax  result. 

In  the  next  section,  we  discuss  various  assumptions  on  which  our  algorithms  are  based 
and  describe  the  proposed  maximum  current  estimate.  After  that,  we  present  the  iMax  algo¬ 
rithm  in  detail  in  Section  4.3.  Experimental  results  on  several  benchmark  circuits  using  iMax 
are  also  provided  in  this  section.  The  signal  correlation  problem  is  described  in  Section  4.4. 
This  ib  followed  by  a  discussion  of  possible  methods  that  can  be  used  to  resolve  the  signal 
correlations  in  Section  4.5.  in  Section  4.6,  we  present  the  partial  input  enumeration  algorithm 
along  with  extensive  experimental  results  on  several  benchmark  circuits. 

4.2  Maximum  Current  Estimates 

We  consider  sync!  ronous  digital  circuits  consisting  of  combinational  blocks  separated  by 
latches  (Fig.  4.1)  such  that  all  the  inputs  to  each  combinational  bloc:k  switch  simultaneously. 
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Figure  4.1  An  example  of  a  latch  controlled  synchronous  digit;  i  ci  ;u  t. 
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when  they  do.  This  effectively  eliminates  the  time  domain  uncertainty  about  the  input  transi¬ 
tions.  This  assumption  has  also  been  used  by  all  the  previous  approaches. 

We  assume  that  the  delay  of  each  gate  in  the  circuit  is  specified.  Different  gates  can 
have  different  delays.  Further,  we  assume  that  each  time  the  output  of  a  gate  switches,  a  pulse 
of  current  approximated  by  a  triangular  waveform  is  drawn  from  the  supply  lines,  as  shown  in 
Fig.  4.2.  The  duration  of  this  pulse  is  computed  from  the  delay  of  the  gate  and  its  height 
(peak  value)  is  pre-characterized.  Two  separate  values  for  the  peak  current  may  be  .^.pecified, 
corresponding  to  low  to  high  and  high  to  low  transitions  at  the  gate  output. 

Given  the  specific  clocking  scheme  of  the  synchronous  circuit,  the  maximum  current 
waveforms  from  different  combinational  blocks  can  be  appropriately  shifted  in  time  depending 
upon  the  individual  clock  trigger.  We  will  thus  focus  on  the  analysis  of  a  single  combina¬ 
tional  block  whose  inputs  switch  at  time  zero. 

We  define  excitation  at  a  node  (or  net)  at  any  time  t  as  the  stimulus  (or  signal  value) 
present  at  the  node  at  that  time.  At  any  time,  a  node  in  the  circuit  could  be  either  stable  at 
low  or  high ,  or  could  transition  from  high  to  low  or  from  low  to  high .  Thus,  the  excitation 
could  be  any  single  value  from  the  set  X  =  / ,  /i ,  /z  / ,  /  /i ,  where  I  =  low ,  h  =  high ,  h  i  - 

User  specified  Peak 


Figure  4.2  Model  of  a  gate  current  pulse. 
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Figure  4.3  The  Maximum  Envelope  Current  (MEC)  wav  sform. 
high  to  low  transition  and  //i  =  low  to  high  transition. 

For  each  input  pattern  that  is  applied  to  the  circuit,  a  transient  curr  ;nt  waveform  is  drawn 
from  the  P&G  lines  at  a  contact  point.  As  shown  in  Figure  4.3,  inste  id  of  representing  the 
maximum  current  at  a  contact  point  by  a  single  dc  value,  in  our  approi  i;h,  v/e  represent  it  by 
a  waveform  whose  value  at  any  time  is  the  maximum  current  value  that  the  circuit  can  draw 
at  that  time.  We  call  this  the  Maximum  Envelope  Current  (MDEC)  waveform. 

Suppose  that  a  circuit  under  consideration  has  n  inputs  a^id  when  an  input  pattern 
p  =  (^1,^2,  ■  •  •  .e„),  where  e,  6  X,l<  i  <n,  is  applied  to  che  circuit,  a  transient  current 
waveform  /^(r)  is  drawn  from  the  circuit  at  the  contact  point.  Denote  the  set  of  all  possible 
input  patterns  that  can  be  applied  to  Lhe  circuit  hy  U  -  {(t  -  )  1  e/  e  X ,  1  <  /  <  n  }. 

If  the  value  of  the  MEC  waveform  at  any  tune  t  at  the  cor  tact  point  is  denoted  by  /mec(0. 
then  the  following  equation  applie.s 

~  )  (4.1) 

pe.U 

Clearly,  the  MEC  waveform  is  the  maximum  envelope  c'  all  transient  current  wavefonns 
corresponding  to  all  possible  input  patterns.  There  is  a  aniquc  MEC  waveform  at  eveiy 


contact  point 


If  the  power  or  the  ground  bi’3  of  a  circuit  is  represented  by  an  equivalent  RC  network, 
then  we  have  the  following  result: 

f  -r  ) 

Theorem  1:  The  voltage  drop  ,  calculated  at  any  node  k  in  the  power  or 

ground  bus  when  the  MEC  waveforms  are  applied  at  the  contact  points,  is  an  upper  bound  on 
the  voltage  drop  [^jf(^)]  occurring  at  that  node  when  any  input  pattern  p  is  applied  to  the  cir¬ 
cuit,  i.e.,  >  Vj^(t),  for  all  t. 

4.3  The  Proposed  (iMax)  Algorithm 

The  proposed  pattern  independent,  linear  time  a'gorithm  operates  at  the  gate  level 
description  of  the  circuit.  Unless  specified  by  the  user,  we  assume  that  nothing  is  known  about 
the  specific  excitations  at  the  primary  inputs,  except  that  they  may  transition  (only)  at  time 
zero,  i.e.,  each  primary  input  may  cany  any  excitation  from  the  set  X  at  time  zero.  We  call 
this  an  uncertainty  about  these  input  S’  ;nals.  The  basic  idea  of  the  proposed  algorithm  is  to 
propagate  this  uncertainty  present  at  the  inputs  inside  the  circuit,  so  that,  at  the  output  of 
every  logic  gate,  we  know  the  set  of  ail  possible  excitations  and  their  associated  timing.  From 
this,  the  worst  case  gate  currents  are  computed,  as  explained  below. 

4.3.1  Signal  RepresenU.tion 

Perhaps  the  first  question  that  comes  to  mind  is  what  kind  of  information  one  maintains 
in  order  to  represent  the  signal  uncertainty  about  internal  circuit  nodes.  Ideally,  one  would 
like  to  compute  the  set  ot  a'l  possible  transitions  (along  with  their  liming  inforrnalion)  that 
occur  at  the  output  of  every  gate  in  the  circuit.  That  would  certainly  be  enough  to  estimate 
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the  MEC  waveforms.  However,  due  to  the  uncertainty  at  the  primary  inputs  and  the  general 
gi  ie  delay  model  used,  the  number  of  possible  transitions  at  internal  gates  grows  exponen- 
tia'Iy,  and  quickly  becomes  a  bottleneck.  To  avoid  this  problem,  we  maintain  information,  not 
about  individual  transitions,  but  about  intervals  during  which  the  output  of  the  gate  might 
switch.  Thus,  for  each  of  the  excitations  low,  high,  hi  and  //j,  we  store  a  list  of  intervals 
during  which  a  node  might  carry  that  excitation.  These  intervals,  which  might  overlap,  serve 
to  describe  the  signal  uncertainty.  We  call  these  intervals  uncertainty  intervals. 

Definition  1  (Uncertainty  Set  X„(/)):  The  uncertainty  set  at  time  t  for  a  node  n  defines 
the  set  of  all  excitations  that  the  node  can  possibly  assume  at  that  time.  X^(t)  £  X . 

Definition  2  (Uncertainty  Waveform):  The  uncertainty  v  aveform  describes  the  signal 
uncertainty  present  at  a  node  as  a  function  of  time.  At  time  t ,  the  set  of  values  taken  by  the 
waveform  is  the  uncertainty  .set  for  the  node  at  that  time. 

An  example  of  the  uncertainty  waveform  is  given  in  Fig.  4.4.  In  this  figure,  we  show  an 
uncertain  signal  U{t)  represented  as  four  sets  of  intervals'  ttlong  the  time  axis.  Thus,  if  u{t) 
is  a  logic  signal  that  belongs  to  the  family  U{t),  i.e.,  u{t)  e  U(t),  then  u(t)  will  be  low  up 
to  t|,  wiil  switch  from  low  to  high  sometime  between  tj  and  t,*  "'•il  ^hen  be  high  up  to 
etc.  Thus,  at  any  time  between  f|  and  t2,  the  signal  can  be  either  high  or  low  and  may  switch 
any  number  of  times  between  these  values,  but  will  be  high  at  Notice  that  between  tj,  and 
l-i  the  signal  may  make  any  number  of  low  to  high  and/or  high  to  low  transitions.  At  the 
primary  inputs,  signals  are  represented  by  such  waveforms  with  a  single  point  of  possible 
transition  at  time  0.  As  internal  signals  are  generated,  the  number  of  possible  transition  points 

'Ont‘  sei  if  tniervais  tor  each  /oh.  hi  and  Ih 
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increases.  In  order  to  contain  the  complexity,  we  then  start  to  merge  neighboring  transition 
points  into  intervals.  In  general,  this  strategy  can  be  stated  as  follows:  when  the  number  of 
intervals  associated  with  a  gate  corresponding  to  any  excitation  exceeds  a  certain  user- 
specified  threshold  (Max_No_Hops),  we  repeatedly  merge  closest-neighbor  intervals,  so  as  to 
keep  their  count  below  the  threshold. 

43.2  Independence  Assumption 

While  propagating  infonnation  at  a  logic  gate,  we  know  the  uncertainty  waveforms  at 
each  of  its  inputs  and  we  would  like  to  derive  the  corresponding  waveform  at  its  output. 
However,  one  cannot  do  this  accurately  without  knowing  how  some  of  these  inputs,  if  any, 
aie  correlated.  For  instance,  certain  combinations  of  the  gate  input  excitations  may  not  be 
possible.  Unfortunately,  maintaining  infonnation  about  correlation  between  various  circuit 
nodes  is  very  expensive.  We,  therefore,  use  a  conservative  approximation,  one  that  does  not 
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Figure  4.4  An  illu.stration  of  uncertainty  waveform  at  a  gate  output, 
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underestimate  the  MEC  waveforms,  as  follows.  If  we  assume  that  all  combinations  of  the 


gate  input  excitations  are  possible,  i.e.,  the  gate  inputs  are  independent,  then  the  worst  case 
current  in  that  case  will  be  an  upper-bound  on  the  gate  current  for  the  case  when  the  inputs 
are  dependent.  In  other  words,  the  worst  case  current  over  all  combinations  of  inputs  is  cer¬ 
tainly  an  upper  bound  on  the  worst  case  current  over  some. 

4.3.3  Single  Gate  Simulation 

Given  the  type  of  a  Boolean  gate  and  the  independence  assumption  for  the  uncertainty 
waveforms  at  its  inputs,  we  now  describe  how  the  uncertainty  waveform  at  the  output  of  the 
gate  is  calculated.  This  process  is  divided  into  the  following  two  parts: 

1 .  Calculation  of  the  uncertainty  .set  at  the  output  of  the  gate  at  a  time  t . 

2.  Calculation  of  uncertainty  intervals  at  the  output  of  the  gate. 

4.3.3.1  Calculating  Uncertainty  Sets 

One  can  calculate  ail  possible  excitations  at  the  output  of  the  gate  at  time  /  from  the 
excitations  present  at  the  inputs  at  time  t  -D,  where  D  is  the  delay  of  the  gate.  Let  us 
denote  the  uncertainty  set  at  the  i‘^  input  of  the  gate  at  time  t  ~  D  by  X, ,  Let  us  further  sup¬ 
pose  that  the  gate  lias  m  inputs.  Then  the  set  of  all  possible  input  patterns  that  lead  to  an 
excitation  at  the  output  of  the  gate  at  time  t  can  be  represented  by 
{(x  |,a:2,  ••  •  )  lx,  e  X; ,  1  <  1  <  m\.  For  each  input  pattern,  the  output  of  the  gate  can  be 

easily  determined  from  the  Boolean  equation  of  the  gate.  Thus,  by  calculating  the  output  of 
(he  gate  corresponding  at  each  input  pattern,  the  uncertainty  set  at  the  output  of  the  gate  at 
time  t  can  be  determined.  This  proce.ss  would  require  one  to  generate  and  eviluate 
iXiHXjf  •  ■  I.Y,„  I  input  patterns.  This  worst  case  complexity  can  be  greatly  reduced  by  the 
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following  observations. 


1.  The  above  input  pattern  generation  and  evaluation  process  can  be  stopped  when  the 
uncertainty  set  at  the  output  of  the  gate  becomes  equal  to  X.  Obviously,  trying  out  any 
more  input  patterns  would  not  lead  to  any  improvement  in  the  uncertainty  se^  at  the  out¬ 
put  of  the  gate. 

2.  An  input  to  a  gate  is  completely  ambiguous  at  time  t  if  its  uncertainty  set  at  time  t  is  X. 
If  all  the  inputs  to  the  gate  are  completely  ambiguous  at  time  t  -  D,  then  the  output  of 
the  gate  is  also  completely  ambiguous  at  time  i . 

3.  All  the  gates  in  a  circuit  can  be  divided  into  the  following  two  categories. 

(a)  Gates  who.se  outputs  depend  not  only  on  the  excitations  present  on  its 

inputs  but  also  on  the  total  count  of  input  lines,  e.g.,  XOR,  XNOR  etc., 
gates. 

(b)  Gates  whose  outputs  do  not  depend  on  the  count  of  the  inputs.  The  out¬ 
puts  of  such  gates  depend  only  on  the  specilic  excitations  pre.sent  on  the 
inputs,  e.g.,  NAND,  NOR  etc.,  gates 

For  all  the  gates  which  belong  to  the  second  category  above,  all  the  input  lines  which 
have  the  same  uncertainty  sets  can  be  merged  into  a  single  line.  This  observation  leads 
to  a  reduction  in  the  numbc’  ^  '..put  lines  and  thus  the  number  and  size  of  the  input  pat¬ 
terns. 

All  of  these  observations  lead  to  tremendous  savings  in  the  calculation  of  uncertainty  .sets  at 

the  (.lutput  of  the  gate  and  thus  contribute  to  the  speed  ot  the  algiaithni. 
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43.3,2  Cai^'jilating  Uncertainty  Intervals 


In  iMax,  sihce  .signals  aie  represented  in  the  fonti  ot  uncertainty  intervals  at  the  inputs  of 
the  gate,  therefore,  the  output  of  the  gate  would  rJfo  be  in  the  form  of  uncertainty  intervals. 
An  interv'al  at  the  output  o.f  the  gate  could  begin  or  end  at  time  t  only  if  an  interval  begins  or 
ends  at  any  of  the  inputs  at  time  t  ~  D.  Between  the  time  at  inputs  when  an  interval  begins 
or  ends,  and  th,.  next  interval  begins  or  ends,  the  sets  of  excitations  that  the  inputs  can  assume 
do  not  change  and  therefore  no  corresponding  uncertainty  interval  can  begin  or  end  at  the  out¬ 
put  during  that  time  shifted  by  D .  Based  upon  these  observations,  the  uncertainty  intervals  at 
the  output  oi  the  gate  can  oe  easily  calculated. 

An  example  illustrating  how  uncertainty  waveforms  at  various  circuit  nodes  are  calcu¬ 
lated  is  shown  in  Fig.  4.5. 


Input  Description  :  il,  i2  G  {l,h,hl,lh}  at  time  0. 

Uncertai,nty  Intervals  : 

il,  i2:  0],  hl[0,  Oj,  1(0,  oc),  MO,  oo) 

ul:  1],  IJ,  1(0,  oo),  /i(0,‘oo) 

ol:  /A[2,  2][3,  3),  Mt2,  2|[3,  3],  1(0,  c»),  /i[0,  oo) 
if  HAXJJD JiOPS  =  1  then 

ol:  l/i[2,  3j,  /M[2,  3),  /[O,  oo),  li[0,  x) 

Key  :  Exdtatiori[Interva)  Begin,  Interval  End] 

Figure  4.5  .An  cxtunple  of  uncertainty  waveforms  calcuiation. 
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Figure  4.6  Current  wavefcrm  due  to  an  uncertainty  inten'al. 

43.4  Current  Calculation 

After  the  uncertainty  waveform  at  the  output  of  a  gate  is  laiowri,  its  current  contribution 
is  caJculated  next.  Since  the  output  of  the  gate  could  switch  at  any  time  during  an  uncertainty 
interval,  therefore,  a  triangular  pulse  of  current  could  be  drawn  at  any  time  during  that  inter¬ 
val  (shifted  backwards  by  the  delay  of  the  gate)  from  the  P&G  lines  due  to  these  transitions, 
as  shown  in  Fig.  4.6.  Hence,  by  taking  an  envelope  of  all  possible  triangular  cun-ent  pulses, 
we  get  the  worst  case  current  contribution  of  a  gate  for  an  uncertainty  interval.  At  every  gate, 
there  are  two  types  of  uncertainty  intervals  that  result  in  some  switching  activity  at  the  output 
and  therefore,  there  are  two  possible  current  waveforms,  one  due  to  the  hi  uncenainty  inter¬ 
vals,  called  hlCurrent  and  the  other  due  to  Ik  uncertainty'  intervals,  called  IhCiirrent.  Since 
at  any  time,  the  output  of  the  gate  could  switch  either  from  high  to  low  or  from  low  to  high , 
therefore,  by  taking  an  envelope  of  the  hlCurrent  and  Ihtiiurrent  waveforms,  we  get  the 
maximum  current  contribution  of  the  gate.  Once  all  the  gate  currents  arc  calculated,  the 
current  waveforms  at  the  contact  points  are  calculated  by  combining  the  individual  current ;  of 
those  gates  that  are  tied  to  n. 
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4.3.5  Implementation  Details 


The  above  approach  has  been  implemented  in  a  program  in  C.  In  the  program,  the  cir¬ 
cuit  is  first  leveliied  so  that  the  output  of  a  gate  at  level  j  does  net  feed  any  other  gate  at  a 
level  less  than  or  equal  to  j.  Any  user-specified  restrictions  on  certain  inputs  are  then 
imposed,  while  all  other  inputs  are  assumed  to  take  all  possible  excitations  from  the  set  X. 

this,  the  circuit  is  analyzed  in  a  leve'  by  levc  fashion,  starting  from  the  lowest  level,  by 
propagating  the  uncertainty  waveforms  at  the  inputs  of  every  gate  to  its  output.  As  a  result, 
we  get  a  current  waveform  that  is  a  point-wise  upper  bound  on  the  MEC  waveform  at  every 
contact  point. 

In  order  to  assess  the  quality  of  the  solution  at  every  contact  point,  we  need  to  deter¬ 
mine,  how  close  the  upper  bound  waveforms  obtained  from  iMax  are  to  the  exact  MEC 
waveforms.  One  way  of  doing  this  would  be  to  perform  an  exhaustive  enumeration  over  all 
possible  input  patterns  and  actually  calculate  the  MEC  waveforms  at  every  contact  point. 
However,  this  would  be  very  expensive  <md  practically  impossible  tor  circuits  with  more  than 
about  10  inputs  (Note:  4*®  =  1,048,576).  We,  therefore,,  resort  to  the  following  random  optim¬ 
ization  approach.  We  repeatedly  apply  different  input  patterns,  randomly  selected  from  the  set 
of  all  possible  patterns,  to  the  circuit.  Then  we  use  a  logic  simulator  to  calcu.nuo  the  outputs 
of  various  gates.  From  these  gate  outputs,  the  P&G  cuiTent  wavefonns  at  every  contact  point 
are  easily  calculated.  At  every  contact  point,  by  maintaining  an  upper-bound  envelope  of  the 
current  w’aveforms  obtained  for  different  input  patterns,  we  basit:ally  get  a  low'er  bound  on  the 
MEC,'  vvavcfoim.  Naturally,  ihc  more  patterns  are  simulated  the.  closer  this  waveform  will  get 
In  the  MEC.’  waveform.  Ideally,  one  wmuld  like  to  see  the  upper  bound  obtained  iVorn  iMa.x 
come  as  clt-'se  to  tins  lower -bound  as  possible.  The  prognun  fhul  impleii.onts  this  random 
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optimization  U'chnique  is  called  iLogSim  (Cunent  Logic  Simulator), 

fn  our  expe..niiient,s,  iLogSim  ca,  :uiaics  the  lower  bound  by  trying  out  several  thousand 
I  ily  generated  input  patterns.  However,  such  a  technique  does  not  make  any  use  of  the 
current  waveform  of  a  currently  applied  input  pattern  in  generating  the  next  (hopefully  better) 
input  r>attc:‘n.  By  making  use  of  th;s  information,  we  can  generate  better  input  patterns  and 
thu.s  avoid  wasting  time  in  raridomly  enumerating  the  input  space.  We  have  experimented  with 
one  such  iterative  optimization  sche.nie,  namely  the  simulated  annealing  (SA)  algorithm  [41], 
[41:].  However,  simulated  annealing  algorithm  needs  an  objective  function  to  indicate  how 
well  the  search  is  progressing.  We  have  supplied  the  peak  of  the  current  waveforrr^  as  the 
objective  function  to  the  algorithm.  Therefore,  in  oui  v^AiieihuCiila!  results,  we  compare  the 
peak  of  the  current  v.'aveform  obtained  from  iMax  to  that  obtained  from  SA.  If  one  is 
interested  in  a  comparison  of  the  wavefonns  at  different  or  several  places  (and  not  just  the 
peaks),  the  lower  bound  waveform  should  be  optimized  at  several  places  by  the  SA  algorithm, 
or  the  iLogSim  program  should  be  run  for  a  large  number  of  input  patterns. 

4.3.6  Experimental  Results 

In  all  the  circuit  examples  considered  below,  two  assumptions  are  made.  First,  a  fixet 
number  is  assigned  to  each  gate  as  il.s  delay  value.  This  delay  value  can  be  different  for 
different  gates.  Second,  the  peak  of  the  transition  current  for  every  gate  for  boti!  Ih  and  hi 
transitions  is  taken  to  be  2  units  of  current.  The  results  of  simulated  annealing  were  collected 
after  tiding  about  100, (XX)  input  patterns  for  each  circuit.  Without  any  detailed  layout  infor¬ 
mation,  vve  have  assumed  that  all  the  gates  in  a  circuit  are  connected  to  a  single  contact  point, 
and  we  report  the  peak  values  of  the  current  wa'  eform^  obtained  from  iMax  and  SA. 
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Table  4.1  lii.,rs  the  results  of  running  iMax  and  SA  algorithms  on  eight  small  CMOS  cir¬ 
cuits.  The  number  (10)  next  to  iMax  in  the  table  indicates  the  value  of  Max_No_H!)ps 
parameter.  In  most  of  the  cases,  the  results  of  iMax  are  in  perfect  agreement  with  the  SA 
results. 


Table  4.1  iMax  and  SA  results  for  8  small  circuits 


Circuit 

No.  Gates 

No.  Inputs 

iMaxlO 
-  .  ,  .  _ 1 

SA 

Ratio 

BCD  Decoder 

18 

4 

32.41 

32.41 

1.00 

Comparator  A 

31 

11 

54.47 

54.47 

1.00 

Decoder 

16 

6 

26.75 

26.75 

1.00 

P.  Decoder  A 

29 

9 

44.11 

U.  1 1 

1.00 

P.  Decoder  B 

31 

9 

49.35 

49.35 

1.00 

Full  Adder 

36 

9 

58.00 

55.4-‘t 

1.05 

Parity 

46 

9 

47.57 

47.57 

1.00 

Aiu  (SN74181) 

63 

14 

78.89 

70.87 

1.11 

In  Table  4.2,  we  report  similar  results  on  the  ten  ISCAS-85  benchmark  circuits  [6].  We 
observe  that  for  all  the  ci.''cuits,  the  linear  time  iMax  algorithm  took  only  a  few  seconds  of 
cpu  time  on  a  sun  SPARCstation  ELC  compared  to  to  several  hours  of  time  needed  by  the  SA 
algorithm  (typical  times  needed  for  trying  10,000  patterns  by  SA  are  shown  in  the  table). 
Furthermore,  for  most  of  these  circuits,  the  ratio  of  the  upper  bound  obtained  from  iMax  to 
the  lower  bound  obtained  froir  SA  is  less  than  1.57.  There  are  two  possible  reasons  for  this 
mismatch.  Firstly,  it  is  quite  possible  that  the  lower  bound  obtained  from  SA  is  not  very 
close  to  the  MEC  waveform.  Since  all  the  circuits  have  at  least  .32  inputs,  the  space  of  possi¬ 
ble  input  patterns  is  huge,  and  the  lower  bound  v/aveforrn  obtained  after  trying  about  100,(X)0 
input  patterns  may  not  be  very  close  to  the  MEC  For  smaller  circuits  (Table  4.1)  where  the 
input  space  is  not  so  huge,  lOO.tXX)  input  pattern::  were  enough  to  get  a  lower  bound  estimate 
uf  MEiC  waveform  that  is  quite  close  to  the  actual  waveform,  as  is  ir  fleeted  by  the  re  suit, s. 
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The  second  possible  source  of  misn’.atch  is  our  conservative  independence  assumption  for  the 
signals  at  the  gate  inputs.  One  can  improve  on  this  assumption  by  attempting  to  resolve  the 


signal  correlations,  as  discussed  in  the  following  sections. 


Table  4.2  iMax  and  SA  results  for  10  ISCAS-85  circuits. 

No. 

No. 

Peak  Currents 

i  CPU  Times 

Circuit 

Gates 

Inputs 

iMax  10 

SA 

Ratio 

iMax  10 

SA(lOk) 

c432 

160 

36 

181.9 

162.0 

1.12 

1.2s 

9m  40s 

c499 

202 

41 

247.9 

186.6 

1.33 

1.2s 

10m  46s 

c880 

383 

60 

418.5 

321.4 

3.0s 

26m  32s 

C1355 

546 

41 

633.8 

415.8 

1.52 

4.6s 

36m  18s 

cl  908 

880 

33 

732.1 

445.6 

■EBI 

7.6s 

Ih  34m 

1193 

233 

1169.2 

866.1 

1.35 

7.9s 

2h  14m 

C3540 

1740.4 

866.1 

2.01 

c53l5 

2307 

178 

2312.6 

1558.9 

■E9I 

16.0s 

4h  40m 

C6288 

2406 

32 

3193.2 

1.28 

37.8s 

61h  58m 

c7552 

3512 

207 

4339.7 

2768.6 

1.57 

28.4s 

7h  28m 

We  next  discuss  the  effect  of  varying  the  Max  No_Hops  parameter  on  the  performance 
of  iMax.  Table  4,3  lists  the  peak  values  of  the  upper  bound  waveforms  for  lSCAS-85  circuits 
for  different  values  of  Max_No_Hops.  In  parentheses,  we  also  tabulate  the  cpu  times  (in 
sec.)  needed  by  the  algorithm.  As  the  value  of  Max  No  Hops  increase  ;,  the  number  of 
•ntervals  being  merged  at  every  gate  decreases  and  therefore,  the  cpu  time  needed  by  the  algo¬ 
rithm  increases.  This  also  improves  the  peak  value  of  the  current  waveform,  as  shown  in  the 
table.  However,  with  an  increa.se  in  Max_No_Hops  parameter,  while  the  cpu  time  continues 
to  increase,  the  improvement  in  peak  value  is  not  significant  beyond  Max  No  Hops  =  10. 
.Similar  behavior  is  observed  for  the  entire  current  waveform  and  not  just  for  the  peak  value. 
In  F-ig.  4.7,  we  plot  the  upper  bound  waveforms  for  cl 908  for  three  different  values  ot 
Max  No  Hopjc  The  difteicnce  between  the  upper  bound  waveforms  for  iMax  10  and  iMax<« 
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is  almost  negligible.  Sinnilar  plots  are  obtained  for  all  other  circuits.  Therefore,  we  conclude 
that  a  value  between  5  and  10  seems  to  be  a  good  choice  for  Max_No_Hops  parameter. 


Table  4.3  iMax  results  vs.  Max_No_Hops 

Circuit 

iMax:  Max  No  Hops 

1 

5 

10 

OO 

c432 

236.3  (0.4) 

185.1  (0.9) 

181.9  (1.2) 

181.7  (2.1) 

c499 

292.1  (0.4) 

247.9  (0.8) 

247.9  (1.2) 

247.9  (1.3) 

c880 

550.6  (0.8) 

424.1  (2.0) 

418.5  (3.0) 

415.1  (11.8) 

cl355 

749.9  (1.3) 

"646.8  (2.7) 

633.8  (4.6) 

633.8  (19.5) 

C1908 

908.3  (2.0) 

740.5  (5.0) 

732.1  (7.6) 

724.9  (86.3) 

C2670 

1476.0  (2.3) 

1188.6  (5.4) 

1169.2  (7.9) 

1166.4  (21.7) 

c3540 

^2088.2  (3.7) 

1752.4  (8.2) 

1740.4  (12.7) 

1730.3  (170.0) 

2632.0  (5,1) 

236.5  (11.2) 

2312.6  (16.0) 

2309.0  (109.5) 

c6288 

4134.8  (10.4) 

4098.4  (20.7) 

4096.2  (37.8) 

4096.1  (7086.0) 

c7552 

5163.1  (8.0) 

4401.6  (20.0) 

4339.7  (28.4) 

4325.8  (177.8) 

4.4  The  Signal  Correlation  Problem 

In  general,  signals  at  internal  nodes  of  a  circuit  are  correlated.  This  limits  the  number  of 
transitions  that  can  possibly  occur  at  the  outputs  of  the  gates,  an  effect  that  is  ignored  by  the 
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Figure  4,7  iMax  current  waveform:,  for  different  values  of  Max  No  Hops  parameter. 


iMax  algorithiT..  Two  examples  of  how  signal  correlation  limits  the  number  of  transitions  are 
illustrated  in  Fig.  4.8. 

In  Fig.  4.8(a),  signal  lines  xl  and  a:2  are  correlated  (in  tliis  case,  they  carry  the  same 
signal).  Depending  upon  the  specific  excitation  present  at  x,  only  one  of  the  two  gates  can 
switch  at  a  time.  However,  since  iMax  ignores  the  signal  correlation  present  between  x  I  and 
r2,  it  calculates  the  uncertainty  sets  at  the  outputs  of  the  two  gates  as  shown  in  the  figure  and 
thus  erroneously  concludes  that  Loth  gates  inay  switch  at  tlie  same  time,  and  therefore  adds 
two  triangular  current  pulses  due  to  both  gates  switching  simultaneously  to  the  overall  current 
waveform.  It  is  this  kind  of  approximation  that  contributes  to  a  loose  iMax  upper  bound. 
Similarly,  in  Fig.  4.8(b),  the  output  of  the  inverter  is  correlated  with  its  input  and  therefore, 
the  NAND  gate  may  never  switch.  However,  ignoring  this  correlation,  LMax  concludes  that 
the  NAND  gate  can  switch. 

\s  is  clear  from  these  examples,  the  source  of  the  signal  correlation  problem,  in  general, 
is  a  gate  (or  input)  whose  output  fans  out  to  several  other  gates.  Such  gates  are  called  multiple 


(a)  At  most  one  of  the  two 
gates  can  switch. 


(b)  Corrclat.^d  signals  .at  the  inputs  of  the 
NAND  gate  block  the  transition 


Figure  4  8  Two  examples  to  illustrate  the  signal  correlation  problem. 
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fan-out  (MFO)  gates.  The  general  situation  is  shown  in  Fig.  4.9,  where  a  MFO  gate  G  with 
output  node  n  fans  out  to  nodes  n^,  n2,  ■  ■  ■  ,  that  in  turn  feed  gates  G|,  Gt,  •  •  •  ,  G^. 
In  this  figure,  inputs  to  the  gates  Gj,  G2,  ■  ■  ,  Gi^  (which  are  nj,  ‘  respectively) 

are  correlated.  Due  to  this  correlation,  even  though  the  output  of  each  gate  (G,,  G2,  •  •  •  , 
Gjt)  can  assume  all  possible  excitations  as  calculated  by  iMax,  they  may  not  simultaneously 
carry  their  worst  case  excitations.  As  one  goes  deeper  into  the  circuit,  where  these  correlated 
outputs  reconverge  and  feed  the  same  gate,  the  inputs  of  that  gate  become  correlated  (e.g., 
NAND  gate  in  Fig.  4.8(b)).  Such  gates  are  called  reconvergent  fan-out  (RFO)  gates.  With 
correlated  signals  at  the  inputs  of  a  gate,  the  number  of  transitions  that  can  possibly  occur  at 
its  output  is  reduced.  The  signal  correlations  considered  above,  which  exist  among  various 
nodes  throughout  the  circuit,  aie  called  spatial  correlations. 

Besides  the  spatial  correlations,  there  is  another  set  of  correlations  that  the  iMax  algo¬ 
rithm  completely  ignores.  The  excitation  assumed  by  a  node  at  time  r,  restricts  the  set  of  pos- 


Figure  4,9  An  example  of  a  Multiple  Fan  out  Gate. 
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sible  excitations  that  the  node  can  assume  at  an  earlier  or  a  later  time.  For  example,  if  a  node 
is  low  at  time  t ,  then  it  can  either  stay  at  low  or  switch  from  high  to  low  at  time  t~  and  it  can 
either  stay  at  low  or  switch  from  low  to  high  at  time  r"*".  These  correlations  which  exist  in  the 
time  domain  are  called  temporal  correlations. 

The  iMax  algorithm  completely  ignores  all  spatial  and  temporal  signal  correlations  and, 
therefore,  overestimates  the  supply  currents.  The  advantage  of  ignoring  correlations  in  the 
algorithm  is  its,  very  desirable,  linear  time  performance. 

4.5  Resolving  Signal  Correlations 

The  upper  bound  produced  by  the  iMax  algorithm  can  be  made  exact  by  doing  a  brute- 
force  enumeration  at  the  inputs  of  the  circuit  and  calculating  an  envelope  of  the  current 
waveforms  produced.  In  enumeration,  since  unambiguous  input  pattern.s  are  applied  to  the 
circuit,  there  is  no  "uncertainty"  present  at  the  inputs  and  therefore,  signal  correlations  do  not 
become  an  issue.  In  a  similar  fa.shio"  one  can  improve  the  results  of  the  iMax  algorithm  by 
doing  a  partial  enumeration  at  a  few  selected  nodes  in  the  circuit. 

An  example  of  how  partial  enumeration  helps  improve  the  upper  bound  can  be  seen  from 
Fig.  4.8(a).  In  this  circuit  with  no  enumeration,  iMax  would  assume  that  the  signal  lines  .t  1 
and  x7.  are  mutually  independent  and  therefoie  infer  that  b(>th  NAND  imd  NOR  gates  can 
switen  at  the  same  time.  However,  if  we  do  ptutial  enumeration  at  signal  line  .r ,  then  we 
would  generate  four  cases  coiTesponding  to  when  .x  =  /  ,  x  =  h. ,  x  -  hi  and  v  —  Ih.  When 
=  /  or  /i  / ,  only  the  NOR  gate  switches.  Similarly,  when  x  -  /,  or  !h  ,  only  the  NAND  gate 
switches.  Thus,  by  splitting  the  problem  into  four  sub-problems.,  vve  have  improved  our  result, 
i.e.,  found  that  only  one  of  the  two  gates  may  switch  at  any  given  une. 
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While  enumerating  a  node,  we  only  need  to  process  a  small  subset  of  the  gates.  We 
define  the  COne  of  INfluencc,  COIN{n),  of  a  node  n  as  the  set  of  all  the  gates  that  can  possi¬ 
bly  be  affected  by  a  change  in  excitation  at  the  node.  Thus,  a  gate  is  in  the  COlN(n)  of  a 
node  n  if  it  is  either  directly  fed  by  n  or  is  connected  to  the  output  of  a  gate  that  is  in 
COIN(n).  While  enumerating  a  node,  we  only  need  to  consider  those  gates  that  are  in  its 
COIN. 

One  technique  to  partially  enumerate  the  internal  nodes  of  a  circuit,  called  Multi-Cone 
Analysis  (MCA),  was  reported  in  [43],  The  motivation  behind  such  an  approach  was  to  be 
able  to  enumerate  at  the  outputs  of  the  MFO  gates,  which  are  the  sources  of  the  signal  corre¬ 
lation  problem.  However,  from  the  results  in  [43]  and  Tables  4.6,  4.7  below,  it  can  be  seen 
that  the  MCA  approach  offers  only  a  modest  improvement  in  the  upper  bound.  There  are 
.several  reasons  for  this. 

As  shown  in  Table  4.4,  there  are  usually  several  MFO  gales/inputs  in  a  circuit  and  all  of 
these  nodes  should  be  enumerated  to  properly  re.solve  the  signal  correlation  problem.  From 
our  experience  with  lSCAS-85  benchmark  circuits,  we  have  found  that  the  COINs  of  several 
of  these  nodes  overlap  and  therefore,  to  properly  handle  signal  correlations,  these  nodes 
should  be  enumerated  simuitaneously.  Furthermore,  because  of  the  presence  of  glitches  in  a 
circuit,  signals  at  internal  nodes  span  several  time  points  (i.e.,  signal  transitions  occur  at 
secerai  time  points).  To  take  care  of  the  temporal  correlation  problem,  itic  u  jJc  .  hohld  lie 
enumerated  at  each  of  these  time  points  Simultaneous  enumeration  is  an  extremely  expensive 
process  specially  vvhen  there  are  several  nodes  and  each  node  needs  to  he  enumerated  tit 
several  time  points.  For  example,  to  enumerate  two  nodes  simultaneously,  the  epu  tune 
needed  to  enumerate  each  node  gets  multiplied  To  avoid  this  multipliea'ive  growth  of  epu 


time,  we  made  several  simplifying  assumptions  in  the  implementation  of  the  MCA  algorithm 
[43].  Because  of  these  simplifications,  the  algorithm  led  to  only  mild  improvement  in  iMax 
results. 


Table  4.4  Number  of  MFO  gates/inputs  in  ISCAS-85  circuits. 

Circuit 

No.  Inputs 

No.  MFO 

Circuit 

No.  Inputs 

No.  MFO 

c432 

36 

153 

c2670 

233 

1129 

c499 

41 

170 

c3540 

50 

1647 

c880 

60 

357 

c5315 

178 

2184 

cl355 

41 

514 

c6288 

32 

2384 

cl908 

33 

855 

c7552 

207 

3405 

There  are  usually  several  RFO  gates  in  a  circuit.  If  at  the  output  of  a  RFO  gate,  a  false 
transition  is  predicted  by  iMax  (.such  as  in  Fig.  4.8(b)),  then  as  this  transition  propagates 
through  the  circuit  to  the  output  nodes,  it  causes  .several  other  false  transitions  in  the  circuit 
along  its  way.  Therefore,  locating  and  suppressing  such  transitions  in  iMax  is  important. 
However,  these  fal.se  transitions  can  be  isolated  only  by  doing  a  simultaneous  enumeration  at 
the  source  MFO  gate(s)  that  created  (hem.  In  effect,  one  needs  to  construct  tl;e  supers^ute  [16] 
for  each  RFO  node  in  the  circu't  ’"d  for  each  supergate,  do  a  simultaneous  enumeration  at  its 
MFO  inputs.  However,  these  supergates  can  be  as  big  as  the  entire  circuit  and  therefore 
enumerating  their  inputs  becomes  intractable.  We  have  implemented  an  approximate  algorithm 
based  on  enumerating  primary  stem  regions  [45 J,  I46j.  Our  results  show  rnly  moilesi 
improvement  in  the  upper  hound  waveforms. 

From  our  experience  with  enumerating  internal  nodes  of  a  circuit,  as  explained  above, 
we  infer  that  improving  iMax  upper  bound  wavefoirns  by  enumerating  internal  nodes  is  very 
expensive  and  does  not  offer  a  practieal  solution  for  VL.SI  circuits.  In  the  next  section,  we 
present  an  alternative  partial  input  enumeration  approach  that  significantly  improves  the  iMax 
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results  and  represents  a  good  speed- accuracy  trade-off. 


4.6  Partial  Input  Enumeration  (PIE) 

As  shown  in  Table  4.4,  there  are  usually  many  more  MFO  nodes  than  primary  inputs  in 
a  circuit.  Secondly,  as  stated  in  Section  4.2,  all  the  inputs  to  a  circuit  switch  at  most  once  at 
time  zero.  Therefore,  there  is  only  one  time  point  at  which  a  primary  input  needs  to  be 
enumerated.  This  is  in  contrast  to  an  internal  circuit  node  which  usually  needs  to  be 
enumerated  at  several  time  points.  These  observations,  combined  with  the  fact  that  iMax  is 
an  extremely  fast  algorithm  led  us  to  explore  the  following  partial  input  enumeration  (PIE) 
algorithm  to  improve  the  iMax  upper  bound  at  every  contact  point. 

4.6.1  The  Algorithm 

Let  -V],  .t2,  ■  •  ■  ,  Xf^  be  the  N  primary  inputs  of  a  circuit  under  consideration.  Let  X, 
represent  the  uncertainty  set  for  input  x-  at  time  zero.  The  input  search  space  for  the  circuit 
consists  of  the  set  of  all  valid  input  patterns  that  can  be  applied  to  the  circuit.  Mathemati¬ 
cally,  tne  input  search  space  is  j(e,,t*2,...,<'/v)  I  e  |  g  X],e2  g  X2 . Cyy  g  X/^).  For  brevity, 

we  denote  this  by  (XiX2,...,X;v)-  Suppose,  for  the  purposes  of  this  illustration,  for  a  particu¬ 
lar  input  ,v, ,  X,  =  X.  Then  the  input  search  space  (X,,.Y2....X7v)  ihe  circuit  can  be  divided 
into  four  disjoint  parts,  namely  (.V  ,.,V2....,{  /  ),.-.,X,,; ),  (X  i.X  i,.  .,!// i,..,,X^. ), 

(.V|,X2 . {h  I  }  and  <X  j,.Y  2, ...,{/// We  can  compute  the  nia.ximuni  current 

waveforms  (at  ever;  contact  point)  for  each  of  thc.se  four  parts  by  running  the  iMax  algorithm 
and  in  eacfi  case,  restricting  the  excitation  on  input  .r,  to  the  value  in  its  respective  iincertaintv 
subset.  Since  the  four  parts  combined  together  con,stitute  the  complete  search  space,  by  taking 
an  upper  bound  envelo[X'  of  the  t'oiii  current  v/aveforms  at  every  contact  point,  we  can  still 
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guarantee  an  upper  bound  on  the  respective  MEC  waveforms.  Since,  in  each  of  the  four  runs 
of  iMax,  specific  excitations  are  present  at  input  jc,  ,  signal  correlations  due  to  x,  disappear 
and  the  resulting  current  waveform  should  be  an  improvemen  on  the  original  upper  bound. 
In  a  similar  fashion,  the  upper  bounds  for  the  individual  subcases  can  be  improved. 

The  set  of  inputs  selected  for  enumeration  has  a  direct  influence  on  the  quality  as  well  as 
the  cost  of  the  solution  obtained.  If  all  the  inputs  are  selected  and  enumerated  then  the  upper 
bound  obtained  at  every  contact  point  would  be  exact.  However,  doing  this  is  practically 
impossible  for  most  circuits.  The  extent  to  which  an  input  contributes  to  signal  correlations 
inside  a  circuit  is  different  for  different  inputs.  For  example,  in  Fig.  4.8(a),  enumerating  input 
x  is  more  beneficial  than  enumerating  any  of  the  other  two  inputs.  Hence,  by  selecting  and 
enumerating  inputs  in  an  intelligent  fashion,  we  can  significantly  improve  the  iMax  upper 
bounds,  without  spending  too  much  cpu  time. 

We  have  developed  an  intelligent  be.\t  first  search  (BFS)  algorithm  [47]  that  is  very 
effective  in  selecting  and  enumerating  inputs  and  thereby  improving  the  upper  bound  at  every 
contact  point.  The  algorithm  starts  with  the  iMax  upper  bound  and  some  knowr.  lower  bound 
and  repeatedly  expands  "search  nodes"  (.s_nodes)  in  a  best  first  fashiv,n.  Various  s  nodes, 
generated  during  the  search,  correspond  to  partially  specified  input  states  in  which  some 
inputs  have  specified  excitations  (e.g.,  say  low),  while  others  have  "uncertain"  (or  unspecified) 
excitations.  An  objective  function  is  a.s.soeiated  with  the  search  and  during  the  seaich, 
s  niKles  w  hich  correspond  to  the  best  objective  value  are  repeatcdlv  expanded,  rtiis  function 
will  be  explained  later  in  this  section.  Because  of  this  best  first  strategy,  there  is  a  gradual 
reduction  in  the  value  (at  any  time)  of  the  upper  bound  w'aveform  at  eveiv  contact  point  in  the 
circuit.  This  iterative  iniprovcrnent  is  a  very  important  feature  of  the  algorithm  lor  large 
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circuits  where  an  exhaustive  exploration  of  the  input  space  is  practically  impossible.  The  BPS 
algorithm  can  be  stopped  at  any  intermediate  stage  and  the  current  best  upper  bounds  at  vari¬ 
ous  contact  points  can  still  be  reported. 

The  BPS  search  starts  with  the  initial  uncertain  state  i.e.,  s_node  =  (X^,  X 2,  •  •  ■  ,  Xy) 
and  a  known  lower  bound,  which  is  the  objective  value  for  some  input  pattern.  During  the 
.search,  a  s  node  with  the  highest  objective  value  is  repeatedly  selected  and  its  descendent 
s  nodes  are  generated  by  enumerating  an  input,  as  explained  in  the  following  outline: 

Remark:  List  is  an  ordered  list  of  s_nodes,  arranged  in  decreasing  objective  values. 

1  List  <r-  starting  s_node  (Initial  uncertain  state). 

Upper  Bound  <—  objective  value  of  the  starting  s  node. 

Lower  Bound  <—  objective  value  for  a  specific  input  pattern,  (otherwise  0.0). 

2.  While  Stopping  Criterion  is  not  satisfied,  do 

2. 1  Remove  the  top  s_node  from  the  List. 

2.2  Calculate  next  input  number  to  enumerate  from  the  Splitting  Criterion. 

2.3  Generate  all  (<  4)  children  s  nodes  by  enumerating  the  input  and  calculate  their  objec¬ 
tive  values. 

2.4  If  these  children  are  leaf  s  nodes,  then  update  the  Lower  Bound,  else,  insert  them  in 
List,  after  paining  if  any. 

2.5  Upper  Bound  objective  value  of  the  top  s  node  in  List. 

3.  Report  the  best  upper  bound  current  waveform  at  every  contact  point.  STOP. 

where  a  leaf  .s  node  is  one  which  has  exactly  one  excitation  associated  with  eadi  of  its  uncer¬ 
tainty  sets  (i.e.,  it  corresponds  to  an  input  pattern),  The  tullowirig  functions  are  u.sed  in  the 
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outline  above. 


Objective  Function:  This  function  .specifies  the  quantity  that  is  being  minimized  during 
the  search.  There  are  usually  several  contact  points  within  a  combinational  block  and  we 
would  like  to  minimize  the  (iMax)  upper  bound  waveform  at  each  of  She.se  points,  so  that 
they  arc  close  to  their  respective  MIEC  waveforms.  One  possibility  is  to  minimize  the  peak  of 
a  weighted  sum  of  the  upper  bound  waveforms,  where  these  weights  are  determined  depend¬ 
ing  upon  how  much  "influence "  the  contact  point  has  on  the  overall  voltage  drops.  We  .m* 
currently  working  on  this  problem.  In  the  experiments  reported  below,  we  have  assumed  these 
weights  to  be  unity  and  wc  minimize  the  peak  value  of  the  sum  of  all  the  upper  bound 
waveforms.  This  corresponds  to  minimizing  the  worst  ca.se  total  current  of  the  combinational 
block. 

Sioppiug  Criterion:  We  stop  the  .search  if  any  of  the  following  two  conditions  is 
satisfied. 

a.  Best  Upper  Bound  <  Lower  Bound  x  ETF, 

b.  Number  of  s  nodes  generated  >  User  specified  parameter  (Max  No  Nodes). 

The  Error  Tolerance  Factor  (ETF)  is  a  user-specified  parameter  that  provides  control 
over  the  final  desired  accuracy  of  the  algorithm.  The  value  of  this  parameter  is  always  bigger 
than  ) .  The  first  condition  above  sjjecifies  that  when  the  highest  upper  bound  value  is  within 
ETF  factor  of  some  known  lower  bound,  then  the  search  can  be  tenninated.  In  large  circuits 
where  calculating  an  exact  upper  bound  value  by  running  the  search  to  completion  is 
extremely  expensive,  and  an  overestimation  by  20%  to  may  be  acceptable,  such  a  param¬ 
eter  can  be  extremely  useful.  The  second  condition  puis  a  hard  limit  (Max  No  Nodes)  on 


the  number  of  sjtiodes  that  arc  to  he  generated  in  the  search. 

Pruning  Criterion:  During  the  .search,  if  we  come  across  a  s  node  for  which  the  upper 
bound  satisfies  the  following  condition. 

Upper  Bound  S  Lower  Bound  x  ETF 

then,  such  a  s  node  can  be  deleted  from  the  search  as  its  upper  bound  value  is  already 
acceptable  This  pruning  critenon  deletes  several  ''Unnecessary  s  nodea  during  the  search  and 
thus  keeps  the  n?err»ory  usage  down. 

Splitting  Criterion  (SC>:  TThis  cnicnon  specifics  the  input  which  should  be  enumerated 
next  from  any  s  node  during  tlie  seasch.  In  the  next  subsection,  we  describe  two  heuristic 
functions  for  doing  this. 

TTic  BPS  algorithm  always  processes  s  nodes  which  arc  on  its  current  wavefront,  see 
Fig.  4.10.  A:  the  start,  this  wavefront  consists  of  only  one  s  node,  namely  the  initial  uncer¬ 
tain  state  As  the  search  progresses,  this  wavefront  moves  forward  through  the  input  search 
space.  An  input  pattern  leading  to  the  maximum  objective  value  could  belong  to  any  of  the 
s  nodes  on  the  wavefront.  Tlicrcfore,  when  the  algorithm  tcmiinaies,  at  every  contact  point, 
we  compute  an  envelope  of  the  current  waveforms  of  all  the  s  nodes  on  the  wavefront  and 
report  that  as  an  upper  bound  wavefonn. 

4.6.2  Heuristic  $C  and  Experimental  Results 

We  now  describe  two  heuristics  for  the  splitting  criterion  that  have  shown  good  results  in 
praclic'.;.  The  first  heuristic  .selects  an  input  which  has  the  highest  sensitivity  while  the  second 
one  seietts  an  input  based  upon  the  influence  it  has  inside  the  ciicuit, 


Figure  4.  iO  Exploring  s  nodes  through  the  input  search  space. 

4.6.2. 1  //|  Heurisfic 

Let  us  suppose  that  during  the  search,  we  are  at  a  particular  s_node  n  and  we  select  an 
input  for  enumeration.  If  we  assume  that  the  uncertainty  set  for  x,  at  time  zero  is  X,  then 
by  enumerating  x,,  we  would  generate  four  children  s_nodes,  as  shown  in  Fig.  4.11.  We 
assume  that  the  objective  value  of  s_node  n  is  denoted  by  obj„  and  the  objective  values  of 
the  children  s  nodes  arc  denoted  by  obJi,  obj^,  obj^i  and  objtt,.  If 

A  obJi  -  obj„  -  max^bj,,objf,,obj^t,obj  1^, 

then  by  enumerating  x, ,  we  can  improve  the  objective  value  of  s_node  n  by  an  amount  Aobj  . 
We  can  repeat  this  calculation  for  every  input  and  then  select  an  input  which  gives  rise  to  the 
maximum  improvement  in  the  objective  value. 

However,  if  AobJi  is  zero  for  all  the  inputs,  which  happens  very  often  in  practice,  then 
the  above  selection  process  would  not  work  well.  For  a  specific  input  x, ,  Aobj^  =  0  means 
that  the  objective  value  of  one  of  its  children  s_nodes  is  equal  to  obj„.  However,  for  the 
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Figure  4. 1 1  An  example  to  illustrate  the  // , . 

remaining  children  s_nodes.  the  objective  value  may  not  be  equal  to  obj„  and  this  informa¬ 
tion  can  be  used  in  assigning  credit  (or  relative  importance)  to  x,  .  Based  on  these  obser\'a- 
tions,  we  have  come  up  with  the  following  H_  1  heuristic  function  corresponding  to  an  input 

x,; 

H i  ~  A  X  {obj„  -obj  |)  +  fi  X  (obj„  -obj{)  +  C  x  {abj„  -obj-^)  +  i.obj„  -ohj^) 
where  obj obj 2,  objj  and  obj^  are  the  objective  values  of  the  children  s  nodes,  generated 
by  enumerating  x^ ,  arranged  in  decreasing  order  and  A ,  B  and  C  are  three  constants  such  that 
A  :>  B  »C  »1.  At  any  s  node  during  the  search,  we  compute  the  heuristic  value*  for  every 
input  and  select  an  input  with  the  maximum  associated  heuristic  value.  This  splitting  criterion 
is  called  dynamic  (H splitting  criterion  because  at  each  s_nodc,  it  calculates  the  heuristic 
value  for  every  input  and  then  selects  the  best  input  for  enumeration. 

The  results  of  partial  input  enumeration  using  the  BFS  algorithm  and  using  the  dynamic 
(//|)  splitting  criterion  for  the  nine  small  circuits  are  documented  in  Table  4,5,  For  ail  the 
circuits,  the  algorithm  was  run  to  completion,  i.e.,  till  the  upper  bound  became  equal  to  the 
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lower  bound  (ETF  =  1).  The  results  clearly  show  that  the  PIE  algorithm  is  very  efficient  in 
scanning  the  input  space.  As  an  example,  the  last  circuit  in  the  table  (Alu)  has  14  inputs  and 
therefore,  the  number  of  possible  input  patterns  for  this  circuit  is  4*^  =  268,435,456.  The  PIE 
algorithm  was  able  to  scan  the  entire  search  space  after  generating  just  233  s_nodes.  This  also 
shows  that  the  upper  bound  produced  by  the  iMax  algorithm  i  very  tight  for  these  circuits. 


Table  4.5  Results  of  PIE  for  9  small  circuits. 

Circuit 

Dynamic  {H ^)  Splitting  Criterion 

Static  (H 1 )  Splitting  Criterion 

No.  S  nodes 

Generated 

iMax  mns 

in  SC 

Time 

No.  S  nodes 

Generated 

iMax  runs 
in  SC 

Time 

BCD  Decoder 

17 

40 

\Sts 

17 

17 

1.2s 

Comparator  A 

109 

664 

87.7s 

277 

45 

35.1s 

Comparator  B 

45 

264 

30  6s 

45 

45 

9.7s 

Decoder 

25 

84 

3.6s 

25 

25 

1.8s 

P.  Decoder  A 

37 

O 

00 

11.4s 

37 

37 

.3.9s 

P,  Decoder  B 

37 

180 

11.7s 

37 

37 

4.2s 

Full  Adder 

53 

296 

29.2s 

53 

37 

7.7.S 

Parity 

37 

180 

18.7s 

37 

37 

6.4s 

Alu  (SN74I81) 

233 

1952 

_ 1 

378.2s 

252.5 

57 

.352,7s 

As  can  be  seen  from  Table  4.5,  the  number  of  iMax  runs  needed  in  the  dynamic  splitting 
criterion  far  exceeds  the  number  of  s  nodes  generated.  At  a  s  nodc,  to  calculate  the 
value  for  a  particular  input  x, ,  we  need  to  run  the  iMax  algorithm  IX,  !  times.  If  a  s  node 
has  k  inputs  which  are  possible  candidates  for  enumeration  (i.e.,  thei.»  IX,  i  ;>  1),  then  we  need 
to  run  the  iMax  algorithm  J^il-i'X,  I  number  of  times  to  find  the  best  input  to  enumerate  next. 
For  bigger  circuits,  with  large  number  of  inputs,  this  time  will  be  even  more  dominant  render¬ 
ing  the  PIE  algorithm  prohibitively  expen*;ivt;.  Therefore,  we  have  experimented  with  other 


less  expensive  alternatives. 


Instead  of  calculating  the  heuristic  function  value  (f/  [)  for  every  input  at  every  s  node 
during  the  search,  we  calculate  the  heuristic  value  for  every  input  at  the  beginning  of  the 
search.  All  the  inputs  are  arranged  in  the  decreasing  order  of  their  heuristic  values.  During  the 
BFS  search,  inputs  are  selected  in  this  fixed  order.  This  criterion  is  called  static  (Hi)  splitting 
criterion.  The  amount  of  time  spent  in  the  static  splitting  criterion  is  fixed  and  is  equal  to 
2^,^,  lx,  I  runs  of  the  iMax  algorithm  for  a  circuit  with  N  inputs.  The  results  of  the  PIE  algo¬ 
rithm  using  the  static  (//,)  splitting  criterion  are  also  summarized  in  Table  4.5.  With  static 
splitting  criterion,  the  number  of  runs  of  the  iMax  algorithm  needed  in  the  splitting  criterion 
goes  down,  but  the  number  of  s  nodes  generated  uunng  ihe  scan  h  g-  cs  up  for  some  circuits. 
However,  for  all  the  circuits,  we  observe  an  overall  reduction  in  the  epu  times  needed  for  the 
algorithm  to  complete. 

4.6.2.2  H  2  Heuristic 

The  number  of  gates  that  are  affected  by  a  change  in  excitation  at  an  input  is  a  gc-od 
heuristic  measure  of  how  much  influence  the  input  has  on  the  upper  bound  waveforms.  There¬ 
fore,  inputs  which  affect  more  number  of  gates  (i.c.,  which  have  larger  COINs)  should  be 
enumerated  before  otnes.  This  'cads  us  to  another  (static)  splitting  criterion  H 2,  in  which  we 
calculate  the  size  of  the  COIN  1  assc^^iated  with  each  input  x, .  As  with  Hi,  all  the  inputs  arc 
arranged  in  the  decreasing  order  of  H  2  values  and  during  the  search,  inputs  arc  selected  in 
this  fixed  order.  We  will  show  in  the  next  section  that,  while  both  static  //j  and  H2  give  good 
results  in  practice,  f/2  is  much  belter  in  terms  of  speed  and  has  accuracy  comparable  to  //j. 

The  results  of  partial  input  enumeration  using  both  //,  and  sialic  splitting  criteria  for 
the  ISC  AS- 85  benchmark  circuits  are  shown  in  Table  4  6.  In  the  table,  under  various  IMax, 
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MCA  and  BFS  columns,  we  show  the  ratio  of  the  respective  upper  bound  to  the  lower  bound 
obtained  from  simulated  annealing.  The  numbers  in  parentheses  under  the  BFS  columns  indi¬ 
cate  the  numt/e;  of  s_nodes  that  were  generated  before  stopping  the  search  (i.e.,  the 
Max  No  Nodes  parameter;  Ik  stands  for  1(XX)).  Total  cpu  times  needed  by  the  algorithm  on 
a  sun  SPAFfC  station  ELC  (with  Max_No_Nodes=;l(X))  are  also  shown  in  the  table.  From 
Table  4.6,  v  e  note  that  for  all  the  circuits,  the  ratio  of  the  upper  bound  to  the  lower  bound  is 
at  most  l.'j2  (as  opposed  to  a  worst  ca«e  of  2.02  for  the  simple  iMax  algorithm).  This  ratio 
can  be  further  improved  by  mnning  the  PIE  algorithm  for  lo''ger  durations.  We  emphasize 
that,  since  we  can  only  comp.are  the  upper  bound  to  a  lower  hound,  the  numbers  in  the  table 
are  only  upper  bounds  on  the  error.  I;  is  prohibitively  expensive  to  measure  the  true  error. 


Table  4  6  Results  of  PIE  for  10  lSCAS-85  circuits. 


■ 

MCA 

Static  //  jSC 

Static  H  iSC 

BFS 

(100) 

BFS 

(Ik) 

Time 

(100) 

BFS 

(100; 

BFS 

(Ik) 

rime 

MOO) 

c432 

1.12 

112 

1.08 

“ro5 

5m  1 4s 

“Til  ~ 

"TIT 

1  m  34  s 

c499 

1.33 

1.20 

1  33 

1.33 

4m  40s 

1.33 

1,33 

Im  23s 

c880 

1.31 

1.26 

1.25 

1.22 

17m  16s 

1.28 

1,26 

4m  5s 

cl  355 

1,52 

1.52 

1.52 

1  52 

21m  28s 

1.52 

^ITT' 

OiTi  1  3s 

cl908 

1.64 

1.55 

1.49 

1.46 

33m  i7.s 

1  58 

1.54 

1  i  m  51s 

c2670 

1.35 

1.34 

L29 

1,28 

Ih  57m 

1.35 

1.35 

llm  56s 

c354() 

2.01 

1.95 

1.45 

1.59 

1.37 

17m  3s 

c5315 

1,48 

1.44 

1.42 

1.48 

1.47 

26m  2s 

c6288 

1.28 

1.28 

1.28 

. 

1.27 

2h  5m 

1.28 

1  28  1  57m  28s 

c7552 

!  ,57  1  1 .55 

1.52 

1.50 

6h  21m 

1,53 

1.53 

45nri  4s  j 

While  the  improvement  over  the  original  IMax  algorithm  is  not  large  in  all  the  cases,  in 
those  cases  where  the  iMax  bound  war  veiy  loose,  such  as  c3540,  the  new  PIE  algorithm 
wiin  Hi  or  H2  heuristic  gives  significant  improvement:  the  ratio  of  2.02  (maximum  over- 


estimation  by  1.02)  is  now  1.37  (maximum  over-estimation  by  0.37)  with  H2,  a  reduction  in 
the  maximum  over-esd>nation  by  about  64%. 

We  also  emphasize  the  following  attractive  property  of  the  algorithm;  a  significant 
amount  of  improvement  in  the  upper  bound  occurs  in  the  first  few  s_nodes  (about  50-200)  of 
the  algorithm.  This  is  shown  in  Fig.  4.12  for  c3540,  where  the  ratio  of  the  upper  Ijound  to 
the  lower  bound  is  plotted  as  a  function  of  epu  time  fox  the  first  1000  s  nodes. 

The  epu  time  needed  for  generating  the  input  list  by  the  H2  splitting  criterion  i.s  negligi¬ 
ble  cempared  to  the  time  needed  by  the  criterion.  For  VLSI  circuits  with  several  hundred 
inputs,  where  the  time  needed  by  the  /f,  criterion  may  be  large,  H2  criterion  may  be  used 
instead.  As  can  be  seen  from  Table  4.6  (also  see  Table  4.7),  the  results  produced  by  using 
either  splitting  criteria  are  quite  comparable,  .specially  for  those  circuits  where  iMax  did  not 
produce  a  good  upper  bound. 


Figure  4.12  ‘Upper  Bound  /  Lower  Bound  vs  Time'  plot  for  c3540. 
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Table  4.7  Results  of  PIE  for  10  ISCAS-89  (Comb.)  circuits. 


Circuit 

m 

Static  H  ;SC 

Static  H  2SC 

BPS 

(100) 

BPS 

(Ik) 

MW 

BBl 

w 

si  423 

657 

1.35 

1.32 

1.32 

1.29 

37m  22s 

1.35 

1.34 

7m  43s 

sl488 

653 

2.21 

2.10 

1.40 

1.08 

5m  32s 

1.41 

1.06 

si  494 

647 

2.18 

2.08 

1.37 

1.06 

5m  35s 

1.39 

1.05 

2m  51s 

s5378 

2779 

1.38 

1.37 

1.29 

1.25 

2h  23  m 

1.30 

1,23 

13m  21s 

s9234 

5597 

1.76 

1.74 

1.51 

1.47 

7h  24m 

1.56 

_  .  j 

1.56 

37m  18s 

si  3207 

7951 

1.37 

1.35 

- 

- 

- 

36m  53s 

S15850 

.  9772 

1.81 

1.80 

- 

- 

- 

1.64 

1.57 

s35932 

16065 

1.66 

1.66 

- 

- 

- 

1.56 

1.56 

2h  6m 

S38417 

22179 

1.73 

1.70 

- 

- 

- 

1.72 

1.68 

2h  46m 

S38.584 

19253 

1.45 

1.38 

- 

- 

- 

1.39 

1.37 

2h  15m 

In  order  to  demonstrate  the  applicability  of  the  partial  input  enumeration  algorithm  for 
VLSI  circuits  with  several  thousand  gates,  we  have  also  experimented  with  the  ISCAS-89 
benchmark  circuits  [48].  For  these  synchronous  sequential  circuits,  we  have  extracted  the 
combinational  blocks  by  deleting  the  flip-flops.  The.se  combinational  blocks  have  gate  counts 
ranging  up  to  22,000  and  number  of  inputs  ranging  up  to  1750.  The  results  of  the  PIE  algo¬ 
rithm  on  some  of  the  lSr,A3-89  circuits  (combinational  blocks)  using  both  /¥,  and  //i  split¬ 
ting  criteria  are  sunimari?ed  in  Table  4.7.  It  is  clear  from  the  table  that  eeen  for  circuits  of 
t.'iis  size,  our  algorithms  show  good  speed  and  accuracy  performance. 
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5.  BUS  MODEL  EXTRACTION 

5.1  Introduction 

In  this  chapter,  we  describe  the  work  accomplished  on  the  extraction  bus  models  from 
layout  information.  The  aim  is  to  develop  fast,  accurate,  and  general  techniques  for  extracting 
the  RC  network  models  of  power  busses  in  VLSI  chip  design.  The  models  are  used  to  esti¬ 
mate  electromigration  and  voltage  drop  in  the  busses.  The  main  results  accomplished  on  the 
extraction  subtask  this  past  year,  compared  to  our  previous  work  are:  (1)  the  ability  to  extract 
the  resistive  models  of  bus  partitions  that  do  not  fit  precharacterized  ‘standard’  partitions;  this 
is  done  using  the  Boundary  Element  Method  (BEM);  (2)  improved  modeling  of  the  capaci¬ 
tance  by  the  use  of  finite  element  partitioning  and  model  reduction  techniques  to  reduce  the 
complexity  of  the  extracted  model  without  sacrificing  accuracy. 

5.2  Extraction  of  Resistance  Models 

For  simplicity  of  pre.sentation  and  without  loss  of  generality,  we  wiM  a.ssunie  i  sheet 
resistance  of  1  Ll/square.  The  layout  is  first  partitioned  into  primitive  geometrical  shapes. 
The  idea  behind  partit.’cning  is  to  simplify'  the  computation  since,  if  we  ignore  capacit;uice, 
the  sum  of  resistances  of  the  partitioned  pieces  will  equal  the  resistance  of  the  unpartitioned 
piece.  The  partitioning  strategy  relies  on  partitioning  a  section  of  the  layout  along  equipoten- 
tial  lines  that  will  occur  when  a  voltage  source  is  applied  to  one  of  the  terminals  and  the  oth¬ 
ers  are  grounded.  Usually  cquipotential  lines  that  arc  approximately  straight  are  chosen.  The 
distributions  of  cquipotential  lines  for  some  common  shaf)cs  arc  shown  in  Fig.  5.1.  For  two 
terminal  segments,  like  those  in  Fig.  5.1,  the  cquipotential  lines  will  be  the  same  no  matter 
which  terminal  is  the  grounded  one.  This  is  not  true  for  the  four-terminal  segment  in  Fig, 


102 


5.2,  thus  making  three  or  more  terminal  segments,  such  as  T  and  +  segments,  harder  to  parti¬ 
tion  along  equipotential  lines. 

After  partitioning,  a  resultant  N  terminal  shape  is  modeled  as  an  N  port  resistor  network 
having  no  internal  nodes,  where  each  resistor  connects  one  terminal  to  another,  yielding  N(N- 
l)/2  resistors.  Conceptually,  to  compute  (the  resistance  connecting  terminal  i  to  each 
other  terminal  j),  a  unit  voltage  is  applied  to  terminal  i  and  all  other  terminals  are  grounded. 
Let  Ij  be  the  current  flowing  through  terminal  j.  Then  =  l/Ij. 

If  partitioning  does  not  follow  the  equipotential  partition  constraint,  then  error  is  intro¬ 
duced.  By  applying  tlie  general  multiterminal  model  described  above  to  each  of  the  incorrectly 
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Figure  5.1  Potential  distribution  in  some  common  two-terminal  shapes. 
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Figure  5.2  Four-terminal  resistor. 

partitioned  pieces,  an  equipotential  line  is  forced  to  exist,  by  virtue  of  the  unity  voltage  excita¬ 
tion,  where  it  normally  does  not  exist  when  actual  curr  it  is  flowing.  This  causes  the  error. 
Looking  at  Fig.  5.3,  the  arbitrary  partition  yields  a  i%  error  in  the  total  resistance  even  if  the 
exact  resistance  of  the  individual  pieces  is  known. 

However,  often  a  resistance  method  applied  to  a  large  piece  that  cannot  be  partitioned 
anymore  yields  a  larger  error  in  resistance  than  the  same  method  applied  to  an  incorrect  parti¬ 
tioning.  Thus  relaxing  the  partitioning  constraint  (as  described  in  the  next  section)  so  as  to 
increase  the  partitioning  of  large  pieces  will  increase  accuracy.  For  example,  returning  to  Fig, 
5.3,  if  the  BEM  used  in  JET2  is  applied  to  the  whole  piece,  a  4%  error  occurs.  If,  however. 
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Rexact  =  ?.639fi 
Rapprox.  =  R1+R2+R3 
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Figure  5.3  Nonequipotential  partition. 
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this  method  is  applied  to  the  arbitrary  partitions  and  their  resistances  are  summed,  an  error  of 
2.6%  occurs,  which  is  actually  an  improvement. 

Partitioning  is  even  more  important  for  three  (or  more)  terminal  shapes  where  the  path  of 
current  flow  between  two  terminals  is  shared  by  .  third  terminal.  The  first  two  terminals 
become  weakly  coupled,  corresponding  to  a  large  resistance  between  them.  Also  if  the  dis¬ 
tance  between  two  terminals  is  relatively  greater  than  between  others  (as  is  commonly  the 
case  with  multiterminal  shapes),  weak  coupling  results.  Many  numerical  methods,  including 
the  BEM  used  in  JET2,  are  inaccurate  in  estimating  the  resistance  between  two  weakly  cou¬ 
pled  points.  Each  panitioning  line  in  effect  introduces  an  internal  node  which  splits  the 
segment’s  network  model  into  two  networks.  Let  us  call  these  two  new  networks  A  and  B, 
which  are  effectively  decoupled  from  each  other.  This  eliminates  tiying  to  compute  the  value 
of  the  coupling  between  the  two  weakly  coupled  terminals  as  long  as  one  of  the  terminals  is 
in  A  and  the  other  is  in  B.  Also  the  distance  between  the  new  internal  node  and  any  terminal 
in  network  A  is  smaller  than  the  distance  between  any  terminal  in  A  and  any  terminal  in  B, 
and  vice  versa.  Thus  the  coupling  between  the  internal  node  and  any  terminal  in  network  A  is 
stronger  than  the  coupling  between  any  terminal  in  A  and  any  tenninal  in  B,  and  vice  versa. 
Thus  in  JET2  we  will  relax  the  constraint  of  partitioning  along  known  equipotential  lines,  as 
described  in  detail  at  the  end  of  this  section. 

Previous  approaches  to  extracting  the  resistance  of  the  bus  segments  could  be  classified 
into  table  look-up  and  on-the-fly  numerical  methods.  Table  look-up  methods  are  faster  but 
on-the-fiy  numerical  methods  are  more  general.  Look  up  methods  generally  pattern  match  the 
partitioned  segment  with  the  closest  case  for  which  it  has  an  analytical  formula  or  table  look¬ 
up  values.  In  our  approach,  we  use  a  combination  of  the  two. 
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Many  extractors  use  table  look-up  and  various  heuristics  to  model  bus  segments.  One 
early  heuristic  was  to  model  the  comer  rectangle  of  material  in  an  L-bend  as  0.55  /?,  f49], 
where  R,  is  the  sheet  resistance.  Thus  the  L-bend  shown  in  Fig.  5.1  would  be  nuxleied  as  a 
2.55  O  resistor.  For  w,  =  Wj,  this  approximation  is  fine.  But  for  w,  >  the  resi.stance  of 
the  comer  becomes  large.  Thus  the  heuristic  breaks  down  for  this  case.  It  also  gives  .lo  indi¬ 
cation  on  how  to  model  the  four-terminal  shape  in  Fig.  5.2.  The  method  in  [50]  partitions  bus 
segments  along  equipotential  lines  and  analyzes  the  resulting  expected  simple  polygons  using 
approximate  empirical  fomiulas.  The  method  would  still  have  trouble  with  the  segment  in  Fig. 
5.2,  because  in  the  case  of  multiterminal  shapes,  it  concentrates  on  the  cunent  flow  between 
two  terminals  at  a  time,  ignoring  the  effects  that  the  other  terminals  have  on  the  current  flow, 
and  hence  the  effect  on  the  resistance  value  between  the  two  terminals  under  consideration. 
This  is  caused  by  the  lack  of  unique  equipotential  lines  for  segments  with  more  than  two  ter¬ 
minals,  as  pointed  out  earlier.  Also  as  current  crowding  effects  become  more  pronounced,  the 
heuristics  employed  in  [50]  would  break  down.  These  current  crowding  effects  are  the  same 
effects  that  cause  the  L-bend  heuristic  error  previously  mentioned.  This  is  equivalent  to  hav¬ 
ing  the  approximately  straight  equipotential  lines  that  were  originally  chosen  as  a  partition 
start  to  deviate  significantly  in  location  and  shape  as  h',/w2  changes.  While  shapes  such  as  in 
Fig.  5.2  m.ight  not  be  very  conunon,  for  the  sake  of  completeness,  a  general  resistance  extrac¬ 
tor  shcjld  be  able  to  handle  them  accurately  without  human  intervention,  such  as  modifying 
the  layout  to  suit  the  heuristics  used  in  [51]. 

The  original  JET  program  [51]  also  is  a  table  look-up  program.  Straight  equipotential 
lines  occur  at  places  of  symmetry,  such  as  between  points  A  and  B  in  Fig.  5.1  and  approxi¬ 
mately  .straight  lines  occur  one  square  away  from  bends  or  contacts  as  noted  in  [51]  and  as 
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seen  in  Ftg  5.1  Since  JET  only  haiKllcs  Manhattan  geoirtctncs.  it  patinions  one  square  away 
from  all  5T>ends  or  contacts.  Therefore,  a  segment  must  have  at  Icxsj  two  squares  of  meia! 
between  any  bends  or  contacts  for  further  partitioning.  JET  then  pattern  matches  the  piece 
exactly  to  the  shapes  and  ie,spcciivc  network  irKxk^ls  shown  in  Figs,  .5,4  ta)-(0,  if  possible.  It 
then  uses  table  look-up  vtiiiies  for  tlie  «:sistors.  These  primitive  shapes  are  the  most  com¬ 
monly  occurring  ones  in  a  VLSI  bus.  If  a  partition  docs  not  match  a  stored  pattern,  then  JET 
skips  it,  flags  an  error,  and  then  terminates.  Like  dte  method  in  [50],  JET  cannot  parse  seg¬ 
ments  like  the  one  in  Fig.  5.2.  It  cannot  even  parse  segments  the  method  in  [50]  could  handle, 
such  as  in  Fig.  5.5,  In  Fig.  5.5,  the  requirement  that  comer  bends  be  at  least  two  squares 
away  from  each  other  yields  a  nonprimitive  shape. 

Table  5.1  shows,  for  the  VDD  and  GND  busses  of  a  real  test  chip,  the  relative  frequen¬ 
cies  of  occurrence  of  the  primitives  and  nonprimitives  according  to  the  partitioning  criteria 
used  in  JET.  Note  the  significant  number  of  nonprimitives.  Table  5.2  categorizes  the  multiter- 
rninal  nonprimitives  by  the  num’ber  of  edges.  This  gives  a  measure  of  the  extent  of  partition¬ 
ing.  The  computation  time  for  large  nonprimitives  is  equal  to  that  of  .many  small  noaprimi- 
tives  because  of  the  O(V^)  running  time  of  the  resistance  calculation  algorithm  used  in  JET2. 
These  data  were  compiled  using  a  modification  of  the  new  extractor,  JET2. 
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friiijic  5  1  Occurrerjcc  oif  JET  primitives  in  typical  VDD  and  GND  busses. 
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JET?,  relaxes  rhe  partitioning  criteria  used  in  JET.  If  two  geometric  features  that  cause 
bends  in  the  current  flow  are  two  or  more  squares  away  from  each  other,  then  JET2  does 
what  JET  does.  Otherwise,  if  the  distance  between  the  two  geometric  features  is  at  least  213  of 
a  square,  then  (he  segment  is  partitioned  approximately  half  way  between  the  geometric 
features.  Table  5.3  shows  the  same  data  as  in  Table  5.2  for  the  same  chip,  except  using  this 
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Figure  5.4  jET  primitives. 
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Figure  5.4  (continued) 


new  criterion.  Note  that  the  number  of  nonprimitives  has  increased,  yet  there  are  no  computa¬ 
tionally  expensive  ones  with  a  large  number  of  edges. 


Table  5.3  Edge  frequency  for  layout  using  JET2  partitioning  criteria. 


GND  VDD 


No.  of  edges 

Occurrence 

No.  of  edges 

6 

331 

6 

1 16 

8 

537 

8 

562 

10 

255 

10 

70 

12 

65 

12 

32 

14 

34 

14 

34 

16 

37 

16 

9 

18 

2 

18 

1 

20 

i 

22 

2 

24 

1 

26 

1 

48 

2 

1 

1 


unrecogni2able  bus  segment  (nonprimitive) 


Figure  S.S  JET  partitioning  itself  into  an  unrecognizable  situation. 

S3  On>the*Fly  Numerical  Method:  The  Boundary  Element  Method 

JET2  uses  the  same  concept  of  splitting  the  bus  into  primitives  as,  JET,  but  in  addition  it 
uses  BEM  to  analyze  the  nonprimitives. 

The  potential  u  in  any  general  linear  resistive  region  (Fig.  5.6)  obeys  Laplace’s  equation 
in  tv^o  dimensions  (assuming  constant  thickness)  in  its  interior 


Fq  (insulating  boundary) 
Fu  (terminals) 


Figure  5.6  Domain  for  Laplace’s  equation. 


=  0 


(5.1a) 


the  Dirichelet  boundary  condition  on  its  conducting  boundary  (the  terminals)  is: 

u  =  it  (5.1b) 

and  the  Neumann  boundary  condition  on  its  insulating  boundary  is: 

<7=0  (5.1c) 

where  q  Buldn  and  n  is  the  unit  outward  normal  to  the  boundary.  This  follows  from 

0  =  7n  =<5En  =-o-~  (5. Id) 

on 

r„  and  need  not  be  continuous,  but  F  =  r„  +  F^.  The  current  through  any  conducting 
subboundary  F,  is  given  by 

/  =  ^  Ir  (5,2) 

where  is  the  sheet  resistance. 

There  are  a  i, umber  of  standard  methods  for  solving  Laplace’s  equation  [52],  [53].  We 
have  chosen  the  BEM  [54],  [55]  because  of  its  flex  bility,  generality,  and  computational 
efficiency. 

For  convenience,  from  now  on  the  Laplacian  operator  is  represented  by  A.  The  arc  length 
element  d  T  and  the  area  element  dil  will  be  omitted  ^'rom  integrals  since  they  are  under- 
stood.  Consider  the  same  domain  U  and  boundary  F.  Green’s  second  identity  [52]  (which  is 
derived  from  the  Divergence  Theorem)  states  that  for  any  arbitrar/  functions  u  and  v  w'hich 
have  continuous  .second  partial  derivatives  in  Q,  and  continuous  first  partial  derivatives  in  the 
combined  domain  Q  -f  F. 
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C5.3) 


Let  u  be  the  solution  to  our  PDE  as  before,  while  we  still  have  a  degree  of  freedom  in  choos¬ 
ing  V .  Using  Am  =  0  and  rearranging  (5.3),  we  obtain 


.  r  dv  f  du 

+  1  u- - j  V— 

•'r  an  •'r  an 


=  0 


Now  choose  a  function  v  such  that 


(5.4) 


Av=-8(P)  (5.5) 

where  -5(P)  is  the  2-D  Kronecker  delta  function  located  at  point  P  anywhere  in  Q.  This  is 
called  Green’s  function  (or  fundamental  solution),  otherwise,  known  as  a  delta  source.  In  2-D, 
Green's  function  is 


« 

u 


_L 

2jt 


In  r 


r  -  ^{x-XpYTJy-y^ 


(5.6a) 

(5.6b) 


where  r  is  the  absolute  value  of  the  distance  between  the  observ'ation  point  at  (x,\)  and  the 


delta  source  at  {xp,yp).  Substituting  (5.5)  into  (5.4)  yields 


,  r. ,  f  dv  r  dw 

M(P)  +  j  M- - J  V—  =0 

•'r  on  ■'r  on 

P  e  Q 


(5.7) 


This  equation  is  expressed  in  terms  of  only  the  unknown  potential  and  outward  nornurl  flux  at 
the  boundary,  except  for  the  term  u{P).  To  eliminate  m(P)  we  choose  to  take  the  point  P  of 
the  delta  source  on  F  instead  of  in  Q.  If  P  is  on  a  smooth  part  of  the  boundiiry  (i.e.,  not  on  a 
comer),  then  it  can  be  shown  [61J  that  the  area  integral  reduces  to  0.5w(/M.  If  P  sits  on  a 
corner,  then  it  reduces  to  cu{P)  where  c  dept-nds  on  the  kx;al  geometry.  It  will  be  shown 
later  that  this  constant  does  not  have  to  be  calculated. 
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If.  for  convenience,  we  replace  du/dn  ~  q,  as  we  have  done  hefore,  let  v  =  u*  (to  rem¬ 
ind  us  that  V  is  a  potential  as  well,  but  that  of  a  hypothetical  delta  source),  and  analogously 
let  du'/dn  =  we  obtain  as  a  final  contour  integral  equatioi'; 

cu  (F  )  +  uq*  ~  qu*  =  0  (5.8) 

Per 

c  -  0.5  if  the  delta  source  is  on  a  smooth  boundary. 

llie  boundary  element  equation  (5.8)  can  be  applied  as  follows:  once  we  choose  the  form 
of  the  approximating  function  for  u  and  q  on  the  bc/undary,  the  equation  gives  a  constraint  on 
the  choice  of  coefficients  of  the  approximating  function  so  as  to  minimize  the  error  between 
the  approximations  and  the  exact  answer.  If  the  coefficients  are  put  in  leims  of  values  of  u 
and  q  at  discrete  points  on  'he  boundary,  then  the  constraint  becomes  one  on  these  u  and  q 
values.  Given  enough  constraints,  we  can  form  a  system  of  equations  to  solve  for  these  u 
and  q  values. 

5.4  Discretization  of  BEM  Equations 

To  di.scretize  (5.8)  into  a  system  of  linear  equations,  we  construct  an  approximation  to  u 
and  q  by  first  dividing  the  boundary  into  N  segments  a-s  shown  in  Fig.  5.7.  Then  we  designate 
the  u  and  q  values  of  the  endpoints  of  each  segment  as  variables.  Note  that  dividing  the  boun¬ 
dary  into  N  .segments  yields  N  endpoints  and.  hence,  2N  variables  since  each  endpoint  j  has  a 
Uj  and  qj  variable  asscK'iated  with  it  (there  is  an  exception  that  will  lx;  dealt  with  later).  How¬ 
ever,  for  every  point  on  (not  ju.st  the  endpoints),  u  is  known  and  q  is  unknovv'ii.  The  exact 
opposite  occurs  for  T^.  Thus  there  will  be  only  .N  unknowns,  with  the  other  N  known  values 
combining  to  form  a  source  term.  Next,  on  each  segment,  we  approximate  the  ItKal  u  solution 
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on  the  segment  by  a  function  of  the  u  value  of  the  endpoints.  Similarly,  the  local  q  solution 
is  approximated  by  the  q  value  of  the  endpoints.  These  functions  can  be  expressed  by  interpo¬ 
lation  functions; 


u  =  Tjq(s)uj  +  Tn(s)Uj  +  I  for  segment  Jq 

q  =  TjQis)qj  +  Tji(s)qj  +  1  for  segment  Jq 

f  1  if  j  =  o] 

^>0^ lo  if5  =  1 


jo  it  J  =  0 

=  ]!  if  ^  =  1 


(5.9) 

(5.10) 

(5.11) 

(5.12) 


where  s  is  the  normalized  arc  length  along  the  seginent  as  shown  in  Fig.  5.8.  For  example,  if 
u  and  q  are  to  linearly  vary  between  the  endpoint  values,  we  obtain 


Tjq^I-s  (5.13) 

Tjx^s  (5.14) 

u  -{[—s)Uj  for  segment  y’o  (5.15) 

q  =  {l~s)qj  +sqj+x  for  segment  Jq  (5.16) 

Tills  is  called  a  linearelernent  and  is  a  common  approximating  function  used  on  boundaries 

with  Manhattan  geometries.  Higher-order  polynomial  approximations  are  possible  and  require 

k  points  per  segment  for  an  order  k  approximation. 


Figure  5,7  Division  of  boundary  into  segments. 
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Fig.  5.8  Typical  Segment 

Note  that  both  q  and  u  use  the  same  order  of  interpolation  functions  in  (5.9)  and  (5.10). 
This  is  not  necessary,  but  is  done  because  it  is  convenient.  Note  also  in  looking  at  Fig.  5.7 
that  each  point  j  i‘'  shared  between  segment  j  and  segment  j-l.  Thus  each  point  has  two 
interpolation  functions  associated  with  it. 


Next  we  choose  to  place  the  delta  source  at  point  i.  Substituting  (5.9)  and  (5.10)  into 
(5.8),  we  obtain  the  following  equation,  which  is  a  matrix  equation  since  i  can  be  any  of  the 
N  endpoints: 


£  =  i 

;=l 


H,/  =1 


rj,.  T„q,dT^ 


1 


c  + 


*  /* 


*=o 


I 

j*i 

1 

X  fink 

j=i 

L  f^ijk 

k-O 


(5.17) 


(5.18) 


(5.19) 


where  c  is  defined  in  (5.8).  As  mentioned  before,  calculation  of  c  and  hence  H,,  is  more 
complex  if  point  i  is  at  a  comer.  We  can  avoid  the  calculation  of  c  and  calculate  directly 
by  noticing  that  in  a  physical  resistor  of  any  shape,  forcing  the  boundary  condition  u  =  1  on 
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r„  makes  u  =  1  and  q  =  0  everywhere  on  F.  Thus  substituting  Uj  -  1  and  qj  =  0  into  (5.17) 
yields 

N 

=  -  X  (5.20) 

i=\,j*i 

As  mentioned  before,  because  half  of  the  2N  variables  are  known,  we  can  formally  write  the 
final  linear  algebraic  system: 

Ax  -  b 

A  =  {a,  -  I  a^i  =  Hij  if  point  j  €  r„  else  a,  =  -  G,  } 

,  T  (5-20 

X  =  {x^  I  jc,  =  u,  it  point  i  e  else  jc,  =  <?,  } 
b  =  —  H  z  where  z  =  {z,  /  !  z,  =  if  point  i  e  F^  else  z,  =  0}^ 

As  mentioned  before,  the  known  u  variables  multiplied  by  their  corresponding  column  of  H 

coefficients  are  transfeited  to  the  right-hand  side  to  form  the  source  lenn.  The  known  q  vaii- 

ables  are  not  similarly  treated  because  all  of  the  known  q  variables  are  equal  to  zero  by  the 

Neumani.  boundary  condition. 

As  described  before,  to  find  the  P(P-l)/2  resistors  for  our  representation  of  the  P  port 
resistor,  P  different  Dirichlet  boundary  conditions  are  needed.  Hence,  P  source  terms  are 

-1 

required.  Using  LU  factorization,  approximately  — +  PN^  op)erations  are  necessary  to  solve 
this  system. 

Again  the  advantage  of  the  BEM  approach  is  that  only  the  boundary  of  the  shape  is  bro¬ 
ken  up,  not  the  interior.  Thus  the  generated  matrix  sizes  are  much  smaller,  esjrecially  for  com¬ 
plex  shape.s,  than  those  for  finite  element  method  (ITiM)  [52].  Also  generation  of  grid  points 
on  a  boundary  is  much  simpler  than  the  generation  of  triangular  meshes  m  the  interior,  as  is 
done  in  FEM.  Also,  the  outward  norma!  tlu.x  q  is  solved  for  directly  in  BtiM,  whereas  in 
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traditional  F’EM  one  has  to  subtract  two  nearby  voltage  values. 

5.4 J  Extension  to  Multiply-Connected  Domains 

So  far  we  have  considered  simple  domains,  as  in  Fig.  5.6,  where  the -boundary  i.s  con¬ 
tinuous  and,  hence,  the  contour  used  for  the  contour  integral  is  straightforward.  Fcr  more 
complex  regions  with  internal  contacts  or  holes,  as  shov/n  in  Fig.  5.9(a),we  can  visualize  the 
contour  integral  path  (dotted  line)  starting  at  the  outside  boundary,  covering  part  of  it,  then 
crossing  the  interior  to  cover  the  two  holes  in  the  manner  shown,  then  doubling  back  on  itself 
to  cover  the  rest  of  the  outer  boundary.  Those  segments  where  the  path  crosses  over  itself 
have  outward  normals  that  point  in  opposite  directions  as  indicated.  Looking  back  at  the  con¬ 
tour  integrals  in  (5,9),  we  see  that  these  segments’  contributions  cancel  out,  resulting  in  the 
contours  shown  in  Fig.  5.9(b).  Note  that  the  two  internal  contours  are  counterclockwise  Vv-hile 
the  outer  one  is  clockwise.  This  will  always  Lc  the  case.  Thus,  the  net  contour  integral  of  a 
disjoint  boundary  equals  the  sum  of  the  contour  integrals  of  the  disjoint  boundaries  with  each 
internal  boundary  evaluated  in  a  direction  opposite  to  the  outer  hoimdary. 

5.4.2  Special  Case  of  Points  with  Two  Unknowns 

Earlier  it  was  mentioned  at  the  beginning  of  Section  5.4  thaf  each  point  has  one  unk¬ 
nown  and  one  known  variable  associated  with  it.  The  exception,  shown  in  Fig.  5.10,  occurs  at 
comers  in  the  Dirichlet  boundary,  u  is  known,  but  there  are  two  unknown  q  variables 
since  there  are  two  outward  normals  and  the  outward  norma!  flux  is  not  zero  for  both  of  them. 
The  authors  of  [55 J  describe  two  ways  of  dealing  with  this  problem.  The  first  is  the  partially 
discontinuous  element.  In  this  method,  two  separate  but  close  nodes  are  used  to  model  the 
corner.  This  preserves  the  atle  that  each  point  has  only  one  unknown  associated  v  ith  it.  A 


(a) 


(b) 


Figure  5.9  Multiply  connected  regions  and  their  contours. 

(a)  total  contours  (b)  net  contours 

more  accurate  method,  adopted  in  [55]  and  adopted  by  us,  is  the  double  node  where  the 
comer  node  is  treated  as  having  two  unknowns  and  qjQ.  The  G^j  coefficient  of  a  normal 
node  is  made  of  the  sum  of  the  two  contributions  from  the  two  segments  to  its  immediate  left 
and  right,  as  in  (5.19);  here,  the  coefficient  for  is  only  the  contribution  from  the  segment 
j  1  and  similarly  for  qjQ.  In  other  words. 


^lyl 


(5.22) 


(5.23) 


Because  of  the  extra  unknown,  an  extra  equation  will  be  needed.  This  can  be  formed  by  plac¬ 
ing  the  delta  source  anywhere  else  on  the  boundary,  except  at  the  segment  endpoints  (where 
ail  the  previous  delta  .sources  have  been  placed).  As  in  [55],  we  choose  to  place  it  at  the  mid¬ 
point  of  segment  j  1 .  Unlike  the  other  N  equations,  the  potential  and  flux  at  the  delta  source 
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Figure  5.10  Double-node  condition. 

point  are  not  considered  unknowns  to  be  solved  for  since  by  the  very  definition  of  the  linear 
approximation  we  are  using,  these  quantities  are  only  the  arithmetic  mean  of  those  at  the  end¬ 
points  of  the  segment.  Thus  the  //,,  coefficient  does  not  have  the  same  meaning  for  the  extra 
equation  as  for  the  other  normal  N  equations.  Because  the  delta  source  is  on  a  smooth  sur¬ 
face,  the  c  tenn  in  (5.8)  is  0.5.  Since  the  delta  point  has  no  unknown  of  its  own,  and  hence 
no  coefficient  of  its  own  to  add  this  c  correction  term,  we  can  add  it  to  the  H  coefficient 
associated  with  one  of  the  segment  endpoints,  i.c.,  //,y  or 

5.4.3  Properties  of  the  Resulting  Linear  System 

The  flux  and  potential  produced  by  a  delta  source  placed  at  one  piiit  of  the  boundary  are 
nonvanishing  at  every  other  part  of  the  boundary.  Thus,  any  point  couples  to  every  other 
point  and  the  resulting  matrix  A  defined  in  (5.21)  is  full.  The  matrix  is,  in  general,  asym¬ 
metric  since  the  //  and  G  matrices  of  which  it  is  compri’^ed  are  asyronietric  as  well,  because 
the  H  and  G  coupling  coefficients  between  two  Doint.s  do  not  depend  only  on  the  distance 
between  the  points  (which  is  the  same  no  matter  which  point  has  the  delta  source  placed  on 
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it).  The  coefficients  also  depend  on  the  length  and  orientation  of  the  segments  the  points  are 
on.  If  these  aie  not  symmetric,  then  neither  is  the  matrix.  Also  the  system  is  not  guaranteed  to 
converge  using  relaxation  methods,  thus  requiring  the  full  LU  factorization  to  solve  the  sys¬ 
tem. 

Thus,  even  though  the  BEM  system  of  equations  is  significantly  smaller  than  the 
conesponding  FEM  or  the  finite  difference  method  (FDM)  [52]  system  which  yields  the  same 
accuracy,  if  does  not  have  the  advantages  of  sparsity  and  symmetry  that  the  FEM  system 
does.  Also  the  computation  of  the  coefficients  in  BEM  are  more  expensive  than  in  FEM  or 
FDM.  Thus  the  question  remains  which  system  is  faster  to  solve.  Comparison  in  [55]  of  the 
computation  time  between  both  BEM  and  FDM  using  the  same  shapes  has  shown  that  the 
BEM  system  is  still  faster,  by  factors  as  large  as  3  to  7.  The  more  complicated  the  shape,  in 
fact,  the  larger  the  savings.  A  simple  FEM  program  written  as  part  of  this  work  has  shown 
that  FHM  with  uniform  grids  to  be  slower  than  BEM  as  well  for  the  same  accuracy,  even  with 
sparsity  and  symmetry  of  the  system  taken  advantage  of. 

In  [55],  for  simple  shapes  most  of  the  computation  time  was  due  to  computing  the  rriatrix 
coefficients.  The  authors  evaluate  the  integrals  in  (5.18)  and  (5.19)  numerically,  probably  by 
some  special  techniques  because  the  kernels  of  the  integrands  becouie  singuiai'  at  the  delta 
source  point.  In  our  approacii,  we  try  to  decrease  computation  time  and  improve  accuracy  by 
calculating  the  coefficients  analytically. 

5.4.4  Analytical  C'aiculation  of  Matrix  H  and  Matrix  G  Coefficients 

To  calculate  the  contribution  of  a  segment  to  //,,  and  ,  without  loss  of  generalization, 

we  shift  the  coordinate  system  so  that  the  delta  source  point  i  is  at  the  origin  and  roiate  the 
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axes  so  that  the  segment  with  point  j  as  its  endpoint  is  vertical  and  assume  that  the  contour 
integration  direction  is  the  positive  y  direction  (c/T  =  dy).  We  also  assume  that  the  outward 
norma!  n  is  Point  j  is  located  at  _y|  and  point  j  +  1  at  ^2  with  x  as  the  directed  distance 
l)etween  the  delta  source  and  the  segment  (Fig.  5.11).  If  the  actual  contour  direction  is  in  the 
negative  y-direction,  then  and  G,y  are  negated.  The  //y,-  is  also  negated  if  the  actual  out¬ 
ward  normal  in  the  layout  is  . 

The  interpolation  function  Tjq  in  (5.35)  and  (5.36)  and  u  *  and  q  *  in  terms  of  this  new 
coordinate  system  are 
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(5.25) 

(5.26) 


Bn  dx  2rcin  (x'-fy  ^■) 

Substituting  the  above  into  the  contributions  //,yo  anti  G,y()  of  (5.18)  and  (5.)  9),  respectively, 

yields: 
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where  terms  have  been  grouped  as  much  as  possible  to  reduce  the  number  of  evaluations  of 
logarithmic  functions  first  and  then  the  number  of  floating  point  divides  and  multiplies.  The 
singularities  of  the  kernel  are  integrated  out. 

Note  that  for  certain  zero  values  of  x,  y 2  oc  yi  the  above  functions  can  be  simplified.  To 
save  computation,  we  check  for  these  conditions  and  use  intermediate  variables  to  eliminate 
redundant  operations. 

The  contribution  of  the  other  segment  (with  points  j  and  j~l  as  its  endpoints)  yields 
and  G,yi  and  is  calculated  much  the  same  way,  except  that  Tji  is  used  as  the  interpola¬ 
tion  function  since  point  j  now  fulfills  the  same  role  for  this  segment  as  poi.it  j  +  l  did  for  the 
first  segment.  In  fact,  if  one  designates  the  coordinate  of  point  j—l  as  yg  redefines  x  with 
respect  to  the  new  segment,  one  finds  that 

^^7 1(^  ')  =  ,y  1  .yo)  (5  29) 

Gij  I .y  o.y  1)  “  -y  i  -y  o)  (5-30) 

Thus  one  function  can  be  used  to  evaluate  both  parts  of  and  similarly  for  G,y . 
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Figure  5.11  Setup  for  analytical  calculation  of  H,j  and  G,y  coefficients. 
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5.4.5  Current  Calculation 


After  the  system  is  solved,  all  of  the  q’s  are  known  and  can  be  used  to  solve  for  the 
cunent.  The  current  through  any  tenninal  equals  the  sum  of  the  currents  through  the  segments 
that  compose  that  terminal.  To  find  the  current  through  the  segment  in  Fig.  5.11,  let  <71  and 
^2  be  the  solved  values  at  yi  and  respectively.  Then  on  the  segment 

9= -  (5.31) 

Substituting  this  in  the  definition  of  /  in  (5.2)  yields 
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In  summary,  in  this  section  we  have  derived  the  boundary  element  method  which  solves 
Laplace’s  equation  by  first  transforming  it  from  a  statement  about  the  potential  over  an  area 
domain  into  an  equation  on  the  boundary  of  that  domain  Then  this  resulting  contour  equation 
is  approximately  solved  by  discretizing  it  into  a  linear  algebraic  system  using  linear  approxi¬ 
mations  to  the  potential  and  flux  on  the  boundary.  The  method  was  further  generalized  for 
multiply  connected  regions.  Solving  the  linear  system  yielded  the  outward  normal  flux  at 
discrete  points  on  the  boundary,  which  were  then  used  to  calculate  the  current  flowing  through 
the  terminals  in  the  original  problem. 

5.5  Comparison  with  Finite  Element  Method  (FEM) 

In  this  section  we  compare  BEM  with  a  FEM  method  that  directly  computes  resistance. 
The  FEM  will  be  used  later  in  conjunction  with  BEM  to  extract  the  capacitance  model. 


124 


Again,  in  this  section,  we  assume  normalized  =  1  and  normalized  =  1 . 

In  [56]  a  method  of  replacing  a  triangle  of  an  FEM  mesh  with  an  equivalent  three- 
terminal  delta  RC  network  and  a  rectangle  with  its  four-terminal  RC  equivalent  is  derived. 
The  rectangle  is  of  interest  to  us  since  the  layout  shapes  considered  in  this  thesis  are  all  rectil¬ 
inear  due  to  the  Manhattan  geometry.  The  rectangular  mesh  model  and  its  equations  are 
shown  in  Fig.  5.12. 


R 1  =  2x/y 

(5.33a) 

R2  =  2y/x 

(5.33b) 

C=xy 

(5.33c) 

Replacing  all  of  the  mesh  rectangles  with  their  equivalent  delta  networks  yields  a 
detailed  RC  network  model  of  the  segment  being  modeled.  Considering  the  resistance  network 
alone,  by  collapsing  all  nodes  on  one  terminal  into  a  single  terminal  node  and  using  delta-Y 
transformation  of  to  eliminate  internal  nodes  yield  the  equivalent  star  model.  In  Section  5.6 
below,  we  show  how  this  method  can  be  used  in  conjunction  with  BEM  to  extract  the  capaci¬ 
tance  of  bus  segments.  However,  for  now,  we  consider  it  to  be  representative  of  FTM 
methods  to  do  a  comparison  with  BEM  for  resistance  extraction. 


Figure  5.12  RC  equivalent  of  rectangular  mesh  element. 
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now,  we  consider  it  to  be  representative  of  FEM  methods  to  do  a  comparison  with  BEM  for 
resistance  extraction. 

Shown  in  Fig.  5.13  are  shapes  for  which  the  author  of  [56]  uses  the  above  FEM  method 
with  an  adaptive  mesh  to  calculate  the  resistance.  The  resultant  absolute  percentage  enror  in 
resistance  vs.  number  of  nodes  is  plotted  next  to  each  shape.  The  performance  of  BEM  using 
a  simple  uniform  grid  is  also  plotted.  Clearly  BEM  outperforms  FEM  since  for  every  desired 
level  of  accuracy,  BEM  requires  a  significantly  fewer  number  of  nodes,  even  without  the 
benefit  of  an  adaptive  mesh. 

5.6  JET2  Capacitance  Mode!  Extraction 

Iti  JET2,  a  bus  segment  is  simply  modeled  as  a  parallel  plate  capacitor  to  determine  the 
total  capacitance  to  ground.  JET2  improves  over  JET  in  distributing  this  capacitance  among 
the  terminal  nodes  of  »he  segment’s  RC  model.  The  next  section  describes  how  a  combination 
of  the  finite  element  method  and  a  special  node  reduction  technique  determines  the  distribu¬ 
tion  of  capacitance  among  the  terminals  so  as  to  accurately  model  the  transient  response  of 
the  bus  segment.  For  primitives  shown  in  Fig.  5.4,  the  resulting  capacitance  to  ground  at  each 
node  is  stored  in  a  look-up  t.able.  F'or  nonprimitives,  the  capacitance  model  extraction  is  done 
on-the-fly. 

We  first  divide  the  bus  segment  into  many  rectangles.  Then  the  FEM  method  described 
in  Section  5.5  is  used  to  replace  each  rectangle  with  a  simple  RC  model  (Fig.  5.12).  This 
yields  a  detailed  RC  network  for  the  whole  bus  segment  which  accurately  models  its  distri¬ 
buted  RC  effects.  Note  that  since  each  rectangular  element  has  a  total  capacitance  to  ground 
equal' to  its  parallel  plate  capacitance,  the  sum  of  all  the  capacitances  m  the  overall  segment’s 
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RC  model  is  equal  to  its  parallel  plate  capacitance  as  well. 


The  FEM  method  usually  produces  a  very  large  RC  model  for  the  segment,  depending 
on  the  number  of  grids  chosen  to  model  the  segment.  We  would  like  to  simplify  the  network 
into  a  delta  model  while  still  retaining  the  essence  of  its  transient  response.  Reference  [56] 
presents  an  RC  network  node  reduction  method  based  on  Elmore  time  constants  [57]  that  does 
just  that.  The  final  simplified  network  has  the  following  properties: 

1)  The  Elmore  delay  between  any  two  terminal  nodes  is  the  same  as  that  of  the  original 
network. 

2)  The  terminal  resistance  parameters  of  the  network  are  preserved. 

3)  No  additional  coupling  capacitances  are  introduced. 

•>)  The  total  capacitance  to  ground  of  the  final  network  is  equal  to  that  of  the  original  net¬ 
work. 

Thus  the  transmission  properties  of  the  detailed  RC  network  are  accurately  preserved.  The 
node  elimination  formulas  [56]  for  networks  that  contain  no  coupling  capacitance  are 


8‘u  =  Sij  +  — ^ 
^  Smk 


(5.34a) 


Cl,  Sij, 

C',  =  C.  +  ~ 


Smk 
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(5.34b) 


where  the  above  is  for  the  elimination  of  node  k.  All  of  the  above  four  properties  are 
preserved  after  each  node  elimination  and  hence  are  preserved  in  the  final  network.  After 
elimination  of  all  internal  nodes,  a  delta  model  is  obtained  with  all  the  capacitances  still  con¬ 
nected  to  ground. 
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Shown  in  Fig.  5.14  is  a  comparison  of  different  RC  models  for  the  unit  step  response  of 
the  T-junction  primitive  of  W2/W,  =5  (see  Fig.  5.24  below  for  meaning  of  W|  and  1V2)  for 
zero  initial  conditions.  The  input  terminal  is  one  of  the  symmetric  terminals  with  the  other  ter¬ 
minals  floating.  The  output  terminal  is  the  asymmetric  one.  The  most  accurate  model  is 
"fern"  v/hich  uses  a  rectangular  mesh  of  400  nodes.  Thus  we  want  the  other  two  simpler 
models  to  approximate  the  response  of  this  first  model.  The  other  two  models  are  "eimore," 
which  uses  the  above  node  reduction  technique  to  reduce  the  model  to  three  terminals,  and 
"simple,"  also  reduces  the  model  to  three  terminals  with  the  same  resistance  values  as 
"eimore,"  bs't  with  the  total  capacitance  C  equally  distributed  to  each  node  as  C/3.  The  figure 
shows  that  the  Elmore  model  is  very  accurate. 

Because  the  Elmore  node  reduction  technique  does  not  alter  the  multiport  resistance 
parameters  of  the  network,  one  can  determine  them  from  some  other  method  that  converges 
faster  to  the  correct  resistance  values,  such  as  BEM.  Hence  the  resistances  from  the  simplified 
network  are  ignored,  and  the  resistance  values  obtained  by  BEM,  which  are  more  accurate, 
are  retained.  The  lumped  capacitances  at  each  teiininal  of  the  simplified  network  are  kept  and 
used  for  the  capacitances  of  our  model.  This  hybrid  approach  is  used  in  JET2,  as  explained 
next. 

5.7  RC  Models  of  Primitives 

The  rectangular  segment  and  the  L  primitive  have  analytical  resistance  models.  The 
resistances  for  the  other  primitives  are  stored  in  tables.  For  ail  of  the  primitives,  except  for  the 
rectangular  segment,  the  fraction  of  ground  capacitance  distributed  to  a  terminal  node  is  stored 


in  tables. 


5.7.1  Rectangular  Segment 


The  model  of  Fig.  5.15,  where  R  =  L/W  and  C  is  the  total  ground  capacitance,  has  the 
same  Elmore  delay  as  a  detailed  distributed  RC  model  [56]. 

The  added  improvement  in  this  model  is  that  C  now  takes  into  account  fringing  capaci- 
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tance  as  well  as  parallel  plate  c^ipacitance.  Using  the  model  in  Fig.  5.16,  from  [58]  we  have 
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where  t  thickness  of  the  metal,  h  =  height  above  the  ground  plane,  and  /  =  length  of  the 
resistor,  and  the  metal  is  surrounded  by  a  material  of  dielectric  constant  £.  The  first  term  in 
the  bracket  is  the  normal  parallel  plate  contributio  i  and  the  second  term  is  due  to  the  fringing 
field.  The  above  formulas  are  accurate  except  whtn  t<Ji.  If  t  and  h  are  not  specified,  JET2 
uses  only  the  parallel  plate  model  and  .  For  the  case  of  metai2,  which  might  be  embedded 
in  a  different  dielectric  than  metal  1  (Fig.  5.17),  the  fringing  term  in  (5.35)  is  multiplied  by  a 
correction  factor  of  (l+e2''£i)/2  [13]. 


T  T 

C/2  C/2 


Figure  5.15  Rectaigular  segment  model. 
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Figure  S,16  Ground  capacitance  model  for  rectangular  segment. 
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Figure  5.17  Differing  dielectrics. 


5.7.2  L  Primitive 


In  the  model  shown  in  Fig.  5.18,  /?  is  given  analytically  in  [59]  as 


2 


a 


2  ,  ,  4a  .  a^-l  „i .  u  ~i  . 

—  In  — )  + - cos  — ) 

K  7W  aVi 


niax(w  j.vvo) 

a  -  — - >  1  (5.36) 


mj«(Wi,W2) 

where  width  ratio  a  reflects  /?(h'i,W2)  — /?(w2,w  j).  The  resistance  R  versus  a  is  plotted  in 
Fig.  5.19.  The  total  capacitance  to  ground  C  for  this  and  the  rest  of  the  primitives  is  derived 
from  the  parallel  plate  model.  In  this  case,  the  total  capacitance  is: 

C  -  C„{wf  +  +  w,W2)  (5,37) 

The  plot  of  the  normalized  Elmore  capacitance  C|/C  of  the  smaller  width  terminal  versus 
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Figure  5.20  Cl/C  vs.  widtli  ratio  (L  bend). 

integer  values  of  a  is  sliovn  in  Fig.  5.20.  They  are  computed  by  using  the  RC  network  node 
reduction  method  mentioned  abov,?  applied  to  a  fine  rectangular  mesh.  These  normalized 
capacitances  ai'c  stored  in  a  look-up  table  indexed  by  the  width  ratio.  Normalized  capaci¬ 
tances  for  other  width  ratios  are  determined  by  Linear  interpolation  and  extrapolation.  This 
table  look-up  method  is  also  used  for  the  rest  of  the  primitives  in  determining  the  Elmore 
capacitance.  Tlie  capacitances  of  the  other  terminals  in  tfie  primitives  are  determined  by  sym- 
metiy  and  the  fact  that  all  of  the  Elmore  capacitances  add  up  to  C . 


5.7 J  T  Primitive 

The  thrce-teirninal  model  is  shown  in  Fig.  5.21.  Note  that  because  of  symmetry,  only 
two  conductances  have  to  bi  computed.  Pie  conductances  are  computed  for  integer  values  of 
and  (plotted  ic  fug.  5.22)  using  BEM  witL  a  very  high  numbe'-  of  nodes.  These 
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conducta^ces  are  stored  in  a  look-up  table  indexed  by  the  width  ratio.  Conductances  for  other 
ratios  are  linearly  interjiolated  or  extrapo'atcd  from  table  values.  This  look-up  table  method  is 
also  used  for  the  rest  of  the  primitives.  The  total  capacitance  is 

C  =  Cojt(2  wf +w|  +W1W2)  (5.38) 

and  the  aormalized  Elmore  capacitance  for  the  symmetric  temiinal  is  plotted  in  Fig.  5.?.3  for 
integeis  W2/W1  and  wilw^.  'rhis  priEnitivc  is  an  example  of  the  hybrid  BEM-FEM  approach 
since  the  conductances  are  precomputed  by  BEM  and  the  capacitances  are  precomputed  by 
FEM. 


D.OO  100  4  00  4.00  f.OO  10  00 

JRJL 

Figure  5.22  T-section  conductance  vs. width  ratio. 
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Figure  5.23  Cl/C  vs.  width  ratio  (T  junction). 

5 J.4  Four-Way  Junction  Primitive 

'^The  four-terminal  model  is  .shown  in  Fig.  5.24.  Because  of  symmetry,  only  three  conduc¬ 
tances  have  to  be  computed.  Actually  in  the  look-up  table  only  two  conductances  have  to  be 
stored  (G1  and  G2)  because  G3{w2,w,)  =  G2(wj,w2)  due  to  further  symmetry.  The  conduc¬ 
tance  versus  width  ratio  is  plotted  in  Fig.  5.25.  The  capacitance  is 

C  =  wf  +  -t-  WiW2)  (5.39) 

and  the  normalized  Elmore  capacitance  of  the  smaller  terminal  is  plotted  in  Fig.  5.26. 
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Figure  5.24  Four-way  junction  model. 


Figure  5.25  Four-way  conductance  vs.width  ratios 


Figure  5.26  Cl/C  vs.  width  ratio  (four- way). 

S.7J  Contact  Primitive 

The  three-terminal  model  is  shown  in  Fig.  5.27.  Because  of  symmetry,  only  two  conduc¬ 
tances  have  to  be  computed.  Besides  being  dependent  on  the  width  ratio,  the  conductance  is 
also  a  function  of  how  far  the  contact  is  off  center,  measured  by 
shift  =  nun(a/(w-n'^),l-a/(w-Wc)),  where  tlie  meaning  of  a  is  shown  in  Fig.  5.27.  Shift 
thus  varies  from  0  to  0.5.  The  conductance  versus  width  ratio  and  shift  is  plotted  in  Fig. 
5.28.  The  total  capacitance  is 

C  =  C^x (5-40) 
and  the  normalized  Elmore  capacitance  of  the  center  terminal  is  plotted  in  Fig.  5.29.  Like  the 
conductance,  the  normalized  Elmore  capacitance  is  a  function  of  the  shift  parameter  as  well. 


138 


Figure  5.29  Cl/C  vs.  width  ratio  (coiitacf  prinitive). 

5.8  Algorithms  and  Data  Structures  of  JET2 
5.8.1  Structure  of  JET2 

The  basic  structure  of  JET2  is  shown  in  Fig.  5.30  and  a  more  detailed  algorithm  is  given 
below; 
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Figure  5.30  Basic  structure  of  JET2. 


Main  JET2  algorithm 

{Read  technology  parameters; 

read  in  metal  1 ,  metal2.  diffusion  contact,  vias 

duplicate  via  list 

for  each  metal  layer  n 

(break  metal  where  overlaps  via  and  diffusion 
merge  contacts  and  metal  lists 
put  layout  in  maximal  vertical  strip  format 
partition  layout  horizontally 
put  layout  in  maximal  horizontal  strip  format 
partition  layout  vertically 
assign  node  number  to  equipotential  lines 
identify  primitives  and  output  RC  model 
for  each  nonprimitive 
{determine  boundary 
do  rough  FEM  analysis 
do  BEM  cn  boundary 
output  RC  model 
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A  file  containing  a  list  of  rectangles  specifying  metal  1,  meial2  and  vias,  which  comprise 
the  bos  as  well  as  the  diffusion  contacts  that  connect  the  bus  to  transistors,  is  read  into  a  data 
structure  of  rectangles.  Vertical  and  horizontal  equipotenii.al  lines  are  selected  by  examining 
each  rectangle  to  see  if  it  can  be  split  according  to  the  partitioning  criterion  given  in  Section 
5.2.  During  this  partitioning  step,  horizontal  and  vertical  simple  resistors  are  identified,  since 
unlike  other  primitives,  they  are  composed  of  only  one  rectangle  and  we  already  examine 
each  rectangle  individually  in  the  partitioning  step.  A  horizontal  resistor  has  horizontal  current 
flow  and  is  bounded  by  vertical  equipotential  lines.  A  vertical  resistor  has  vertical  current  flow 
and  is  bounded  by  horizontal  equipotential  lines.  Node  numbers  are  then  assigned  to  the  edges 
of  rectangles  that  coincide  with  these  partitions.  Next,  starting  with  the  first  rectangle  i  in  the 
data  structure,  we  start  accumulating  rectangles  that  touch  this  rectangle  and  each  other.  Accu¬ 
mulation  is  stopped  when  the  group  of  abutting  rectangles  topologically  match  a  primitive  or 
when  it  is  impossible  for  the  group  to  match  any  primitive.  Both  possibilities  are  determined 
by  a  decision  tree.  If  a  group  does  match  a  primitive,  it  is  removed  from  the  data  structure 
and  has  its  RC  model  outputted.  Regardless  of  matching,  we  go  to  rectangle  i  +  1  and  repeat. 
Once  all  primitives  have  been  identified,  the  data  stmcture  is  scanned  again  for  groups  of  rec- 
tangle.s  that  touch  each  other  and  are  bounded  by  equipotential  or  insuiating  edges.  The  boun¬ 
dary  of  this  group  is  extracted  and  put  into  a  data  structure  of  points  and  the  rectangles  irre 
removed  from  the  data  strucmre.  The  structure  of  points  is  analyzed  by  a  boundary  element 
routine  which  outputs  the  structure’s  RC  model.  This  is  done  for  each  group  of  rectangles 
until  the  original  data  stmcture  has  no  more  noncontact  rectangles  left. 

Once  this  whole  proce.ss  is  done  for  metal  I,  it  is  repeated  for  rnetal2  The  two  circuit 
netlists  for  metal  1  and  metal2  are  connected  by  the  node  numbers  of  the  vias,  wliicb  arc  given 
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in  the  input  file  to  JET2.  Thus  the  vias  are  implicitly  treated  as  perfect  nodes;  hence,  current 
crowding  effects  aiound  the  vias  are  ignored.  Accurate  modeling  of  \  ias  still  has  to  be  done. 

The  main  data  structures  and  algorithms  used  in  the  program  will  now  be  presented. 

5.8.2  Rectangular  Data  Structure 

The  layout  data  for  a  given  layer  are  read  into  a  doubly  linked  list  of  rectangles  whose 
elements  are  of  the  form: 

rtruct  rectlist 

( int  rx,  uy,  lx,  dy; 
char  right,  left,  up,  down; 

short  node  up,  node  down,  node_right,  node  left; 

chai'  type; 

boolean  visited; 

struct  rectlist  *next,  *prev; 

} 

rx  and  uy  are  the  x  and  y  coordinates,  respectively,  of  the  upper-right  comer  of  the  rec¬ 
tangle.  lx  and  dy  are  coordinates  of  the  lower-left  comer,  righi,  left,  up  and  down  mark  each 
of  the  four  edges  of  the  rectangle  as  being  in  1  of  3  states:  INSULATOR,  EQUIPOT,  or 
TOUCH.  INSULATOR  means  the  edge  has  no  current  flowing  through  it,  i.e.,  it  does  not 
touch  any  other  rectangle.  EQUIPOT  means  the  edge  is  an  cquipotential  line.  This  occurs 
when  either  the  edge  is  completely  engulfed  by  a  contact  or  was  created  by  partitioning. 


Figure  5.31  Edges  with  stattis  of  TOUCH. 
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TOUCH  refers  to  any  edge  not  covered  by  the  previous  states  (Fig.  5.31).  A  "TOUCH"  edge 
indicates  that  current  bends  as  it  passes  through  that  edge. 

node  up,  node  down,  node  rlght,  node  left  apply  to  their  respective  edges  if  that 
edge’s  status  is  EQUIPOT.  They  specify  the  node  in  the  final  RC  model  that  their  edge 
corresponds  to.  If  the  edge’s  status  is  not  EQUffOT,  then  the  value  of  node  **  is 
NO_NODE_NUM.  These  node  numbers  are  either  already  assigned  to  the  edge  during  data 
structure  initialization  if  their  rectangle  is  listed  as  a  diffusion  contact  in  the  input  file  or  they 
are  assigned  after  partitioning  as  mentioned  in  Section  5.8.1. 


type  takes  on  the  values  CONTACT,  SIMPLE,  and  REGULAR  which  indicate  if  the 


whole  rectangle  is  either  a  contact,  simple  resistor,  or  neither.  This  is  necessary  since  contact 


rectangle  lists  are  merged  with  a  metal  list,  visited  (TRUE  or  FALSE  indicates,  during  a 


recursive  traversal  of  all  the  rectangles  in  the  list,  whether  the  rectangle  has  been  processed 


already,  next  and  prev  are  pointers  to  the  ne.xt  and  previous  rectangles  in  the  list.  The  first 


and  last  elements  in  the  list  do  not  correspond  to  any  rectangle.  They  arc  just  a  header  and 


tail.  'Fhe  data  structure  field  values  for  a  shape  are  shown  in  Fig.  5.32. 


IIIH 

m 

G 

■ 

R  model 
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V  rx4008  right  TOUCH 

\uy  508  left  EQUIPOT 

lx  500  up  INSULATOR 

dy  4000  down  INSULATOR 

node_up  NO_NODE_NUM 
node_down  NO_NODE_NUM 
nodejeft  5 

node_right  NO_NODE..NUM 
type  REGULAR  ” 


©■ 


Figure  5.32  Values  of  ficld.s  in  realist  data  structure  for  typical  rectangle. 
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S.8.3  Maximal  Vertical  and  Horizontal  Strip  Formats 

To  horizontally  partition  the  layout  with  vertical  equipotential  lines,  it  is  necessary  to 
first  put  the  metal  layout  into  maximal  vertical  strip  (MVS)  format.  The  two  properties  of 
MVS  format  are:  1)  all  the  rectangles  of  one  layer  are  as  vertically  tall  as  possible  (i.e.,  no 
rectangle  touches  any  other  rectangle  of  the  stime  layer  on  the  top  or  bottom  edge)  and  2) 
these  vertical  rectangles  are  as  wide  as  possible.  Figure  5.33(a)  shows  a  metal  layout  in 
MVS  format.  Figure  5.33(d)  shows  the  same  iayout  in  MVS  format.  The  algorithm  for  put¬ 
ting  an  arbitrary  layer  of  rectangles  (given  by  the  data  structure  in  the  previous  section)  into 
MVS  format  [51]  is  given  below  with  the  intermediate  steps  shown  in  Fig,  5.33(a)-(d): 


(a)  InilJally  not  MVS  (b)  split  veitically  and  horizontally  (c)  combine  vertically  (d)  combine  horizontally 

Figure  5.33  Step;;  in  MVS  algorithm. 


Procedure  MVS 

1.  For  each  rectangle,  split  current  rectangle  into  left  and  right  pieces  if  either  left  or  right 
edge  of  another  rectangle  is  within  the  left  a.id  nght  liimts  of  current  rectangle. 


2.  For  i^ach  rectangle,  split  current  rectangle  into  upper  and  lower  pieces  if  either  upper  or 
lower  edge  of  another  rectringle  i^^  within  upper  and  lower  limits  of  current  rectangle. 

3.  Combine  rectangles  in  y  -direction  as  much  as  possible. 

4.  Combine  rectangles  in  x  -direction  as  much  as  poss.ibIe. 


(a)  Improper  (b)  Proper 

Figure  5.34  Exception  to  Step  4  of  MVS  algorithm. 


Property  1  maximizes  the  vertical  cross-sectional  area  that  the  horizontal  cunent  flows 
into  and  ensures  that  the  top  and  bottom  edges  are  not  touched  by  any  other  rectangle  of  the 
same  layer.  Property  2  ensures  that  the  left  and  right  edges  of  the  rectangle  are  bordered 
either  by  an  insulating  or  equipotential  region  or  a  rectangle  of  the  same  layer  that  causes  a 
bend  in  the  current  flow.  The  rnly  problem  is  that  a  bend  in  the  current  flow  can  also  be 
caused  by  a  contact,  which  is  a  rectangle  in  another  layer,  as  .shown  in  Fig.  5.34(a).  Thus 
step  4  of  the  MVS  algorithm  where  the  rectangles  are  combined  horizontally  if  possible  has  to 
be  modified  to  take  care  of  this  case.  The  correct  partition  is  shown  in  5.34(b)  Now  the  top 
and  bonon.  edges  of  each  rectangle  are  gu:iranteed  not  to  have  a  status  of  TOUCH. 

After  M  VS  foimatting,  .my  rectangle  with  a  top  and  bottom  edge  status  of  INSULATOR 
is  ;*  candidate  for  partitioning  by  vertical  lines  as  the  various  cases  of  Fig.  5.35  show. 
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Figure  5.35  Partitioning  of  variou.?  cases. 

Vertical  partitioning  uses  maximal  horizontal  strip  (MHS)  fonnat,  which  is  completely 


analogous  to  MVS.  Tlie  only  difference  is  that  since  MHS  formatting  is  done  iifter  MVS  for- 
ntatting  and  horizontal  partitioning,  it  is  not  allowed  to  split  or  combine  iiny  horizoiral  resis- 
t‘  rs,  thus  undoing  previous  work. 


5.8.4  Primitive  Identification 

As  mentioned  before,  primitive  identification  is  done  on  the  fly  as  cdjacent  rectangles  are 
accumulated.  This  saves  time  by  avoiding  the  further  search  of  adjacent  rectangles  when  the 
accumulated  group  is  clearly  not  a  primitive  any  longer.  Searching  for  a  rectangle  that 
touches  the  current  one  is  an  O  process  in  the  average  case  for  a  layout  of  N  rectangles 
and  is  an  OiN)  process  in  the  worst  case  and  should  be  avoided  when  possible.  A  binary 
decision  tree/aJgorithin  does  this  identification  for  the  noncontact  primitive  topologies  in  MHS 
format  given  in  Fig.  5.36.  Each  added  rectangle  causes  a  branch  in  the  tree  to  be  taken  and 
possible  matches  to  certain  primitives  to  be  eliminated.  If  a  match  is  found,  the  group  of  rec¬ 
tangles  is  removed  from  the  list  and  matching  continues  with  the  next  rectangle  on  the  list. 
'FIjC  contact  primitive  topologies  shown  in  Fig.  5.37  are  handled  in  a  totally  separate  search  of 
tne  data  structure  and  have  their  own  decision  trees  for  the  sake  of  simplicity  and  debugging. 


Figure  5.36  MHS  form  of  all  possible  topologies  of  noncontact  primitives. 
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Figure  5.37  MHS  form  of  ail  possible  topologies  of  contact  primitir’e. 

5,8  J  Contour  Data  Structure 

The  contour  of  the  remaining  nonprimitives  is  stored  in  the  following  data  structure, 
which  is  visually  represented  in  Fig.  5.38: 

stnict  contour  {Boolean  is_^contact; 
int  portnum; 

struct,  point  firstpoint,  lastpoint; 

} 

struct  point  { int  x,  y; 

short  nodcnum,numprevious,numnexl; 

boolean  comer,  nextq; 
signed  char  fluxdir; 
f 

As  seen  in  Fig.  5.38,  the  total  contour  is  a  liaked  list  of  individual  contovirs,  with  the 
first  contour  being  the  external  boundary  and  tlie  remaining  contours,  if  any,  being  holes  and 
internal  contacts.  Each  contour  is  a  circularly  linked  list  of  points,  with  the  extemai  boundary 
points  stored  clockwise  and  other  points  stored  counicrciockwise.  The  difference  in  orienta¬ 
tion  is  necessary  lor  proper  evaluation  of  the  contour  integrals.  Is  contact  is  TRUE  if  the 
contour  is  a  contact,  firstpoint  and  lastpoint  point  to  the  first  and  last  points  of  the  contour, 
portnum  indicates  how  many  teiminals  the  contour  has. 


Figurs  5.38  Contour  data  structure  for  a  non-primitive. 

X  Jind  y  are  the  coordinates  of  the  point.  node_flum  is  the  node  number  in  the  final  RC 
model  of  the  node  that  the  point  belongs  to.  It  equals  NO_NODE_NUM  if  the  point  belongs 
to  an  insulating  edge.  When  BEM  is  applied,  a  system  of  linear  equations  is  formed  with 
each  point  corresponding  to  an  unknown,  ntimnext  is  the  index  of  the  unknown  belonging  to 
that  point,  For  double  nodes  (Section  5.4.2)  with  two  unknowns,  numnext  is  the  index  for 
and  nurripreYiou.s  of  (Fig.  5.9).  corner  is  TRUE  if  the  point  is  at  a  comer,  nextq  is 


TRUE  if  the  contour  segment  between  the  point  and  the  next  point  is  a  conducting  edge,  i.e., 
q  is  unknown  on  that  segment.  Huxdir  =  1  if  the  outward  normal  n  on  the  segineni  between 
the  point  and  the  next  point  is  or  j\,;  otherwise,  fluxdir  =  -1.  next  points  to  the  next  point 
on  the  contour.  The  contour  of  a  group  of  abutting  rectangles  is  produced  by  the  following 
algorithm  from  [60]: 

1.  Place  vertical  edges  of  all  rectangles  in  a  list  (an  edge  is  two  points  with  the  first  point 
pointing  to  the  second  point  and  the  second  point  being  clockwise  on  the  rectangle  to  the 
first). 

2.  Sort  points  in  the  list  according  to  increasing  y-coordinate. 

3.  Sort  groups  of  points  with  same  y-coordinate  according  to  increasing  x-coordinate. 

4.  Link  point  2k  -  1  with  point  for  1  <  /  <  N/2  (N  is  the  number  of  points)  to  form  the 
horizontal  edges. 

This  algorithm  is  0{NlogN)  since  sorting  is  a.symptotically  the  most  expensive  step. 

5.8.6  FEM  Analysis 

To  do  a  rough  finite  element  analysis  on  the  nonprimitives,  the  abutting  group  of  rectan¬ 
gles  in  the  contour  extraction  method  of  is  removed  from  the  main  list  and  placed  ii  to  their 
own  list.  Then  they  are  broken  into  smaller  rectangles  according  to  steps  1  and  2  of  the  MVS 
procedure.  After  each  comer  is  numbered,  the  FEM  method  of  Section  5.5  is  used  to  form  an 
RC  model,  which  is  then  reduced  using  the  Elmore  time  constant-preserving  node  reduction 
technique  of  Section  5.6. 
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5.8.7  BEM  Algorithm 


Given  below  is  the  BEM  algorithm,  executed  after  the  contour  formation  step: 


Procedure  BEM 

if  number  of  terminals  P  >2 

{multiply  all  coordinates  in  the  contour  by  3; 
divide  each  contour  segment  into  3  segments; 
clear  matrix  A  (N  x  N)  and  source  matrix  B  (N  x  P); 
for  each  delta  source  p  at  a  node  on  the  contour 
for  each  observation  point  q  at  a  node  on  the  contour 
{ form  matrix  A  coefficient; 
form  contribution  to  source  vectors  on  left; 

} 

solve  matrix  equation  AX=B; 

determine  and  output  conductances  from  matrix  X; 

} 

The  contour  coordinates  are  multiplied  by  3;  since  they  are  stored  as  integers,  dividing 
boundary  segment  into  three  pieces  must  yield  integer  coordinates  also. 
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6.  SUMMARY  AND  CONCLUSIONS 


In  this  report,  we  have  described  our  progress  on  the  development  of  current  estimation 
techniques  for  reliability  analysis  of  VLSI  circuits. 

In  Chapter  2,  we  have  described  our  progress  in  developing  and  implementing  the  proba¬ 
bilistic  simulation  approach  for  estimating  average  and  variance  current  waveform  drawn  by 
digital  CMOS  gates  at  contact  point  to  the  power  bus  for  electromigration  analysis  and  for 
average  power  estimation.  Improvements  in  both  computational  speed  and  accuracy  have 
been  made.  In  particular,  a  new  method  of  calculating  the  stochastic  properties  of  the 
equivalent  edge  of  a  channel-connected  subcircuit  during  elimination  has  been  derived,  which 
improves  the  accuracy  of  the  estimation.  The  calculated  expected  value  for  all  the  examples 
we  simulated  was  within  3%  of  the  expected  value  calculated  through  exhaustive  SPICE-like 
simulation,  while  our  previous  method  was  at  times  50%  off.  For  the  variance  our  method 
was  within  12%  of  the  exhaustively  calculated  value,  while  the  previous  method  was  an  order 
of  magnitude  away  in  the  worst  case.  The  calculation  of  the  stochastic  quantities  of  the 
equivalent  edge  during  switching  is  performed  in  a  different  fashion.  The  new  method  is 
straightforward  and,  as  the  examples  show,  more  accurate  than  the  previous  one.  In  addition, 
the  calculation  of  the  delay  through  a  gate  is  treated  from  a  different  perspective,  producing 
more  accurate  results  compared  to  the  previous  method.  Finally,  with  the  introduction  of  the 
correlation  coefficient  method  into  probabilistic  simulation  the  important  issue  of  signal 
dependencies  is  addressed  and,  although  some  error  remains,  execution  time  and  memory 
usage  are  kept  within  reasonable  bounds,  which  permits  the  simulation  of  large  circuits,  The 
approach  has  been  implemented  in  program  iProbe-c.  We  are  currently  working  on  extending 
the  approach  to  include  current  estimation  in  sequential  circuits,  gates  with  pass  transistors. 
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and  mixed-signal  designs. 


In  Chapter  3,  we  have  described  the  application  of  probabilistic  simulation  techniques  for 
estimating  hot-carrier  induced  degradation  within  the  transistors  of  VLSI  CMOS  digital  cir¬ 
cuits.  Since  HCE  in  a  transistor  depends  on  the  current  flowing  through  it  when  the  transistor 
is  operating  in  the  saturation  region,  which  in  turn,  is  a  function  of  the  input  signal  slev/  rate, 
we  have  extended  the  probabilistic  approach  to  include  signal  slew  rate.  We  have  imple¬ 
mented  the  method  in  another  version  of  iProbe,  called  iProbe-d.  We  have  also  proposed  a 
redesign  strategy,  based  on  the  simulation  results,  that  calls  for  interchanging  signal  connec¬ 
tions,  depending  on  the  average  switching  frequencies  of  the  signals,  at  inputs  to  transistors 
connected  in  series  to  reduce  overall  HCE  degradation.  Since  HCE  degradation  affects  the 
timing  of  a  circuit,  we  have  also  combined  HCE  degrrdation  estimation  with  critical  path 
analysis  The  aim  is  to  analyze  HCE  caused  degradation  effects  on  timing  along  critical  paths 
in  the  design  and  recommend  redesign  strategies  to  op  imize  long  term  circt  it  performance. 
More  work  needs  to  be  done  in  this  area. 

In  Chapter  4,  we  have  described  a  linear  time  alg  xithm  (iMax)  that  computes  maximum 
currents  in  the  supply  lines.  Most  of  the  previous  alg  )nthms  on  maximum  current  estimation 
suffer  from  exponential  complexity  and  are  not  at  equate  for  large  circuits.  Our  approach 
avoids  exponential  complexity  by  adopting  a  patten  independent  approach.  The  results  pro¬ 
duced  by  the  algorithm  are  within  acceptable  bcunds  for  most  circuits.  We  have  also 
presented  a  new  partial  input  enumeration  algorithm  that  partially  resolves  the  signal  correla¬ 
tions  and  further  improves  the  upper  bound  obtained  from  iMax.  The  algorithm  is  based  on 
the  best  first  search  (BPS)  technique  and  represepts  a  i;ood  time-accuracy  trade-off.  'I'he  P.?E 
algorithm  involves  a  search  procedure,  but  this  search  teed  not  be  carried  too  deep  to  obtain 
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good  results  The  algorithm  is  quite  applicable  to  VLSI  circuits,  as  is  demonstrated  by  the 
experimental  results  on  circuits  with  up  to  22,000  gates.  In  our  future  research,  we  plan  to 
extend  the  study  to  include  better  gate  delay  and  current  models  and  to  identify  troublesome 
voltage  drop  sites  in  supply  lines,  using  RC  models,  from  the  maximum  current  estimates. 

In  Chapter  5,  we  explain  the  methods  used  in  the  RC  bus  extractor  JET2.  Unlike  our 
previous  extractor  on  which  it  is  based,  JET2  does  not  require  human  intervention  to  edit  sec¬ 
tions  of  the  layout  to  complete  the  extraction.  It  does  this  by  using  the  BEM  with  a  uniform 
grid  of  three  linear  boundary  elements  per  edge  to  model  the  multiport  conductances  of  a  bus 
segment.  It  u.ses  finite  element  analysis  and  an  Elmore  time  constant-preserving  node  reduc¬ 
tion  technique  to  model  the  distributed  capacitance  of  bus  segments.  For  primitives,  the  lesult- 
ing  lumped  capacitance  at  each  node  as  a  percentage  of  the  total  capacitance  to  ground  is 
stored  in  a  look-up  table.  The  lumped  capacitance  for  nonprimitives  is  done  on-the-fiy  using 
FEM  with  a  coarse  grid.  These  methods  lead  to  more  accurate  models  for  both  the  primitives 
and  nonprimitives. 

JET2  can  still  be  improved.  Most  of  the  runtime  on  large  layouts  is  spent  on  the  BEM 
and  specifically  on  the  computation  of  the  boundary  element  coefficients  rather  than  geometric 
processing  of  the  layout  and  primitive  identification.  This  can  be  reduced  by  either  having  a 
larger  static  table  of  primitives  or  storing  shapes  in  a  dynamic  table  as  they  are  encountered  to 
build  a  set  of  primitives  particular  to  that  layout.  Another  solution  would  be  to  improve  the 
calculation  time  of  the  coefficients.  Numerical  integration  via  Gaussian  quadrature  might  be 
effective  in  trading  some  speed  for  accuracy.  Another  improvement  is  the  handling  of  non- 
Manhattan  geometries,  most  probably  those  with  a  fixed  number  of  slopes.  This  would  require 
major  changes  in  the  data  structure  and  geometric  preprocessing  algorithms.  The  BEM  is 
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however  well-suited  to  nonrectilinear  geometries  and  would  require  minor  changes  in  the 
functions  that  compute  the  linear  system  coefficients.  The  rough  FEM  analysis  however  would 
need  a  full-blown  triangularization  algorithm.  Finally  as  feature  sizes  scale  down,  accurate  3- 
D  capacitance  models  will  be  needed,  especially  as  coupling  capacitance  become  more  impor¬ 
tant.  JET2  has  only  a  crude  3-D  model  for  an  interconnect  that  can  be  modeled  as  a  simple 
resistor.  This  model  does  not  take  coupling  capacitance  into  account. 

In  summary,  this  work  has  focused  mainly  on  analysis  and  modeling  issues,  which  are 
the  first  steps  needed  for  including  physical  reliability  issues  in  large  circuit  design.  Even 
though  more  work  still  needs  to  be  done  on  improving  the  analysis  and  modeling  techniques, 
the  next  step  is  to  integrate  the  results  into  a  design  system  so  that  reliability  measures  can  be 
included  at  all  levels  in  the  design  process,  together  with  other  design  objectives,  such  as  tim¬ 
ing,  area,  and  system  performance. 
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